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ABSTRACT 

Synchronous, permanent magnet, and induction machines 
are modeled using computer programs. The computer programs 
incorporate an optimization algorithm which converges on 
lowest weight, volume, and inefficiency. Machine designs 
for high and low rpms are performed, with a varying number 
of pole-pairs. The machine designs are analyzed to find the 
optimum combination of generator and motor for inclusion in 
a naval ship propulsion system. 

Three ships are used for the system study. a baseline 
mechanical transmission ship, a ship retaining the same sub- 
division as the baseline but with the electric machinery, 
and an electric transmission ship with subdivision and 
machinery box arrangement chosen to benefit from the in- 
herent arrangeability of electric transmissions. 

Two generator/motor combinations are used in the final 
ship analysis. Both employ a 3600 rpm, six-pole synchronous 
generator, which turns at the shaft speed of the prime 
mover. One combination uses a 180 rpm, direct-drive, 16- 
pole synchronous motor, and the other uses an 1800 rpm, 
geared, 8 -pole synchronous motor. Power converters are used 
in both combinations to control motor speed. 

The geared combination in the rearranged ship 
demonstrated the best endurance speed efficiency, reducing 
the endurance fuel load by 18%, while maintaining the maxi- 
mum and sustained speed of the baseline ship. The savings 
in ship volume translated to an additional twenty Tomahawk 
missile cells in the rearranged ship. When the fuel load 
was held at the tonnage of the baseline ship, endurance 
range Increased as much as 25%. 

Permanent magnet machines were not competitive in this 
study due to their high weight and volume, even though their 
individual machine efficiency was the highest of all types. 
Induction machines were not used as propulsion generators 
because of the inherent difficulties in control. The induc- 
tion machine motor candidates were not competitive because 
of of f -design-point inefficiency. 
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Chapter One: Introduction 



The use of electric drive as a propulsion method for 
naval ships brings to the ship design process improved ar- 
rangeability 5uid efficiency, though electric machines may 
increase the weight of the plant. Water-cooled electric 
machines are being studied today for ship transmissions^ ; 
these are smaller and lighter, for the same power output, 
than air-cooled machines. They promise reduced overall 
weight for the ship through more economic prime mover load- 
ing, as less fuel will be needed on board. A review of the 
literature has found no work comparing various types of con- 
ventional motors and the effect of each type on the overall 
ship system when used as a propulsion method. St. John [1] 
showed the effects of a superconducting generator /motor 
transmission on the design of a DD963 destroyer hull. Many 
simulations of electric motors and their transients have 
been done. There have been several papers written on naval 
ship integrated electric propulsion systems2 . Also, much 
effort has been expended in the area of electric motor 
design and optimization. Herein, various kinds of conven- 
tional electric machines are modeled. Those machines were 
used in ship designs to find the sensitivity of the designs 
to their use. 



1. Greene, Mole, Welch, and Seng, “Analysis of a High-Power 
Water-Cooled Electric Propulsion System," SHAME Trans., 

Vol 86, 1978, pp 140-162. 

2. Ames, "Marine AC Generation Systems," LSE Journal, 

Vol 12(1), pp 13-29. 
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1.1. Review of electric drive 



Navy ships operate in almost every salt water location 
in the world, including the Black Sea and the Indian Ocean. 
Regular deployments are made to the Mediterranean Sea, North 
and South Atlantic Ocean, and all areas of the Pacific. 

Ships transit the Suez and Panama Canals, the Saint Lawrence 
Seaway, and operate in the Arctic and Antarctic Oceans. 

These operations are made under greatly varying enviroxunen- 
tal conditions, ranging from the sub-freezing temperatures 
of the high latitudes to the hot, dusty conditions of the 
Middle East. All this variety requires ships (and men) 
capable of sustained and efficient operation under any known 
condition. To that end, naval ships are tremendously com- 
plicated systems. 

It is impractical to equip naval ships for every con- 
tingency, but ship designers try to include as much as pos- 
sible when deciding what systems, equipment, and spares to 
put aboard a ship. Once a ship is designed, or as part of 
the design process, political acceptance of the product and 
its purchase is required. Since ships cost tax dollars, 
there are usually limitations on the size aund complexity of 
the design. Still, the designers try to work within the 
given constraints and produce an acceptable and survivable 
(in both the battle and political senses) design. Usually 
this results in a ship that has very small margins in avail- 
able weight and volume. 3 

Weight (displacement) is a semi-direct measure of ship 
cost. Volume is required to place desired systems aboard a 
ship. Therefore, any design change that results in less 
required weight or volume with no decrease in ship 



3. To have small margins is to be limited in the quantity of 
additional system weight and volume that can be added to the 
ship over its lifetime. If a ship is limited in this 
fashion, the flexibility one has in backfitting new systems 
as the ship ages is significantly decreased. 
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effectiveness, is usually a welcome one. 

There are several different kinds of ship propulsion 
systems now in use. They include steam, nuclear, diesel, 
and gas turbine, and there are two principle drive systems: 
mechanical and electric drive. Below is a crude comparison 
of the various propulsion systems, for the purpose of plac- 
ing the thesis in perspective. 

Steam plants burning coal or oil have been in use for 
over a hundred years. They require large amounts of prime 
ship volume, in the center part of the ship.< Steeun is 
produced in boilers and used to power turbines that rotate 
the shafts and propellers through reduction gears. Steam 
plants are very reliable mechanically, but are not terribly 
economic. The large size of the system components demands 
that the boilers and turbine machinery be placed in the cen- 
ter of the ship. This necessitates long runs of shafting 
from the center to the stern of a ship. Shafts typically 
are 18 “-24*' hollow steel cylinders of two to four inch 
thickness; they are heavy and their required placement and 
length makes valuable volume unavailable for other use. All 
propulsion plants except for electric drive have this ar- 
rangement and shafting disadvantage. Steam plants are used 
in all sizes of ships, from the 3000-ton displacement 
frigates to the 50,000-ton battleships to the 80,000-ton 
aircraft carriers. See Table 1 for a list of ship types and 
principle dimensions. 

Nuclear plants are steam plants with a different heat 
source. They also produce steam to power turbines and suf- 
fer the same volume disadvantages as conventional steam 
plants. They are also very heavy and very costly. Mauinlng 
requirements are more stringent, since personnel levels are 



4. The center part of the ship is the most useful ship 
region to place and arrange systems. Ship designers try to 
keep free as much center ship volume as possible. This al- 
lows much greater flexibility in arranging systems that have 
large objects, such as boilers and turbines. 
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Table 1. Typical principle dimensions for various ship types 



Ship types Total 


shaft horsepower 


Length 


Tonnage 


Patrol Hydrofoil 


18,000 shp 


145 ft. 


240 


Frigate 


37,417 


441 


3880 


Destroyer 


78,555 


510 


6395 


Cruiser (non-nuclear) 


82,462 


556 


8872 


Cruiser (nuclear) 


77,460 


590 


9487 


Carrier (non-nuclear) 


280,000 


1050 


80,000 


Carrier (nuclear) 


270,000 


1100 


73,000 


Battleship 


212,000 


887 


58,000 


Submarine (nuclear, attack) 

15,000 


292 


4640 


Submarine (nuclear, strategic missiles) 

16,000 


410 


7880 



Source: Jane’s Fighting Shins 1985-86 . Jane’s Publishing 

Company, Ltd. , London, edited by John Moore. 

These numbers represent the geometric mean of several 
classes of ships within a type and should not be taken to be 
those of a particular ship. Their value lies in the ap- 
preciation of the differences between various ship types. 

An apt weight comparison would be that a typical forty foot 
sailboat might displace fifteen tons. 

The U. S. Navy has other ship types besides those listed 
above. They include amphibious warfare ships, submarine and 
destroyer tenders, and fleet oilers, and supply ships. 
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rigidly controlled and the operators of the propulsion plant 
must be nuclear trained. Nuclear plants are used on 
cruisers, submarines, and aircraft carriers. 

Diesel propulsion plants have high weight to volume 
ratios, making their use costly in weight dollars, but 
cheaper in volume dollars.® They are very noisy, ruling out 
their use in antisubmarine warfare ships. Although very 
reliable mechanically, they require much tinkering and 
tuning. Their specific fuel consumption® is among the 
lowest of all the plants. Medium speed diesel engines are 
not commonly manufactured in the 25,000 hp range, which 
means diesel plants cannot be used in ships requiring high 
shaft horsepower. They are typically used in smaller ships 
as cruise engines and in amphibious warfare ships as main 
propulsion. (Amphibious warfare ships typically have lower 
top speeds than frigates or destroyers . ) 

Gas turbine propulsion plants seem to have many good 
points. They are reliable, quiet, relatively low weight, 
and come in power ranges that are useful in ships ranging 
from 300 ton hydrofoils to 8000 ton destroyers. They 
require large amounts of volume for intake and exhaust duct- 
ing, but this is not a great disadvantage. Their fuel 
economy is not as good as other plants, but this is not 
intrinsic to the gas turbine engine. It is a fault of the 
operating method; gas turbine plants have mostly been built 
with mechanical transmissions. Usually one or two engines 
are coupled to each shaft. If the ship is proceeding at 
high speed, the gas turbines are operating at their full 



5. When designing a ship, total ship cost is monitored by 
the use of marginal cost factors. Every additional cubic 
foot or ton of weight added has a marginal cost associated 
with it. When a ship’s total cost is constrained, design 
changes that add cost are discouraged or must be offset by 
the reduction of other systems’ weight or volume. 

6. Specific fuel consumption is the ratio of pounds of fuel 

burned per horsepower-hour. SFC=lbf /hp-hr . 
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load design point and are relatively fuel economic. Good 
fuel economy is not usually realized, however. In the main, 
the ship proceeds at a cruise speed, using one gas turbine 
for each shaft, and the gas turbines operate at about half 
power. Specific fuel consumption rises rapidly as gas tur- 
bine power level drops, which makes for inefficient opera- 
tion. 

Usually, the above propulsion plants have mechanical 
transmissions. This means the main engines, whether they 
are diesels, gas turbines or steam turbines, are mechani- 
cally connected to the shafts and propellers. There is 
usually a reduction gear between the engine and shafting. 
These gears are large, very heavy, and expensive. To 
provide an idea of size, the largest, or "bull’' gear in a 
typical locked-train double-reduction gear is about seven 
feet in diameter. The reduction gear must be placed in-line 
with the shafting, thereby using more of that prime ship 
volume. Some mechanical transmissions have cross- 
connections between shafts, but this is not common. 

Electrical transmissions are characterized by prime 
movers of any type providing power to generators. The out- 
put electricity is conditioned and sent to propulsion motors 
via a distribution network. Cross-connection is done with 
switches and breakers. There can be a mechanical reduction 
gear if it is desired to operate the propulsion motors at 
higher than propeller rotational speeds. Direct-drive 
motors may also provide the desired propeller rpm, e.g. , by 
controlling the field current in a DC motor. The propulsion 
motors can be very near the propeller, i.e. aft, eliminating 
the long runs of shafting associated with a mechanical 
transmission. 

Naval ship propellers are of two types, controllable or 
fixed pitch. The pitch of a propeller is the distamce the 
ship moves forward in the water for one turn of the propel- 
ler. A fixed pitch propeller has this characteristic dis- 
tance the scune at all times. A controllable pitch propeller 
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can vary this distance by changing the angle of attack (the 
angle at which the blade slices through the water) of the 
propeller blades, including reversing the blade so that the 
ship moves astern. Controllable pitch propellers are prac- 
tically required for propulsion plants that have non- 
reversing main engines, such as diesel and gas turbine 
plants."^ Steaun plants can reverse their propellers and 
shafts by use of an astern turbine, albeit much more slowly 
than a ship with a controllable pitch propeller system. 

Fixed pitch propellers have a slight advantage in efficiency 
(1-3%) over the controllable ones. This is mostly due to 
the large propeller hub required for varying the blade angle 
of a controllable pitch propeller. Quick reversal of shaft 
direction or propeller pitch means quick ship braking and/or 
ship reversal. This ship quickness is mandatory for an- 
tisubmarine operations and safe navigation. For example, 
the ability to stop "on a dime" may be important in a 
crowded sea lane, where a small wooden sailboat has the 
right of way over a powered naval vessel . 

Electric drive seems to combine the best of all the 
propulsion plants. It has all the advantages of a conven- 
tional prime mover plus the advantage of electrical cross 
connection and better arrangements. In the cruise scenario 
above, the electric drive ship could have both shafts 
operating from one gas turbine engine. That engine would be 
coupled to an electric generator which would produce enough 
power to run the motors that turn each shaft. The drive 
motors would be placed at the stern of the ship, near the 
propellers, on the same level. The long runs of shafting 
would be replaced by electric cable, which is smaller, 
weighs less, and can be placed in non-prime real estate. 



7. Ships with non-reversing prime movers can also have a re- 
versible reduction gear with a fixed pitch propeller instead 
of a controllable reversible pitch propeller. This is new 
technology for the United States aind only the latest naval 
ship design, the DDG51, has a reversible reduction gear. 
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Cable is also in many cases cheaper than shafting, espe- 
cially to repair. Since one gas turbine would provide power 
to both shafts, it would operate at a higher power level and 
would be therefore more fuel economic. A typical propulsion 
plant might consist of three gas turbines with three gener- 
ators. Most ships have an even number of prime movers be- 
cause mechanical shaft cross connection between shafts is 
not often used and each shaft in a mechanical transmission 
ship requires the same number of prime movers to balance 
loading at high power levels. The extra prime mover 
requires a lot of weight and volume. An advantage of 
electric drive is that it becomes possible and perhaps 
desirable to use an odd number of prime movers . Each of the 
two shafts would have one propulsion motor. The heavy 
reduction gears could be replaced by the motors, which would 
have an infinitely variable reduction ratio. The control- 
lable pitch propeller system so far required by this gas 
turbine ship would be replaced by the cheaper, slightly 
smaller, and far less complicated fixed pitch propeller, ad- 
ding a small efficiency gain. The hydraulic system used to 
vary blade angle would be eliminated. Ship braking and 
reversal would be accomplished by electrically controlling 
the motor rotation direction, combined with energy dissipa- 
tion through the use of resistor banks. 

A disadvantage of this arrangement would be the high 
weight of the propulsion motors. They would be special 
designs and have a high capital cost. Hopefully, the high 
weight of the motors would be offset by the reduction in 
shafting and fuel weight and the possible elimination of the 
reduction gears. The high cost would be made palatable by 
the savings in fuel over the life of the ship. A Life Cycle 
Cost comparison of various propulsion plants, including 
electric drive, is available in reference two. 

The change to electric drive would likely be accom- 
panied by an overall decrease in propulsion plant weight and 
volume. The ship could be smaller and lighter, and would 
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require less onboard fuel for the same endurance range (the 
distance the ship can travel without refueling) . Since less 
volume would be required for the fuel, the ship could be 
smaller and lighter. Since the ship would then be smaller 
and lighter, less horsepower would be required to achieve 
the same top speed. Since less horsepower would be 
required, the ship could be smaller and lighter. This is an 
example of the design spiral that would result in a smaller, 
lighter, cheaper, more risk-free ship. An example of this 
type of ship improvement is given in reference one. There 
is a limit on ship improvement, usually due to the non- 
propulsion systems or payload. One cannot make an ocean 
crossing missile ship the size of a small yacht. 

So why are not all Navy ships electric drive? They are 
not largely because the technology has not existed in a 
usable, fully developed, and manageable form. Because of 
the high cost of naval ships (a small one may cost $350 
million) and the lack of experience with current electric 
drive technology, the Defense Department is reluctant to 
build large electric drive ships. There have been electric 
drive ships, including five battleships with 21 MW shaft 
output and two aircraft carriers with 135 MW shaft output. 
Over 160 escort vessels were built during the Second World 
War with turboelectric or diesel-electric drives ranging 
from about 4.5 to 9 MW.® A new class of ocean surveillance 
ship, the T-AGOS 19, is being built with diesel-electric 
propulsion, but it is only a 3500 shaft horsepower (shp) 
ship. 

Electric drive was replaced by conventional mechanical 
transmission plants after World War Two because of the com- 
petition afforded by Improved gear cutting methods. Double- 
reduction locked-train gear transmissions became the 
standard. Since the electric drive ships all had non- 



8. Doyle, T. J. and Harrison, J. H. , "Navy Superconductive 
Machinery Program," Trans. SNAME, 1978, p. 20-1. 
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superconducting, air-cooled motors and generators, they had 
higher weight and increased space requirements and suffered 
in comparison with the mechanical transmission ships. 9 The 
importance of the improvement in power electronics must be 
mentioned. World War Two ships did not have the advantages 
afforded by those electronics. 

Integrated electric drive propulsion must also be men- 
tioned. This propulsion plant is the same as any of those 
discussed above, except that ship service electrical power 
is derived from the main propulsion plant, usually by taking 
power off the reduction gears or main propulsion generator. 
Power conditioning equipment, such as a cycloconverter, is 
needed to "clean up" the power and change it to fixed 
frequency for use in other equipment. Variable speed con- 
stant frequency equipment and concepts embody the integrated 
electric drive concept. The U. S. Navy has investigated 
this in some detail. lo 

Some requirements of electric drive may be viewed as 
disadvantages. The power from the electric generators has 
to be conditioned to provide frequency control of the 
propulsion motors. The power conditioners add weight and 
volume to the overall system, as well as reducing the 
efficiency of the transmission. Braking resistors, used to 
dyneunically and quickly slow the propulsion motors, add more 
weight and volume to the system. There also may be a high- 
frequency radiated noise signature associated with alternat- 
ing current systems that may be deleterious to the mission 
of the ship. 

The research done to date has not explored specific 
motor types in detail. How can it be decided whether to put 
a synchronous , inductive , permanent magnet or other motor in 



9. Ibid. 

10. Robey, Stevens, and Page, "Application of Variable Speed 
Constant Frequency Generators to Propulsion-Derived Ship 
Service," Naval Engineers Journal, May 1985. 
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the electric drive system? What makes the ship system 
"best"? The effect of each motor type on the ship system 
has not been analyzed. The electric transmissions used in 
the current version of the Advanced Surface Ship Evaluation 
Tool (ASSET), a ship design computer program written for the 
United States Navy by Boeing Computer Systems, Inc. , are 
generic combinations of AC and DC motors and generators, 
using rough estimates of weight, volume, and efficiency. 

Ship designs more involved than the feasibility level need 
detail on just those items. 

Motor design is a well known subject and there are many 
texts on the subject. The use of motors and pertinent tech- 
nologies in a ship as a propulsion method is discussed in 
Greene, Powell, and Gripp [3]. The advantages of electric 
drive include flexibility of arrangement, controllability, 
variable reduction ratio, reliability, and provision of ship 
service power from the main bus. Jolliff and Greene [4] go 
on to propose a specific water-cooled Advanced Integrated 
Electric Propulsion Plant (AIEPP) for a frigate/destroyer- 
sized ship. They discuss the essential characteristics of 
such a plant, establish the feasibility of the drive system 
and identify the method to technically demonstrate the sys- 
tem. Acker, Greene, and Jolliff [5] present several model- 
ing techniques and scaling relationships that allow estima- 
tion of volumes and weights of propulsion motors and gener- 
ators, solid state power conditioners, electrical 
switchgear, and associated electrical propulsion systems 
components as functions of propeller shaft power. A case 
study of AIEPP is given in the paper by Kastner, Davidson, 
and Hills [6]. 

Simulation of electric motors and associated systems is 
a popular topic. Many persons have done work in this area, 
from the micro-consideration of high frequency inductance 
changes to the more macro-consideration of hunting tran- 
sients, etc. Smith, Stronach, and Tsao [7] model a complete 
electromechanical marine drive system while Smith, Stronach, 
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Tsao, and Goodman [8] concern themselves more with a marine 
power system, including pump drives. Nonlinearities and 
operational transients are addressed. 

1.2. Optimization 

Optimization is the process of making a system, sub- 
system, or idea the best it can be. "Best" is defined by an 
"objective function," a measure of what is optimum. For ex- 
ample, an optimal manufacturing process may produce the max- 
imum number of units at the lowest cost. The objective 
function would combine units-produced with cost in an equa- 
tion that could be analyzed to find the proper production 
level. The output of the objective function is a single 
scalar measure of "goodness." It may be difficult to repre- 
sent complicated processes with only one number. 

Optimization can be performed on a global or subor- 
dinate basis. The optimum motor might be the one that has 
the highest efficiency, even if that efficiency was achieved 
by designing a very large and heavy motor. The sub-system 
(the motor) has efficiency as its objective function. The 
ship in which the motor is to be placed may be optimum when 
its overall weight and volume are the lowest (ignoring cost, 
for example). A large, heavy motor, then, may not be op- 
timum for the ship, even if it is very efficient. A good 
case study of motor optimization is the EPRI report authored 
by Fuchs, et al . [9] 

"Optimization" can be an ill-defined term but there are 
fairly well defined methods of achieving it. Linear 
programming, Markov modeling, and Monte Carlo schemes are 
examples of these methods. The accessibility of high-speed 
digital computers has made multiple random excursions in a 
multi-variable space a much easier way of finding the 
"optimum" solution, provided an objective function and con- 
straints can be devised to describe the problem. This 
method of random excursions (Monte Carlo scheme) and ex- 
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amples of it are among the methods described in references 
ten, eleven, and twelve. 

Monte Carlo schemes take their name from the action of 
the roulette wheel in the gambling casinos of Monte Carlo. 
Around and around the wheel goes, stopping on random num- 
bers. If the computing power is available, this is an ac- 
ceptable method of exploring a large variable space. It can 
be much quicker than looking at every possible permutation 
of all variables. 

The steepest-descent scheme is so named because op- 
timization moves down the sharpest gradient of the objective 
function. From a valid design point, random steps are taken 
in every variable and the design point is moved "downhill" 
toward the objective function over the steepest slope. This 
is different from the "drunkard’s walk," where the random 
steps are only evaluated on whether or not the objective 
function’s output has improved, not if it is improving at 
the fastest possible rate. 

Optimization Is almost always subject to constraints. 

In the previous manufacturing process, warehouse space may 
be limited, so only a certain number of units may be stored. 
This could act to limit production. For motor design, con- 
straints include maximum rotor tip speed, maximum current 
density, minimum power output, etc. All constraints should 
be combined with the objective function to yield a 
"constrained optimization." 

1.3. The objective function 

Optimization is not possible without an objective func- 
tion. It may be very difficult to devise a good objective 
function for a complicated system such as a ship. It may be 
even harder to find one for a sub-system of that complicated 
system. There are very many characteristics that could be 
optimizing variables, and assembling them into one objective 
function with all the constraints is not easy. 
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Even with a properly defined objective function, it may 
be difficult to choose among designs that result from the 
constrained optimization. For example, a low-volume and 
low-weight motor may have poor efficiency. A very efficient 
motor may also be large and heavy. If all three elements 
are important, which is the best motor? Deciding between 
competing designs has been the subject of various papers, 
one of which is by Schweppe and Merrill [13]. In that 
paper, the authors suggest the use of "knee curves,” saying 
that the essential characteristics of a multiple attribute 
tradeoff can be plotted on a series of x-y graphs. Uncom- 
petitive designs are easily discerned and discarded. The 
decision process can be limited to only those designs that 
are competitive. 

Table 2. Optimizing characteristics for ship and transmis- 
sion 



Weight 

Volume 

Efficiency 

Cost 

Reliability 

Maintainability 

Commonality 

Manning 



The above table lists many of the possible optimizing 
characteristics for the ship and its propulsion sub-system. 
Manning estimates are typically based on historical data and 
do not indicate that the baseline and variant ships will 
require a significantly different number of personnel. Com- 
monality measures the use of the same equipment in other 
ships . Since there are no other electric drive ships at the 
power levels used in this thesis, commonality is not an 
issue. It is very difficult to discern maintainability and 
reliability differences between designs that are as close as 
the machine designs of this thesis, so these two charac- 
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teristics were not made part of the objective function. 

Cost is a measure that should be part of this sort of op- 
timization. The cost of the ship system, quantified in the 
ship displacement, was used in the final recommendations for 
the transmission sub-system. In the case of permanent mag- 
net machines, the relative costs of magnet material and mag- 
net steel were included in the objective function. 

Weight, volume, and efficiency were made an explicit 
part of the objective function for the computer design of 
the machines . An Effective Weight was calculated for every 
machine design. The design with the lowest effective weight 
was the "best" within its class. The generic objective 
function is 

Effective Weight = weight + ke*(l-ef f iciency ) + kv*volume 

where ke and kv are weighting factors for efficiency and 
volume, respectively. The weighting factors were obtained 
from changes in ship displacement for marginal changes in 
efficiency and volume of the transmission. They were 
modified to reflect the actual designs resulting from the 
process . 
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Chapter Two. General Considerations 



Only steady-state behavior of electric machines was 
modeled. The modeling of dynamic behavior is very difficult 
and was not viewed as being within the design problems posed 
by this thesis. The changes in machine design necessary to 
solve dynamic instabilities, etc., are much smaller than the 
approximate nature of the algorithms used here. 

All derivation work was performed without specifying 
the number of winding turns or the number of rotor or stator 
slots. The only exception to this was the case of induction 
machines, where an arbitrary number of rotor slots was 
selected. This selection was necessary for the calculation 
of the equivalent circuit components. The number of rotor 
slots chosen, 71, was a number designed not to Induce pole 
harmonics. Since no turn numbers were specified, units in- 
clude volts-per-turn, ampere-turns, and impedance-per-turns- 
squared. Power is measured in watts. 

All machines used as their synchronous frequency the 
maximum allowed by the particular combination of pole pairs 
and shaft rpm. Developmental work showed that the optimiza- 
tion algorithm converged to the highest frequencies, so the 
algorithm now starts at the highest possible frequencies. 

Up to six pole pairs were used in the higher rpm 
machines and up to 25 pole pairs in the 180 rpm machines. 
Diminishing improvements in volume, weight, and efficiency 
show up at half these limits. 

The random number generator was taken from Kelley and 
Pohl [14], with one change. After every run of each 
program, the random number generator seed is stored. This 
means the sequence of pseudo random numbers is not repeated 
until the full range possible has been used. For the 
machines, it gives differences at every run and means the 
multidimensional variable space is more fully explored. 
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2.1. Optimization method 



The chosen optimization method is a combination of the 
Monte Carlo and steepest-descent schemes. A design point is 
established by randomly selecting machine geometric 
parameters, subject to constraints. Ten random steps are 
taken around the design point, in all variable directions. 
The effective weight of each random step is evaluated and 
compared with that of the design point. The best of the 
eleven is designated the new design point. More random 
steps are taken, and the process continues until no more im- 
provement is seen in effective weight. At that point, the 
size of the random steps is halved, and the process repeats 
itself, with the step size continually halved (up to ten 
times). The best effective weight is the index to the best 
design. The number of original design points used in any 
particular run of a program is under user control. If ten 
original design points are desired, the algorithm will look 
at over a thousand designs. 

The purpose of having original design points is to 
start the optimization process in different sectors of the 
multidimensional variable space. In this fashion, the op- 
timization process zeroes in on several local "best" points. 

The variables that are randomly selected include stator 
current density, rotor radius, air gap dimension, stator 
slot space factor, and rotor slot space factor. The back 
iron dimension (the iron behind the stator teeth) is sized 
to handle a saturation level of flux. The stator slot depth 
is originally sized as a random fraction of the back iron 
dimension, and the rotor slot depth is originally a random 
fraction of the rotor radius. 

Only steady-state behavior was modeled in this algo- 
rithm. Dynamic modeling may or may not show different op- 
timum configurations for machine types. 
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2.2. Constraints 



The constraints placed on the optimization process are 
listed in the following table. The most difficult con- 
straint to satisfy while still achieving a valid design was 
the rotor current density constraint in induction machines. 
Only a few valid designs were achieved in induction machines 
using the above algorithm, leaving some question about the 
application of the algorithm in the case of induction 
machines. Only those induction machine designs in which 
there was reasonable confidence were included in the thesis 
analysis . 



Table 3. Optimization constraints during machine design 



Minimum air gap flux density 
Maximum (saturation) flux density 
Maximum rotor radius 
Maximum rotor tip speed 
Stator and rotor space factor 
Maximum rotor slot depth 
Maximum synchronous reactance: 
synchronous machines 
permanent magnet machines 
Power factor 



1.05 tesla rms 

1.5 tesla rms 
2 . 0 meters 

200.0 meters/sec 
0.35 

33% of rotor radius 

2.0 per unit 

3 . 0 per unit 
0.8 



The magnet steel chosen was 26 gauge M19. Its magnetic 
properties were found in USX technical data [15]. It has 
been observed that saturation flux levels in electric 
machines occur first where the area perpendicular to the 
flux path is the smallest. If the back iron dimension is 
made appropriately large, this saturation will first occur 
in the teeth, as is desirable. Accordingly, the back iron 
dimension was set to 

Br r 

dcore = 

Bsat p 

where Br is the radial air gap flux density, r is the rotor 
radius, Bsat is the saturation flux density, and p is the 
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number of pole pairs in the machine. This equation is 
derived from Gauss’ Law. 

If electric machines are to be installed in a ship, th 
78 ey obviously will need to fit into the space designated 
for them. An electric machine with a two meter rotor radius 
is at least thirteen feet in envelope diameter. This is a 
very large machine to install in a machinery space where 
volume is already at a premium. The rotor radius limit de- 
scends from the physical ability to fit an electric machine 
in a ship. 

The tip speed limitation represents the physical limit 
on material strength with regard to the rotor conductors. 
Rotor conductors may break free from the rotor at higher 
tangential velocities than this limit. The magnitude of the 
limit was taken from a tip-speed-limited, 3600 rpm, two-pole 
turbogenerator, and was verified against standard Navy 
design practices. Several as-built electric machines were 
analyzed and this number seemed to fit their characteristics 
well. The tip speed limit arises when choosing a rotor 
radius (given a particular frequency and number of poles), 
and is less stringent a constraint than maximum rotor 
radius . 

Thermal considerations are often extremely important in 
machine design. The heat build-up in electric machines has 
led to many cooling schemes over the years, including 
natural convection, forced air cooling, and hydrogen cool- 
ing. One of the latest methods is liquid cooling of the 
stator and rotor conductors through cooling passages through 
the copper itself. This has been made possible by better 
de-ionizing methods for cooling fluid and better rotating 
seals for the rotor. Naturally, the cooling passages and 
insulation limit the amount of copper area in a slot cross 
section. The copper area in a typical conductor bar was 
measured and found to be about thirty-five percent of the 
bar cross section. This number was used for the stator and 
rotor slot space factors. 
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Other thermal considerations must be made for permanent 
magnet materials, which suffer a degradation in flux as tem- 
peratures rise. Flux loss rises slowly with increasing tem- 
perature until about 100* C. Above 100® C, flux loss is more 
rapid. An assumption in the design of these machines is 
that there will be sufficient cooling in the operating space 
to limit ambient temperature to about 80® C. This, combined 
with the machine liquid cooling, should keep flux loss to a 
minimum. Transmission lines were asstimed to function satis- 
factorily at the same temperature. 

Insulation also has a thermal rating. No insulation 
class was specified in this thesis but a typical insulation 
used in electric generators by the Navy is Class F. For 
this class, a permissible rise of 100® C over an ambient tem- 
perature of 50® C is standard, but lesser insulation classes 
must run cooler. If a machine must be designed with a les- 
ser insulation class, the consequent lower temperatures may 
result in a quieter machine and longer machine and insula- 
tion life. It probably will be larger than a machine with a 
greater class of insulation. The temperatures quoted above 
are at hot spots , not in the bulk of the machine . 1 1 

Along with the reduction in copper area for liquid 
cooling, a maximum current density was Imposed. The copper 
losses, in the form of heat, have to be removed by the cool- 
ing fluid. There is a tradeoff between the size of the 
cooling passages, the allowable current density, and the 
rating of the machine. Twelve million amperes per square 
meter equates to forty amperes in a twelve gauge copper 
wire . 

Rotor slots were constrained to no more than one third 
of the rotor radius. Some reasonable shaft diauneter is 
required to transmit the mechanical torque. Stator slots 
were allowed to grow as needed to meet the stator current 



11. Private communication, D.F. Schmucker, Naval Sea Systems 
Command, Code 56Z31, 11 December 1986. 
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density limit. 

Synchronous reactance limits were taken from as-built 
machines. A power factor of 0.8 was used for all machines, 
though one researcher has indicated a power factor of 1.0 
might be best for permanent magnet machines. 12 

2.3. Geometric considerations 

End turns were modeled as described in Appendix B. 

While not exactly as machines are constructed, this model 
gives reasonable results. A length allowance equal to one 
rotor diameter on either end of the active length of the 
machine was made to allow for containment of the end turns. 

Fractional slot pitches were not considered. A stator 
winding pitch of 0.8 was assumed, resulting in the elimina- 
tion of the fifth harmonic from the steady-state output 
waveforms . 

Once weight and volume were calculated, an extra ten 
percent was added to allow for the frame and foundation of 
the machine. The calculated weight included an additional 
three percent of the rotor weight to allow for bearings. It 
is the final envelope weight and volume that were used in 
the decision process. 

2.4. Efficiency and losses 

The general equation for efficiency is 

minpwr 

efficiency = 

minpwr + ph + pe + i2r + i2rr 

where minpwr is the minimum mechanical power expected of the 
machine, ph is the hysteresis loss, pe is the eddy current 



12. Robey, H.N., "Permanent Magnet Machine Technology 
Assessment," DTNSRDC Report TM-27-80-87, September 1980. 



29 



loss, i2r is the stator copper loss, and i2rr is the rotor 
copper loss. This formulation is for a motor, but the ef- 
ficiency calculated will not be significantly different if 
the machine is a generator. 

Hysteresis and eddy current losses arise from currents 
circulating within the magnet steel that forms the rotor and 
stator. They are two different mechanisms and depend on the 
metallurgy of the steel. 

Eddy currents are a result of the time-varying magnetic 
fields within the machine, and they oppose the change in 
flux density within the machine. Eddy current losses in- 
crease as the square of the electrical frequency of the 
machine and also as the square of the peak flux density. 

One method of lessening eddy current losses is to use thin 
laiminations to build up the rotor and stator. If the var- 
nish used on the laminations is sufficiently insulating, the 
eddy currents are limited to azimuthal circulation. Axial 
circulation is practically zero because of the small lamina- 
tion thickness. 

Hysteresis losses are inherent to magnetic materials, 
and are proportional to the total volume of the material, 
the area of the hysteresis loop, and the machine electrical 
frequency. 

USX has developed equations to calculate eddy current 
and hysteresis losses. They are 



ph 






and pe 



0^4818_n_Bm* t^ f^ 
rho D 



In these equations, (i is the hysteresis loss factor (the 
ratio of the actual hysteresis losses to the area of a 
square hysteresis loop passing through Br and He), f is the 
frequency in Hertz, Br is the residual induction in 
kilogauss. He is the coercive force in oersteds, D is the 
density in grauns per cubic centimeters, rho is the electri- 
cal resistivity of M19 in microhm-cm, fl is the anomalous 
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loss factor. The losses are in watts per pound of material. 
The numerical factors at the beginning of both equations 
were altered to reflect the use of SI units. The factors 
n, Br, He, rho, and D change with the type of magnetic 
material used. The thickness of the laminations, t, is 
0.014 inches for 26 gauge steel. These equations do not 
reflect some variations caused by differences in silicon 
content and differences in processing treatments leading to 
variations in grain size and crystallographic texture. 
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Chapter Three . Synchronous machines 



Synchronous machines operate because of an interaction 
between stator and rotor flux waves. The rotor flux wave is 
developed by a field winding. The stator flux wave is 
developed by an armature winding. These two waves try to 
align themselves, which is how the machine action is 
produced. In the case of a motor, the armature wave is 
"rotating" around the periphery of the stator bore because 
of the 120® separation between the three phases. A rotating 
action ensues. In the case of a generator, the rotation is 
provided by a prime mover, such as a gas turbine, and volt- 
age is induced in the stator phases. 

Synchronous machines operate at a steady-state shaft 
speed specified by the number of poles and the electrical, 
or synchronous, frequency. This synchronous speed is main- 
tained despite changes in load. This feature makes 
synchronous machines attractive for applications where speed 
control is important. Shaft rpm = ( 120* frequency )/(poles ) . 

A derivation of the equations of synchronous machines 
and the computer modeling program are presented in 
Appendix B. 

3.1. Assumptions 

The rotors of synchronous machines may exhibit 
saliency, or may be smooth cylindrical rotors. The dif- 
ferences in properties and parameters among salient and 
round-rotor machines amount to only a few percent! 8 . The 
approximate nature of the modeling means the saliency ef- 
fects will not be important. Therefore, no special provi- 
sion for salient rotors were made. 

The ships used in this thesis have displacements of 



13. Fitzgerald, Kingsley, and Umans, Electric Machinery . 
McGraw-Hill, New York, 1983. 
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about 5000 LT ( 1 LT = 2240 lbs). Their fuel load is deter- 
mined by the required range and transmission efficiency at 
endurance speed. It was first estimated that a one-percent 
increase in transmission efficiency would reduce total ship 
displacement by about 89 LT. The efficiency factor (ke = 
90,000 kg/percent) corresponds to this 89 LT change in ship 
full load displacement. (89 LT x 2240 Ibs/LT x 2.205 kg/lb) 
When this produced machines with efficiencies about 95%, it 
was doubled to 180,000. Obviously, this factor may be ad- 
justed to any level. The volume efficiency factor (kv = 
1286.1 kg/m3 ) corresponds to the density of a LM-2500 gas 
turbine module. These factors were used throughout the 
thesis . 

3.2. Machine description 

Machines with shaft speeds of 1800, 2400, 3000, 3600, 
and 7200 rpm, with the number of pole pairs varying from one 
to six, were modeled. Also, 180 rpm machines using from one 
to twenty-five pole pairs were modeled. This provided a 
good coverage of the variable space. 

3.2.1. Efficiency 

Synchronous machine efficiency at full load was about 
98.5% for the higher rpm machines, while the 180 rpm 
machines hovered around 93% to 94% efficiency. The number 
of pole-pairs seemed to have little effect in the high rpm 
machines, but there was an "arch" in the efficiency curve of 
the 180 rpm machines, peaking at 97% with 36 poles. Though 
not fully understood, the 24- and 26-pole 180 rpm machines 
had very low efficiencies. Generally speaking, machine ef- 
ficiency was higher when rpm was higher. This was an ex- 
pected result. Of f -design-point efficiency was good for the 
higher rpm machines, but bad for the 180 rpm machines. (See 
the discussion in Chapter Seven. ) 
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3.2.2. Weight and volume 



Weight and volume increased almost linearly with the 
number of pole-pairs in the 180 rpm machine. This is almost 
certainly a function of tip-speed limitations, as the maxi- 
mum tip speed in a machine is a function of the number of 
poles in the machine. When the rotor radius is limited, the 
machine must grow in length to develop enough torque. The 
weight and volume of these machines were much higher than 
for the higher rpm machines . 

The higher rpm machines saw significant decreases in 
weight and volume when the number of pole-pairs increased 
from one to two. There also was an observable increase in 
weight and volume as the number of pole-pairs further in- 
creased. As rpm increased, the machines became smaller and 
lighter. 
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EFFICIENCr 



Figure 1. High rpm synchronous machine efficiency 



SYNCHRONOUS MOTOR EFFICIENCY 




35 



EFFICIENCY 



Figure 2. 180 rpm synchronous machine efficiency 



SYNCHRONOUS MOTOR EFFICIENCY 



DIRECT-DRIVE, 180 RPM. 25775 HP 
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VOLUME (cubfc motors) 



Figure 3. High rpm synchronous machine volume 



SYNCHRONOUS MACHINE VOLUME 
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VOLUME (M^3) 



Figure 4. 180 rpm synchronous machine volume 



SYNCHRONOUS MOTOR VOLUME 



DIRECT-DRIVE, 180 RPM, 25775 HP 
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WEIGHT 

(Thousands) 



Figure 5. High rpm synchronous machine weight 



SYNCHRONOUS MOTOR WEIGHT 



25775 HP 
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WEIGHT (KG) 
(Thousands) 



Figure 6. 180 rpm synchronous machine weight 



SYNCHRONOUS MOTOR WEIGHT 



DIRECT-DRIVE, 180 RPM, 25775 HP 




40 



3.3. Machine characteristics 



The following tables give machine characteristics for 
many of the higher rpm and the 180 rpm synchronous machines 
designed for this thesis. 

The stator slot factor tended to increase to the limit 
of 0.75, while the rotor slot factor moved around the value 
0.58 quite a bit. This demonstrates a partial limit on the 
depth of the stator slots (to the same dimension as the back 
iron depth). The overall diameter was limited and more 
stator slot area was needed to develop the required power. 
Similarly, the stator current density converged to its coo.^-: 
ing limit. 

The longest higher rpm machine was 5.42 meters, while 
the largest overall diameter was 0.85 meters. Machines of 
this size will cause no difficulties when placed in the 
machinery spaces of most ships. The 180 rpm machines are 
larger, typically less than 6 meters long (discounting one 
16 meter machine) and 1.8 meters in diameter. They are also 
good candidates for ship systems. 
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Table 4. Characteristics of 1800 rpm synchronous machines 



nuiber of pole pairs 
poNer, hp 
efficiency factor 
voluie factor 
shaft rpi 
stator current density 
synchronous frequency 
rotor radius 
gap diaension 
back iron depth 
stator slot depth 
rotor slot depth 
stator slot factor 
rotor slot factor 
envelope volute 
envelope weight 
hysteresis loss 
eddy current loss 
stator copper loss 
full load efficiency 
active length 
full load current density 
no load current density 
xs/turns'squared, p.u. 
internal volts/turn, p.u. 
overall length 
overall diaaeter 



1 


2 


3 


25775 


25775 


25775 


180000 


180000 


180000 


1286. 1 


1286.1 


1286.1 


1800 


1800 


1800 


8.00E+06 


1.20E+07 


1.20E+07 


30 


60 


90 


0.1720 


0.2077 


0.2603 


0.0158 


0.0440 


0.0250 


0.1204 


0.0727 


0.0607 


0.0772 


0.0722 


0.0607 


0.0401 


0.1859 


0.0747 


0.653 


0.750 


0.750 


0.650 


0.447 


0.687 


2.781 


1.522 


1.509 


20835.40 


9827.96 


8766.09 


16926.49 


10677.39 


17111.27 


3033.2 


3826.7 


9198.8 


275834.4 


374646.5 


293335.5 


0.985 


0.980 


0.984 


4.670 


1.794 


1.498 


1.19E+07 


1.06E+07 


1.17E+07 


4.40E+06 


6.39E+06 


7.00E+06 


1.99 


0.86 


0.87 


2.71 


1.66 


1.67 


5.42 


2.80 


2.64 


0.77 


0.79 


0.81 



4 


5 


6 


25775 


25775 


25775 


180000 


180000 


180000 


1286.1 


1286.1 


1286.1 


1800 


1800 


1800 


1.20E+07 


1.20E+07 


1.20E+07 


120 


150 


180 


0.2905 


0.3267 


0.3189 


0.0124 


0.0090 


0.0092 


0.0508 


0.0457 


0.0372 


0.0508 


0.0457 


0.0372 


0.1333 


0.0907 


0.0289 


0.749 


0.750 


0.748 


0.659 


0.449 


0.628 


1.541 


1.697 


1.709 


9186.03 


8661.17 


8853.21 


19469.83 


31370.98 


43892.02 


13955.6 


28107.7 


47191.5 


164289.8 


173116.1 


269901.7 


0.990 


0.988 


0.982 


1.513 


1.348 


1.741 


4.58E+06 


8.34E+06 


1.89E*07 


2.43E+06 


4.24E+06 


1.19E+07 


1.11 


1.20 


0.78 


1.89 


1.97 


1.60 


2.72 


2.69 


3.05 


0.81 


0.85 


0.80 
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Table 5. Characteristics of 2400 rpm synchronous machines 



nuiber of pole pairs 
power, hp 
efficiency factor 
voluae factor 
shaft rpi 
stator current density 
synchronous frequency 
rotor radius 
gap diaension 
back iron depth 
stator slot depth 
rotor slot depth 
stator slot factor 
rotor slot factor 
envelope voluae 
envelope weight 
hysteresis loss 
eddy current loss 
stator copper loss 
full load efficiency 
active length 
full load current density 
no load current density 
xs/turns-squared, p.u. 
internal volts/turn, p.u. 
overall length 
overall diaaeter 



1 


2 


3 


25775 


25775 


25775 


180000 


180000 


180000 


1286.1 


1286.1 


1286.1 


2400 


2400 


2400 


1.15E+07 


1.20E+07 


1.19E+07 


40 


80 


120 


0.1563 


0.1986 


0.2388 


0.0273 


0.0164 


0.0089 


0. 1094 


0.0695 


0.0557 


0.0929 


0.0695 


0.0557 


0.0416 


0.0215 


0.0088 


0.621 


0.750 


0.750 


0.676 


0.749 


0.750 


1.611 


1.112 


1.133 


11619.57 


7204.42 


6677.15 


11981.17 


14511.75 


21941.41 


2862.6 


6934.5 


15727.2 


444605.4 


427167.9 


482231.4 


0.977 


0.977 


0.974 


2.397 


1.709 


1.552 


2.11E+07 


3.11E+07 


5.39E+07 


7.75E+06 


1.28E+07 


2.11E+07 


2.00 


1.69 


1.83 


2.72 


2.43 


2.56 


3.13 


2.57 


2.54 


0.77 


0.71 


0.72 



4 


5 


6 


25775 


25775 


25775 


180000 


180000 


180000 


1286.1 


1286.1 


1286.1 


2400 


2400 


2400 


1.20E+07 


1.20E+07 


1.07E+07 


160 


200 


240 


0.2890 


0.3141 


0.3397 


0.0135 


0.0109 


0.0066 


0.0506 


0.0440 


0.0396 


0.0505 


0.0434 


0.0396 


0.0563 


0. 1047 


0.1040 


0.748 


0.750 


0.750 


0.432 


0.453 


0.327 


1.329 


1.436 


1.633 


6733.44 


7106.43 


7727.01 


27615.01 


32769.75 


47103.39 


26391.9 


39147.9 


67525.7 


232958.2 


143185.6 


114790.8 


0.985 


0.989 


0.988 


1.152 


1.144 


1.223 


1.73E+07 


7.83E+06 


7.99E+06 


9.64E+06 


4.55Et06 


4.32E+06 


1.01 


0.92 


1.07 


1.80 


1.72 


1.85 


2.36 


2.44 


2.61 


0.81 


0.82 


0.85 
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Table 6. Characteristics of 3000 rpm synchronous machines 



nusber of pole pairs 
poHer, hp 
efficiency factor 
voluie factor 
shaft rpi 
stator current density 
synchronous frequency 
rotor radius 
gap diiension 
back iron depth 
stator slot depth 
rotor slot depth 
stator slot factor 
rotor slot factor 
envelope voluse 
envelope weight 
hysteresis loss 
eddy current loss 
stator copper loss 
full load efficiency 
active length 
full load current density 
no load current density 
xs/turns-squared, p.u. 
internal volts/turn, p.u. 

overall length 
overall diaieter 



1 


0 

i. 


3 


25775 


25775 


25775 


180000 


180000 


180000 


1286.1 


1286.1 


1286.1 


3000 


3000 


3000 


1.20E+07 


1.20E+07 


1.20E+07 


50 


100 


150 


0.1494 


0. 1890 


0.2304 


0.0390 


0.0172 


0.0101 


0.1046 


0.0662 


0.0538 


0.0862 


0.0628 


0.0537 


0.0904 


0.0879 


0.0372 


0.750 


0.750 


0.750 


0.598 


0.494 


0.724 


1.205 


0.963 


0.972 


8815.49 


6498.77 


5801.99 


9574.79 


14072.41 


21089.10 


2859.6 


8405.7 


18895.4 


380285.7 


212641.9 


183124.5 


0.980 


0.988 


0.989 


1.670 


1.657 


1.363 


1.53E+07 


1.13E+07 


1.38E+07 


6.03E+06 


5.23E+06 


6.09E+06 


1.80 


1.41 


1.53 


2.53 


2.16 


2.27 


2.42 


2.48 


2.32 


0.76 


0.67 


0.70 



4 


5 


6 


25775 


25775 


25775 


180000 


180000 


180000 


1286.1 


1286.1 


1286. 1 


3000 


3000 


3000 


1.20E+07 


l.lOE+07 


1.18E+07 


200 


250 


300 


0.2432 


0.2872 


0.2974 


0.0052 


0.0068 


0.0088 


0.0426 


0.0402 


0.0347 


0.0422 


0.0402 


0.0347 


0.0151 


0.0784 


0.0427 


0.750 


0.750 


0.750 


0.730 


0.255 


0.750 


0.991 


1.203 


1.236 


5725.40 


6176.39 


6143.64 


31331.61 


41076.93 


46178.46 


37429.9 


61340.0 


82749.7 


204027.7 


140446.5 


134484.2 


0.986 


0.988 


0.986 


1.590 


1.306 


1.310 


2.42E+07 


1.40E+07 


1.06E+07 


9.63E+06 


7.44E+06 


6.96E+06 


1.78 


1.11 


0.69 


2.51 


1.89 


1.52 


2.58 


2.48 


2.53 


0.67 


0.75 


0.75 
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Table 7. Characteristics of 3600 rpm synchronous machines 



nuAber of pole pairs 
power , hp 
efficiency factor 
volute factor 
shaft rpt 
stator current density 
synchronous frequency 
rotor radius 
gap ditension 
back iron depth 
stator slot depth 
rotor slot depth 
stator slot factor 
rotor slot factor 
envelope volute 
envelope weight 
hysteresis loss 
eddy current loss 
stator copper loss 
full load efficiency 
active length 
full load current density 
no load current density 
xs/turns-squared, p.u. 
internal volts/turn, p.u. 

overall length 
overall diateter 



1 


2 


3 


25775 


25775 


25775 


180000 


180000 


180000 


1286.1 


1286.1 


1286.1 


3600 


3600 


3600 


1.16E+07 


1.20E+07 


1.20E+07 


60 


120 


180 


0.1413 


0.1824 


0.2191 


0.0139 


0.0120 


0.0123 


0.0989 


0.0638 


0.0511 


0.0492 


0.0638 


0.0486 


0.0496 


0.0541 


0.0572 


0.748 


0.750 


0.750 


0.733 


0.409 


0.528 


1.233 


0.814 


0.870 


9385.14 


5377.73 


5216.22 


14030.95 


15587.02 


22251.57 


5028.6 


11172.5 


23924.3 


190131.9 


228605. 1 


165930.4 


0.989 


0.987 


0.989 


3.258 


1.493 


1.372 


9.15E+06 


1.94E+07 


1.31E+07 


3.38E^06 


7.41E+06 


6.99E+06 


1.99 


1.89 


1.10 


2.71 


2.62 


1.88 


3.88 


2.27 


2.30 


0.61 


0.64 


0.66 



4 


5 


6 


25775 


25775 


25775 


180000 


180000 


180000 


1286.1 


1286. 1 


1286.1 


3600 


3600 


3600 


1.19E+07 


1.20E+07 


1.20E+07 


240 


300 


360 


0.2736 


0.2846 


0,2877 


0.0178 


0.0209 


0.0124 


0.0479 


0.0398 


0.0336 


0.0479 


0.0398 


0.0336 


0.0298 


0.0595 


0.0385 


0.750 


0.750 


0.750 


0.679 


0.691 


0.654 


1.065 


1.127 


1.102 


4979.26 


5177.33 


5190.73 


31042.24 


36284.98 


48590.52 


44501.0 


65021.1 


104486.5 


251863.9 


196473.3 


184008.6 


0.983 


0.985 


0.983 


0.890 


0.975 


1.166 


2.47E+07 


1.45E+07 


1.73E+07 


1.61E+07 


1.12E+07 


1.28E+07 


0.71 


0.43 


0.48 


1.54 


1.30 


1.35 


2.06 


2.20 


2.37 


0.77 


0.77 


0.73 



Table 8. Characteristics of 7200 rpm synchronous machines 



nuRber of pole pairs 
power, hp 
efficiency factor 
voluie factor 
shaft rpi 
stator current density 
synchronous frequency 
rotor radius 
gap diiension 
back iron depth 
stator slot depth 
rotor slot depth 
stator slot factor 
rotor slot factor 
envelope voIuie 
envelope weight 
hysteresis loss 
eddy current loss 
stator copper loss 
full load efficiency 
active length 
full load current density 
no load current density 
xs/turns-squared, p.u. 
internal volts/turn, p.u. 

overall length 
overall diaieter 



1 


2 


3 


25775 


25775 


25775 


180000 


180000 


180000 


1280.1 


1286.1 


1286.1 


7200 


7200 


7200 


1.20E+07 


1.20E+07 


1.16E+07 


120 


240 


360 


0.1159 


0.1710 


0.1958 


0.0096 


0.0093 


0.0053 


0.0812 


0.0599 


0.0457 


0.0405 


0.0452 


0.0457 


0.0329 


0.0189 


0.0105 


0.750 


0.750 


0.749 


0.511 


0.617 


0.661 


0.711 


0.544 


0.525 


5341.21 


3476.08 


2971.92 


17305.86 


22097.70 


28944.03 


12404.5 


31678.5 


62239.7 


149625.0 


159505.1 


162532.1 


0.991 


0.989 


0.987 


2.863 


1.212 


0.972 


1.65E+07 


2.71E+07 


3.92E+07 


6. 14E+06 


1.16E+07 


1.46E+07 


1.97 


1.59 


1.97 


2.69 


2.33 


2.69 


3.36 


1.93 


1.78 


0.49 


0.57 


0.58 



4 


5 


6 


25775 


25775 


25775 


180000 


180000 


180000 


1286.1 


1286.1 


1286.1 


7200 


7200 


7200 


1.20E+07 


1.20E+07 


1.20E+07 


480 


600 


720 


0.2044 


0.2219 


0.1932 


0.0093 


0.0046 


0.0055 


0.0358 


0.0311 


0.0225 


0.0358 


0.0311 


0.0225 


0.0714 


0.0138 


0.0292 


0.750 


0.750 


0.750 


0.633 


0.468 


0.570 


0.544 


0.572 


0.564 


3136.32 


2950.31 


3427.07 


30808.73 


50465.33 


64447.34 


88332.6 


180863.2 


277168.2 


74543.2 


143253.5 


117812.4 


0.990 


0.981 


0.977 


1.081 


1.082 


1.952 


7.81E+06 


3.70E+07 


1.72E+07 


5.00E+06 


1.99E+07 


1.28E+07 


0.74 


1.08 


0.48 


1.56 


1.86 


1.34 


1.94 


1.99 


2.75 


0.57 


0.58 


0.49 
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Table 9. Characteristics of 180 rpm synchronous machines 



nuiber of pole pairs 


1 


2 


3 


4 


5 


power, hp 


25775 


25775 


25775 


25775 


25775 


efficiency factor 


180000 


180000 


180000 


180000 


180000 


Yoluie factor 


1286.1 


1286.1 


1286.1 


1286.1 


1286.1 


shaft rpa 


180 


180 


180 


180 


180 


stator current density 


1.18E+07 


1.12E+07 


1.20E+07 


9.25E+06 


1.20E+07 


synchronous frequency 


3 


6 


9 


12 


15 


rotor radius 


0.2641 


0.3738 


0.4588 


0.5869 


0.5507 


gap diaension 


0.0331 


0.0978 


0.0723 


0.0317 


0.0198 


back iron depth 


0.1849 


0. 1308 


0.1071 


0.1027 


0.0771 


stator slot depth 


0.0592 


0.1308 


0.1071 


0. 1027 


0.0580 


rotor slot depth 


0.0500 


0.1867 


0.1934 


0.0565 


0.1079 


stator slot factor 


0.741 


0.750 


0.750 


0.742 


0.750 


rotor slot factor 


0.558 


0.571 


0.422 


0.535 


0.287 


envelope voluae 


16.305 


9.329 


9.048 


11.378 


10.306 


envelope weight 


119110.5 


56300.0 


48896.8 


56398.0 


56899.5 


hysteresis loss 


10465.3 


6915.1 


9833.5 


18267.4 


23999.1 


eddy current loss 


187.5 


247.8 


528.6 


1309.4 


2150.3 


stator copper loss 


2.25E+06 


2.06E^ 


1.76E+06 


1.34E+06 


1.18E+06 


full load efficiency 


0.895 


0.903 


0.916 


0.934 


0.941 


active length 


14.913 


3.135 


2.590 


2.376 


3.708 


full load current density 


1.53E+07 


1.25E+07 


1.33E+07 


1.91E+07 


1.42E+07 


no load current density 


5.61E^06 


6.14E^06 


7.25E+06 


8.93E+06 


7.27E+06 


xs/turns'squared, p.u. 


2.00 


1.28 


1.04 


1.38 


1.19 


internal volts/turn, p.u. 


2.72 


2.04 


1.83 


2.14 


1.96 


overall length 


16.10 


5.02 


4.71 


4.85 


5.99 


overall diaaeter 


1.08 


1.47 


1.49 


1.65 


1.41 
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Table 10. Characteristics of more 180 rpm synchronous 
machines 



nutber of pole pairs 


6 


7 


8 


9 


10 


poMer, hp 


25775 


25775 


25775 


25775 


25775 


efficiency factor 


180000 


180000 


180000 


180000 


180000 


volute factor 


1286.1 


1286. 1 


1286.1 


1286.1 


1286.1 


shaft rpt 


180 


180 


180 


180 


180 


stator current density 


1.20E+07 


1.20E+07 


1.20E+07 


1.20E+07 


1.20E+07 


synchronous frequency 


18 


21 


24 


27 


30 


rotor radius 


0.6176 


0.5729 


0.7152 


0.6943 


0.7296 


gap ditension 


0.0647 


0.0335 


0.0229 


0.013B 


0.0396 


back iron depth 


0.0720 


0.0573 


0.0626 


0.0540 


0.0511 


stator slot depth 


0.0720 


0.0569 


0.0626 


0.0540 


0.0509 


rotor slot depth 


0.1652 


0.0855 


0.2046 


0.0740 


0.1420 


stator slot factor 


0.750 


0.750 


0.750 


0.750 


0.750 


rotor slot factor 


0.368 


0.750 


0.374 


* 0.283 


0.615 


envelope volute 


11.693 


10.486 


12.870 


12.385 


14.273 


envelope weight 


48739.8 


53234.1 


52286.9 


53189.2 


56440.5 


hysteresis loss 


22385.0 


27947.6 


31793.9 


42718.3 


42086.1 


eddy current loss 


2406.8 


3505.7 


4557.9 


6889.4 


7541.6 


stator copper loss 


1.86E+06 


1.24E+06 


7.69E+05 


1.14E+06 


1.11E^06 


full load efficiency 


0.911 


0.938 


0.960 


0.942 


0.943 


active length 


2.225 


3.418 


2.045 


2.549 


2.365 


full load current density 


1.74Et07 


l.llE+07 


7.09E+06 


1.98E+07 


9.69E+06 


no load current density 


1.29E+07 


7.98E+06 


4.18E+06 


1.07E+07 


7.76E+06 


xs/turns-squared, p.u. 


0.48 


0.53 


0.89 


1.08 


0.36 


internal volts/turn, p.u. 


1.35 


1.39 


1.69 


1.86 


1.25 


overall length 


4.95 


5.84 


5.00 


5.38 


5.44 


overall diateter 


1.65 


1.44 


1.73 


1.63 


1.74 
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3.4. Verification 



Data on a large turbogenerator, Big Sandy Unit Two, was 
available in reference sixteen; it was used to verify the 
synchronous machine design program. Big Sandy is rated at 
907,000 kVA, power factor 0.9, at a rated voltage of 26 kV. 
When the machine parameters were input to the synchronous 
design program, it produced a machine very close to Big 
Sandy. It was judged that the design program would yield 
good results. 
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Chapter Four. Permanent Magnet Machines 



Permanent magnet machines are very similar to 
synchronous machines. The main difference lies in the 
method used to produce the field flux wave. Instead of a 
field current causing the wave, permanent magnets provide 
the flux. No exciter is used with these machines and there 
are no brushes . 

4.1. Magnet material 

Many different elements may be used to manufacture per- 
manent magnets. Past designs used ceramics, aluminum- 
nickel-cobalt-iron-titanium (AlNiCo), and samarium-cobalt 
(SmCo). However, ceraunic magnets do not produce sufficient 
residual flux (see Appendix C for an explanation of terms), 
and any magnet based on cobalt is high in cost and may be in 
limited supply. Recently, magnets of neodymium-iron-boron 
(NdFeB), have entered the marketplace. None of the con- 
stituents of the NdFeB magnets are strategic materials; it 
is expected that availability and cost will improve. 

NdFeB magnets have a high Maximum Energy Product (MEP) 
that may be used advantageously by the machine designer. 

Data from Sumitomo Special Metals [17] indicates their 
NEOMAX line have MEPs as high as 37 MGOe, which is higher 
than the 30 MGOe of the SmCo magnets marketed as REC-30 by 
TDK Corporation [18]. High MEP is not the only criteria for 
magnet selection; flux stability, cost-to-perf ormance ratio, 
ease of machine assembly, and other characteristics may en- 
ter the decision process. This study needed the best per- 
formance of its machines, so NdFeB was selected on the basis 
of its high MEP. Thermal stability was assumed to be satis- 
factory if the thermal considerations in Chapter Two were 
met. 

Cost has been mentioned above. Magnet material is sig- 
nificantly more expensive than magnet steel and copper. 
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Magnets in the quantity used by a large production run of 
25,000 hp machines might cost as much as $120 per pound, 
compared with the 58<(; per pound of M19 steel. 1 5 Obviously, 
permanent magnet machines will be more expensive, but the 
cost of magnet materials may be made part of the optimiza- 
tion process. The degree of magnet overhang, discussed in 
Appendix C, also affects cost. 

4.1.1. Magnet cost factor 

A change to the objective function was made to incor- 
porate the cost of magnet material relative to magnet steel. 
The ratio of the above costs was taken and the result called 
the magnet cost factor, km. The objective function was 
modified to 

Elective weight = (weight * kiliagnet weight)) * Iced - edcy) * kvivoluie * kwlaagnet volute)) 

An initial value for km of 170 was used, and several 
machines designed. Then a value of km = 25 was tried. The 
machines with km = 170 indeed had less magnet material in 
them, but at a cost. The change to km = 25 resulted in a 
larger machine (23.5%) with a lower stator current density, 
20% more magnet material, and about a 1.5% increase in 
machine efficiency. That 1.5% translates to a lot of fuel 
aboard a ship, so it was decided to use km = 25. The extra 
magnet material will add about $22,000 per machine. 

4.2. Assumptions 

The largest obstacle to assembling high-power permanent 
magnet machines is their inherent residual magnetism. If 

14. Estimate by Mr. Yokokura, President of Sumitomo Special 
Metals of America. 

15. Book price for 26 gauge M19 steel from Mr. Dagg of USX. 
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the magnets possess all of their properties at assembly, it 
will be extremely difficult to place the rotor (which con- 
tains the magnets) inside the stator. The rotor would be 
strongly attracted to the iron of the stator. Of course, 
the magnets may be magnetized after machine assembly, but it 
may be difficult to achieve MEP without elevated tempera- 
tures inside the machine. The assumption is made here that 
the magnets will be magnetized prior to assembly. The 
detailed design of the machine will have to include con- 
sideration of the jigs and fixtures necessary for assembly. 

The only other assumption worthy of mention is that the 
load line may be modeled as described in Appendix C. 

4.3. Machine description 

The same rpm and pole-pair combinations were used as in 
the synchronous machines. The magnets on the rotor are ar- 
ranged in a cylindrical-wedge configuration, as shown in the 
figure below. 

Figure 7. End view of a permanent magnet machine 




MrtOiNET MflTERI/M. 



MA6NET STe.e\_ 

BiNces. 



The rotor slot factor (Ir) used in the other types of 
machines is called here the magnet slot factor and refers to 
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the width of magnet per pole pitch. It may vary between 25% 
and 75% of the pole pitch, the same as for rotor slots in 
other machines . 

The rotor slot depth, ds, does not exist in this 
machine. In this case, the rotor radius, r, is added to Im, 
the magnet radial dimension, to find the actual width of the 
rotating core. (In the synchronous and induction machines, 
dr was included in the rotor radius.) 

4.3.1. Efficiency 

For the higher rpm machines, efficiency within a par- 
ticular rpm group decreased with an increasing number of 
poles. The most efficient machines, at about 99%, had four 
poles. This is higher than the synchronous or induction 
machines, largely because there are no rotor copper losses. 
The twelve-pole efficiencies were about 98%, which is not 
too large a spread. 

The 180 rpm machines had a fairly flat efficiency curve 
(excepting one anomaly) up to about a 28-pole machine, where 
efficiency started to vary widely. There, the conflict be- 
tween the number of pole-pairs and maximum rotor radius 
started to become significant. The flat efficiency was 
about 96.5%, which is less efficient than the higher rpm 
machines . 

4.3.2. Weight and volume 

The higher rpm machines had a general tendency toward 
lower weight and volume as rpm increased. Within an rpm 
group, the four- and six-pole machines had the lowest weight 
and volume. The smallest machines were larger than the 
synchronous machines. 
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EFFICIENCY 



Figure 8. High rpm permanent magnet machine efficiency 



PERMANENT MAGNET MOTOR EFFICIENCY 
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EFFICIENCY 



Figure 9. 180 rpm permanent magnet efficiency 



PERMANENT MAGNET MOTOR EFFICIENCY 



DIRECT-DRIVE, 180 RPM, 25775 HP 
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VOLUME (cubic meters) 



Figure 10. High rpm permanent magnet machine volume 



PERMANENT MAGNET MACHINE VOLUME 



25775 HP 




NUMBER CF POLE-PAIRS 
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VOLUME (M'^^) 



Figure 11. 180 rpm permanent magnet machine volume 



PERMANENT MAGNET MOTOR VOLUME 



DIRECT-DRIVE. 180 RPM, 25775 HP 
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WEIGHT (KG) 

(Thout>ancls) 



Figure 12. High rpm permanent magnet machine weight 



PERMANENT MAGNET MOTOR WEIGHT 



25775 HP 




NUMBER OF POLE- PAIRS 
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WEIGHT (KG) 
(Tbousonds) 



Figure 13. 180 rpm permanent magnet machine weight 



PERMANENT MAGNET MOTOR WEIGHT 



DIRECT-DRIVE, 180 RPM, 25775 HP 
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The weight and volume of the 180 rpm machines increased 
with the number of pole-pairs. Again, the tip speed limita- 
tion is linked to this increase, as with the synchronous 
machines . There was a wide variation around the general in- 
crease, and these machines are large relative to all others. 
They are not competitive as ship propulsion motors because 
of their size. 

4.4. Machine characteristics 

The following tables give the machine characteristics 
of the permanent magnet machines designed for this study. 

The stator slot factor. Is, tended to the maximum of 0.75 
for the higher rpm and 180 rpm machines. This occurred as 
the optimization algorithm tried to minimize envelope volume 
and weight. The rotor slot factor, Ir, was a constant 0.378 
for all machines. The rotor slot factor was calculated to 
produce load-line (MEP) operating flux, as derived in Appen- 
dix B. With a different magnet material selection (and a 
consequent change in operating point flux), a different Ir 
would have resulted. 

The magnet overhang tended toward the maximum limit, in 
an attempt to achieve the highest flux levels. Permanent 
magnet machines cannot rival the flux level produced by the 
field winding of a synchronous machine, but the optimization 
algorithm did its best. 

The per-unit synchronous reactance-per-turns-squared 
was limited to 3.0 in these machines. This reactance tended 
toward the limit, but was lower with an increasing number of 
pole-pairs in the 180 rpm machines. It was very difficult 
to achieve valid designs with xsroax = 2.0, as in the 
synchronous machines. If an xsmax greater than 2.0 is unac- 
ceptable, these machines will be less competitive. 

The amount of magnet material varied from 50 kg to a 
few hundred kg in the higher rpm machines to 700-4000 kg in 
the 180 rpm machines. The cost of 700 kg of NdFeB is about 
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Table 11. Character i st i cs o-f 1800 rpm magnet machine® 



nuiber of pole pairs 
poNer, hp 
efficiency factor 
voluie factor 
•agnet factor 
shaft rpi 
stator current density 
synchronous frequency 
rotor radius 
gap dimension 
back iron depth 
stator slot depth 
lagnet radial diiension 
stator slot factor 
rotor slot factor 
envelope voluie 
envelope Meight 
hysteresis loss 
eddy current loss 
stator copper loss 
full load efficiency 
active length 
•agnet Neight 
■agnet voluie 
xs/turns-squared, p.u. 
■agnet overhang 
air gap flux density 
overall length 
overall diaieter 



1 


2 


3 


25775 


25775 


25775 


180000 


180000 


180000 


128&.1 


1286.1 


1286. 1 


25 


25 


25 


1800 


1800 


1800 


7 . 49 E +06 


1 . 20 E +07 


7 . 32 E +06 


30 


60 


90 


0.2887 


0.3396 


0.3643 


0.0326 


0.0382 


0.0067 


0.1332 


0.0870 


0.0530 


0.0687 


0.0742 


0.0456 


0.0342 


0.0401 


0.0070 


0.750 


0.750 


0.750 


0.378 


0.378 


0.378 


3.548 


2.558 


3.309 


24280.7 


12528.8 


19366.5 


20581.5 


20298.5 


52315.8 


3688.1 


7274.8 


28124.3 


171150.0 


265272.5 


93554.8 


0.990 


0.985 


0.991 


2.328 


0.876 


2.833 


405.184 


262.376 


157.355 


0.055 


0.035 


0.021 


2.994 


2.339 


2.996 


0.001 


0.298 


0.315 


0.619 


0.687 


0.642 


3.750 


2.548 


4.345 


1.047 


1.078 


0.939 



4 


5 


6 


25775 


25775 


25775 


180000 


180000 


180000 


1286.1 


1286.1 


1286.1 


25 


25 


25 


1800 


1800 


1800 


9 . 14 E +06 


8 . 00 E +06 


1 . 04 E +07 


120 


150 


180 


0.5964 


0.5544 


0.6068 


0.0188 


0.0067 


0.0123 


0.0742 


0.0502 


0.0482 


0.0353 


0.0263 


0.0277 


0.0198 


0.0071 


0.0130 


0.750 


0.750 


0.646 


0.378 


0.378 


0.378 


6.304 


5.847 


6.319 


20910.3 


24826.4 


20948.1 


80505.5 


121495.4 


124368.0 


57704.9 


108857.1 


133717.2 


98867.0 


69323.3 


85961.4 


0.988 


0.985 


0.982 


0.933 


1.888 


1.257 


273.763 


178.613 


246.460 


0.037 


0.024 


0.033 


1.165 


1.611 


0.920 


0.326 


0.331 


0.329 


0.722 


0.671 


0.700 


3.473 


4.160 


3.786 


1.450 


1.275 


1.390 
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Table 12. Characteristics of 2400 rpm magnet machines 



nufiber of pole pairs 
poHer, hp 
efficiency factor 
voluie factor 
•agnet factor 
shaft rpi 
stator current density 
synchronous frequency 
rotor radius 
gap diiension 
back iron depth 
stator slot depth 
■agnet radial diiension 
stator slot factor 
rotor slot factor 
envelope volute 
envelope weight 
hysteresis loss 
eddy current loss 
stator copper loss 
full load efficiency 
active length 
lagnet Height 
■agnet voluee 
xs/turns-squared, p.u. 
■agnet overhang 
air gap flux density 
overall length 
overall diaieter 



1 


2 


3 


25775 


25775 


25775 


180000 


180000 


180000 


1286.1 


1286.1 


1286.1 


25 


25 


25 


2400 


2400 


2400 


12000000 


12000000 


12000000 


40 


80 


120 


0.3594 


0.0005 


0.4018 


0.0371 


0.1883 


0.0136 


0.1640 


0.0420 


0.0625 


0.0697 


0.0249 


0.0341 


0.0390 


0.1978 


0.0143 


0.348 


0.750 


0.750 


0.378 


0.378 


0.378 


4.497 


1.342 


2.727 


26548.1 


8797.0 


12335.3 


34028.8 


18160.1 


46136.0 


8130.4 


8677.9 


33069.4 


234324.1 


217110.2 


127316.4 


0.986 


0.987 


0.989 


1.536 


4.392 


1.291 


384.207 


10.204 


167.444 


0.052 


0.001 


0.023 


2.417 


0.282 


1.973 


O.OIB 


0.330 


0.331 


0.617 


0.635 


0.676 


3.278 


5.938 


3.010 


1.260 


0.511 


1.024 



4 


5 


6 


25775 


25775 


25775 


180000 


180000 


180000 


1286.1 


1286.1 


1286. 1 


25 


25 


25 


2400 


2400 


2400 


10202760 


11959030 


11997423 


160 


200 


240 


0.4408 


0.3361 


0.4760 


0.0088 


0.0032 


0.0165 


0.0504 


0.0290 


0.0379 


0.0325 


0.0234 


0.0300 


0.0092 


0.0033 


0.0173 


0.750 


0.749 


0.678 


0.378 


0.378 


0.378 


3.157 


2.349 


3.438 


13510.2 


13660.4 


11877.1 


68103.9 


86724.5 


91478.9 


65087.5 


103604.3 


131141.0 


90617.2 


112885.8 


91721.4 


0.988 


0.984 


0.984 


1.387 


3.061 


1.129 


128.106 


73.723 


217.676 


0.017 


0.010 


0.029 


1.965 


2.879 


0.751 


0.315 


0.338 


0.324 


0.672 


0.640 


0.692 


3.222 


4.431 


3.168 


1.065 


0.783 


1.121 
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Table 13. Characteristics of 3000 rpm magnet machines 



nuiber of pole pairs 
poHer, hp 
efficiency factor 
volute factor 
tagnet factor 
shaft rpt 
stator current density 
synchronous frequency 
rotor radius 
gap diiension 
back iron depth 
stator slot depth 
tagnet radial ditension 
stator slot factor 
rotor slot factor 
envelope volute 
envelope Height 
hysteresis loss 
eddy current loss 
stator copper loss 
full load efficiency 
active length 
tagnet Height 
tagnet volute 
xs/turns-squared, p.u. 
tagnet overhang 
air gap flux density 
overall length 
overall diateter 



1 


2 


3 


4 


5 


6 


25775 


25775 


25775 


25775 


25775 


25775 


180000 


180000 


180000 


180000 


180000 


180000 


1286.1 


1286.1 


1286.1 


1286.1 


1286.1 


1286.1 


25 


25 


25 


25 


25 


25 


3000 


3000 


3000 


3000 


3000 


3000 


10639824 


9557320 . 


12000000 


11999567 


12000000 


11400662 


50 


100 


150 


200 


250 


300 


0.2480 


0.2394 


0.2239 


0.0003 


0.3116 


0.4837 


0.0227 


0.0075 


0.0035 


0.1603 


0.0027 


0.0375 


0.1121 


0.0510 


0.0317 


0.0175 


0.0267 


0.0428 


0.0405 


0.0369 


0.0235 


0.0074 


0.0224 


0.0312 


0.0238 


0.0079 


0.0037 


0.1683 


0.0028 


0.0394 


0.750 


0.750 


0.750 


0.750 


0.747 


0.681 


0.378 


0.378 


0.378 


0.378 


0.378 


0.378 


2.213 


1.834 


1.771 


2.143 


1.950 


3.643 


15051.7 


12394.7 


12381.0 


8254.8 


11507.2 


9968.3 


22692.8 


36848.2 


57339.4 


45284.3 


91118.4 


95620.6 


6777.4 


22010.1 


51374.8 


54098.3 136066.7 


171347.8 


157512.6 


108787.2 


132086.7 


170766.7 


97680.0 


74849.0 


0.990 


0.991 


0.988 


0.986 


0.983 


0.983 


2.395 


3.719 


5.490 


16.723 


3.005 


0.732 


250.152 


123.693 


91.981 


15.685 


55.765 


338.030 


0.034 


0.017 


0.012 


0.002 


0.008 


0.046 


2.934 


2.993 


2.997 


0.043 


2.991 


0.351 


0.007 


0.009 


0.326 


0.325 


0.311 


0.336 


0.618 


0.619 


0.627 


0.623 


0.637 


0.736 


3.574 


4.738 


6.415 


18.039 


4.274 


2.975 


0.847 


0.669 


0.565 


0.371 


0.727 


1.191 



Table 14. Characteristics of 3600 rpm magnet machines 



nuiber of pole pairs 
poMer, hp 
efficiency factor 
voluie factor 
•agnet factor 
shaft rpi 
stator current density 
synchronous frequency 
rotor radius 
gap diiension 
back iron depth 
stator slot depth 
•agnet radial diiension 
stator slot factor 
rotor slot factor 
envelope voluie 
envelope weight 
hysteresis loss 
eddy current loss 
stator copper loss 
full load efficiency 
active length 
•agnet weight 
•agnet voluae 
xs/turns-squared, p.u. 
•agnet overhang 
air gap flux density 
overall length 
overall diaieter 



1 


2 


3 


25775 


25775 


25775 


180000 


180000 


180000 


1286.1 


1286.1 


1286.1 


25 


25 


25 


3600 


3600 


3600 


12000000 


12000000 


11966619 


60 


120 


180 


0.3044 


0.2596 


0.3006 


0.0228 


0.0088 


0.0121 


0.1454 


0.0574 


0. 0462 


0.0332 


0.0442 


0.0411 


0.0240 


0.0092 


0.0127 


0.750 


0.553 


0.750 


0.378 


0.378 


0.378 


2.475 


1.602 


1.427 


15468.5 


9853.6 


7072.7 


29048.1 


36290.0 


37383. 1 


10410.6 


26012.0 


40193.3 


163951.5 


114166.3 


104037.5 


0.990 


0.991 


0.991 


1.394 


2.276 


1.279 


221.175 


113.409 


104.478 


0.030 


0.015 


0.014 


2.960 


2.977 


2.067 


0.336 


0.324 


0.327 


0.664 


0.641 


0.663 


2.799 


3.386 


2.581 


1.012 


0.740 


0.800 



4 


5 


6 


25775 


25775 


25775 


180000 


180000 


180000 


1286.1 


1286.1 


1286.1 


25 


25 


25 


3600 


3600 


3600 


11677983 11998747 11975558 


240 


300 


360 


0.2526 


0.3700 


0.4536 


0.0025 


0.0063 


0.0173 


0.0268 


0.0332 


0.0371 


0.0208 


0.0223 


0.0240 


0.0026 


0.0067 


0.0182 


0.735 


0.750 


0.749 


0.378 


0.378 


0.378 


1.694 


2.063 


2.801 


11238.1 


9644.1 


8977.4 


84962.0 


92749.9 104238.5 


121798.6 


166203.6 224148.8 


95474.6 


73859.8 


66008. 1 


0.985 


0.983 


0.980 


4.322 


1.669 


0.907 


58.252 


91.121 


175.220 


0.008 


0.012 


0.024 


2.981 


1.490 


0.593 


0.333 


0.340 


0.339 


0.630 


0.661 


0.708 


5.352 


3.202 


2.863 


0.605 


0.864 


1.064 
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Table 15. Characteristics of 7200 rpm magnet machines 



nuiber of pole pairs 
power, hp 
efficiency factor 
voluse factor 
•agnet factor 
shaft rpi 
stator current density 
synchronous frequency 
rotor radius 
gap diiension 
back iron depth 
stator slot depth 
iagnet radial diiension 
stator slot factor 
rotor slot factor 
envelope voluie 
envelope weight 
hysteresis loss 
eddy current loss 
stator copper loss 
full load efficiency 
active length 
iagnet weight 
iagnet voluie 
xs/turns-squared, p.u. 
iagnet overhang 
air gap flux density 
overall length 
overall diaieter 



1 


2 


3 


25775 


25775 


25775 


180000 


180000 


180000 


1286.1 


1286.1 


1286.1 


25 


25 


25 


7200 


7200 


7200 


12000000 12000000 8393050 . 


120 


240 


360 


0.2266 


0.1938 


0.2622 


0.0146 


0.0059 


0.0095 


0.1045 


0.0423 


0.0392 


0.0285 


0.0284 


0.0315 


0.0153 


0.0062 


0.0099 


0.749 


0.750 


0.750 


0.378 


0.378 


0.378 


1.247 


0.811 


1.113 


8186.9 


5319.9 


6052.8 


30513.8 


38176.1 


65438.7 


21871.7 


54727.9 140716.0 


90777.3 


71500.7 


38515.8 


0.993 


0.992 


0.987 


1.552 


2.384 


1.621 


110.276 


57.403 


87.115 


0.015 


0.008 


0.012 


2.946 


2.884 


1.258 


0.322 


0.331 


0.307 


0.648 


0.635 


0.647 


2.578 


3.207 


2.747 


0.748 


0.541 


0.685 



4 


5 


6 


25775 


25775 


25775 


180000 


180000 


180000 


1286.1 


1286.1 


1286.1 


25 


25 


25 


7200 


7200 


7200 


12000000 12000000 12000000 


480 


600 


720 


0.2627 


0.2467 


0.2418 


0.0041 


0.0020 


0.0022 


0.0289 


0.0209 


0.0171 


0.0244 


0.0114 


0.0104 


0.0043 


0.0021 


0.0023 


0.750 


0.748 


0.750 


0.378 


0.378 


0.378 


0.931 


1.356 


1.398 


4898.9 


8797.5 


9221.7 


72593.1 


172996.0 218032.2 


208133.7 620002.2 937689.5 


52482.8 


49628.2 


49815.1 


0.983 


0.958 


0.941 


1.542 


3.965 


4.508 


36.486 


42.995 


50.295 


0.005 


0.006 


0.007 


2.246 


1.608 


1.142 


0.328 


0.338 


0.319 


0.650 


0.631 


0.629 


2.627 


4.968 


5.493 


0.640 


0.562 


0.543 
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Table 16. Characteristics of 180 rpm magnet machines 



nuiber of pole pairs 


1 


2 


3 


4 


5 


poMer, hp 


25775 


25775 


25775 


25775 


25775 


efficiency factor 


180000 


180000 


180000 


180000 


180000 


voluie factor 


1288.1 


1288.1 


1288.1 


1288.1 


1288.1 


iagnet factor 


25 


25 


25 


25 


25 


shaft rp§ 


180 


180 


180 


180 


180 


stator current density 


7 . 80 E +06 1 . 18 E +07 


1 . 20 E +07 


1 . 09 E +07 


1 . 20 E +07 


synchronous frequency 


3 


8 


9 


12 


15 


rotor radius 


0.4399 


0.1787 


0.9745 


0.7389 


0.8993 


gap diiension 


0.0413 


0.0038 


0.0201 


0.0188 


0.0214 


back iron depth 


0.1994 


0.0377 


0.1491 


0.0779 


0.0885 


stator slot depth 


0.0932 


0.0282 


0.0327 


0.0478 


0.0830 


iagnet radial diiension 


0.0433 


0.0038 


0.0212 


0.0174 


0.0225 


stator slot factor 


0.483 


0.531 


0.750 


0.750 


0.740 


rotor slot factor 


0.378 


0.378 


0.378 


0.378 


0.378 


envelope volute 


29.375 


38.894 


34.709 


20.521 


22.314 


envelope weight 


211948.9 284972.1 


158312.7 


107071.7 


77934.3 


hysteresis loss 


20311.8 


53970.7 


47089.7 


41032.7 


38928.0 


eddy current loss 


384.0 


1934.3 


2531.5 


2941.1 


3308.5 


stator copper loss 


7 . 90 E +05 2 . 54 E +08 8 . 97 E +05 


7 . 01 E +05 


7 . 11 E +05 


full load efficiency 


0.980 


0.881 


0.983 


0.983 


0.982 


active length 


12.100 


171.774 


3.195 


4.580 


1.888 


iagnet weight 


4072.4 


2079.8 


1941.0 


1048.8 


1087.7 


iagnet volute 


0.550 


0.281 


0.282 


0.142 


0.147 


xs/turns-squared, p.u. 


2.998 


2.990 


2.981 


2.983 


2.870 


iagnet overhang 


0.003 


0.034 


0.332 


0.010 


0.339 


air gap flux density 


0.819 


0.819 


0.874 


0.818 


0.704 


overall length 


14.198 


172.518 


7.258 


7.851 


5.839 


overall diaieter 


1.548 


0.498 


2.353 


1.782 


2.140 
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Table 17. Characteristics of more 180 rpm magnet machines 



nuiber of pole pairs 


6 


7 


8 


9 


10 


poNer, hp 


25775 


25775 


25775 


25775 


25775 


efficiency factor 


180000 


180000 


180000 


180000 


180000 


voluie factor 


1286.1 


1286.1 


1286.1 


1286.1 


1286.1 


■agnet factor 


25 


25 


25 


25 


25 


shaft rpi 


180 


180 


180 


180 


180 


stator current density 


6.03E+06 1.20E+07 


1.20E+07 9.42E+06 


1.20E+07 


synchronous frequency 


18 


21 


24 


27 


30 


rotor radius 


0.6047 


1.2439 


0.9696 


0.9792 


0.9526 


gap diaension 


0.0139 


0.1161 


0.0089 


0.0031 


0.0079 


back iron depth 


0.0426 


0.1117 


0.0546 


0.0469 


0.0430 


stator slot depth 


0.0337 


0.0627 


0.0312 


0.0255 


0.0360 


■agnet radial diaension 


0.0146 


0.1219 


0.0093 


0.0033 


0.0083 


stator slot factor 


0.536 


0.750 


0.748 


0.575 


0.748 


rotor slot factor 


0.378 


0.378 


0.378 


0.378 


0.378 


envelope voluae 


44.578 


53.524 


29.253 


43.668 


26.292 


envelope weight 


313139.6 


98601.8 


120336.5 234296.3 


102286.2 


hysteresis loss 


186214.5 


67518.3 


96060.7 214344.0 


101073.4 


eddy current loss 


20021.3 


8469.3 


13770.9 


34568.5 


18111.9 


stator copper loss 


3.91E+05 7.15E+05 


5.34E+05 3.91E+05 


5.31E+05 


full load efficiency 


0.970 


0.960 


0.968 


0.968 


0.967 


active length 


24.182 


0.651 


3.523 


7.417 


3.167 


aagnet weight 


3800.1 


3431.9 


932.5 


704.3 


714.1 


aagnet voluae 


0.514 


0.464 


0.126 


0.095 


0.097 


xs/turns-squared, p.u. 


0.546 


0.525 


2.350 


2.990 


2.391 


aagnet overhang 


0.010 


0.329 


0.336 


0.340 


0.325 


air gap flux density 


0.619 


0.859 


0.669 


0.644 


0.671 


overall length 


26.715 


6.579 


7.474 


11.360 


7.042 


overall diaaeter 


1.390 


3.069 


2.128 


2.109 


2.079 
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$150,000, rendering the 180 rpm machines less economic to 
build. A machine such as these may cost between four and 
eight million dollars. Only those machines with magnet 
costs that are reasonable with respect to the other material 
costs should be candidates for design. 

In all the permanent magnet machines, stator current 
density went to the maximum. The overall length and overall 
diameter of the more reasonable machines would allow them to 
fit in machinery spaces aboard a ship. 

4.5. Verification 

No high-power permanent magnet machines were discovered 
during the search to find a benchmark. Because of the high 
material cost and the competition afforded by synchronous 
and induction machines, it seems none have been built. 
Several paper studies were found [13, 19, 20, 21, and 22], 
and the parameters resulting from this computer modeling 
seem to agree with them. The machine size is what was ex- 
pected, given the lower air gap flux density. The ef- 
ficiency was higher than the synchronous and induction 
machines. All-in-all, this modeling gave good machines. 
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Chapter Five. Induction machines 



The stator of an induction machine is the same as those 
of synchronous and permanent magnet machines. The rotor is 
significantly different. There is no independent mechanism 
to produce a rotor flux wave. The rotor winding is shorted, 
whether it is wound or cast, so that as the stator flux wave 
passes over the rotor, currents are induced in the winding. 
These currents produce only a small reaction flux, but it 
still tends to align with the stator flux wave. When at 
operating speeds, the rotor speed is a bit slower than the 
stator flux wave speed, and the difference in speeds is 
called slip. Typically, slip is a few percent of the stator 
frequency. The rotor currents are at slip frequency. If 
the rotor and stator speeds were the saune, slip would be 
zero, there would be no tendency to align and torque would 
be zero. Then, the rotor would lag behind the stator until 
current was induced in the rotor winding by the passing 
stator flux wave and torque was again produced. 

If solid bars are used as the rotor winding, they are 
shorted at the ends of the rotor by end rings, to form what 
is called a "squirrel cage" rotor. If actual turns are 
used, the winding may be shorted through external resis- 
tances to affect the starting and torque-slip characteris- 
tics of the machine. Fitzgerald et al [23] and Alger [24] 
discuss induction machine characteristics in some detail. 

5.1. Assumptions 

A squirrel cage rotor was assumed for these machines. 
Copper was designated as the material for the rotor bars. 
However, these machines will be fed from a frequency 
changer, so only one layer of bars was used and the effects 
of magnetic diffusion ignored in the analysis (see 
Appendix D for a derivation of the components of an induc- 
tion machine equivalent circuit). The number of rotor bars 
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was arbitrarily set at 71. This quantity should not cause 
undesirable harmonics, as it will not be an integral mul- 
tiple of the number of poles or stator slots in any machine. 
The number and width of the rotor bars were inextricably 
entwined and could not be separated in the analysis. 

5.2. Machine description 

An attempt was made to design the same rpm and pole- 
pair machines as was performed for the synchronous and per- 
manent magnet models, but problems in limiting rotor current 
density allowed only a few of the machines to be designed. 
For example, no 7200 rpm machines were designed and 180 rpm 
machines could only be designed with up to twelve poles. 

All of the induction machines are listed in the tables 
starting on page . Only medium confidence should be placed 
in the induction machine designs, as there were some conver- 
gence difficulties in slip. (Slip is not listed in the 
tables for that reason. ) 

5.2.1. Efficiency 

The higher rpm machines showed a slight increase in ef- 
ficiency as rpm increased. There was much movement around 
the average value of 97.5%. The movement decreased as rpm 
increased. With only six machines, it is hard to detect a 
trend in 180 rpm machine efficiency. Apparently, efficiency 
did increase with the number of pole-pairs, with all ef- 
ficiencies below 90%. Developmental studies for this thesis 
showed that of f-design-point efficiencies for the 180 rpm 
machines were sometimes below 70% for the endurance speed 
condition. 
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EFFICIENCr 



Figure 14. High rpm induction machine efficiency 
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EFFICIENCY 



Figure 15. 180 rpm induction machine efficiency 



INDUCTION MOTOR EFFICIENCY 



DIRECT-DRIVE, 180 RPM. 25775 HP 
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VOLUME (M^3) 



Figure 16. High rpm induction machine volume 



INDUCTION MOTOR VOLUME 



25775 HP 
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VOLUME (M^3) 



Figure 17. 180 rpm induction machine volume 



INDUCTION MOTOR VOLUME 



DIRECT-DRIVE, 180 RPM. 25775 HP 
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Figure 18. High rpm induction machine weight 
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Figure 19. 180 rpm induction machine weight 



INDUCTION MOTOR WEIGHT 



DIRECT-DRIVE, 180 RPM. 25775 HP 
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5.2.2. Weight and volume 



The best of the high rpm machines rivaled the 
synchronous machines in weight and volume. The worst were 
very bad. Weight and volume decreased with an increase in 
the number of poles, but not necessarily with the increases 
in rpm. For the 180 rpm machines, both weight and volume 
decreased dramatically as the number of poles went from two 
to six, with much lower decreases after that. The 180 rpm 
machines were uncompetitive in the synthesis process. 

5.3. Machine characteristics 

The previously mentioned rotor current density dif- 
ficulty showed in the rotor slot factor, which was at the 
limit of 0.75 for almost every motor. The stator slot fac- 
tor gradually grew with the increase in poles, arriving at 
0.75. The length and diameter of both the 180 and higher 
rpm machines is such that they would fit in machinery 
spaces . 

5.4. Verification 

Induction machines were expected to be close to 
synchronous machines in volume, weight, and efficiency. 

They were, and this comparison served as the verification 
for the induction machine model. Because the confidence 
level in the designed machines is only medium, more work 
would be needed to verify that these machines would have the 
advertised properties if built. 
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Table 18. Characteristics of 1800 rpm induction machines 



nuiber of pole pairs 
poHer, hp 
efficiency factor 
voluie factor 
shaft rpi 
priiary atp-turns 
synchronous frequency 
rotor radius 
gap diiension 
back iron depth 
stator slot depth 
rotor slot depth 
stator slot factor 
rotor slot factor 
envelope voluie 
envelope weight 
hysteresis loss 
eddy current loss 
stator copper loss 
full load efficiency 
active length 
rotor copper loss 
iaxiiui torque 
teriinal volts/turn 
air gap volts/turn 
Rl/turns-squared 
Xl/turns“squared 
Xi/turns-squared 
X2/turns-squared 
R2/turns*squared 
overall length 
overall diateter 



1 


2 


3 


4 


5 


6 


25775 


25775 


25775 


25775 


25775 


25775 


IBOOOO 


180000 


180000 


180000 


180000 


180000 


1286. 1 


1286.1 


1286.1 


1286.1 


1286.1 


1286.1 


1800 


1800 


1800 


1800 


1800 


1800 


6.10E+05 


7.74E+05 


2.44E+05 


1.91E+05 


1.86E+05 


1.79E+05 


30 


60 


90 


120 


150 


180 


0.9236 


1.0197 


0.4638 


0.4352 


0.4060 


0.3945 


0.0262 


0.0039 


0.0036 


0.0029 


0.0049 


0.0051 


0.1894 


0.0173 


0.1073 


0.0093 


0.0405 


0.0460 


2.1467 


0.5423 


0.3375 


0.3490 


0.1902 


0.1388 


0.2051 


0.3399 


0.1546 


0.1451 


0.1353 


0.1176 


0.285 


0.641 


0.355 


0.640 


0.699 


0.749 


0.750 


0.750 


0.750 


0.750 


0.750 


0.748 


148.863 


36.240 


7.558 


5.880 


3.760 


3.164 


376603.7 


80669.2 


29447.0 


21499.7 


15333.1 


13395.9 


180951.2 


19122.1 


68838.7 


39013.4 


42456.7 


48214.8 


32425.7 


6853.2 


37006.8 


27964.1 


38040.3 


51839.3 


1.03E+06 


1.45E+06 


498668.3 


165391.9 


264414.2 


320360.8 


0.939 


0.928 


0.969 


0.987 


0.981 


0.977 


0.191 


0.090 


0.759 


0.930 


1.000 


1.081 


8427.7 


8543.0 


19627.2 


22347.5 


30471.4 


37148.9 


185470 


185470 


185470 


185470 


185470 


185470 


234.31 


45.07 


266.45 


233.56 


193.76 


189.53 


63.29 


32.92 


126.58 


145.45 


145.98 


153.34 


9.23E-07 


8.09E-07 


2.80E-06 


1.50E-06 


2.54E-06 


3.32E-06 


3.70E-04 


3.89E-05 


9.60E-04 


9.53E-04 


6.82E-04 


6.17E-04 


2.51E-03 


4.41E-03 


1.23E-02 


1.28E-02 


6.23E-03 


5.15E-03 


2.03E-05 


3.66E-05 


3.76E-04 


6.52E-04 


7.79E-04 


9.08E-04 


4.50E-08 


9.97E-09 


3.89E-07 


5.37E-07 


6.62E-07 


8.49E-07 


3.990 


4.184 


2.629 


2.682 


2.644 


2.680 


6.572 


3.166 


1.824 


1.593 


1.283 


1.169 
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Table 19. Characteristics of 2400 rpm induction machines 



number of pole pairs 


2 


3 


4 


5 


6 


power, hp 


25775 


25775 


25775 


25775 


25775 


efficiency factor 


180000 


180000 


180000 


180000 


180000 


voluie factor 


1286. 1 


1286.1 


1286.1 


1286.1 


1286.1 


shaft rpi 


2400 


2400 


2400 


2400 


2400 


priiary aip-turns 


5.94E+05 


3.78E+05 


2.45E+05 


1.78E+05 


1.93E+05 


synchronous frequency 


80 


120 


160 


200 


240 


rotor radius 


0.7541 


0.7958 


0.4803 


0.4126 


0.3841 


gap diiension 


0.0069 


0.0020 


0.0094 


0.0087 


0.0081 


back iron depth 


0.0952 


0.0013 


0.0834 


0.0575 


0.0448 


stator slot depth 


0.6169 


0.5475 


0.2755 


0.2341 


0.1561 


rotor slot depth 


0.2514 


0.1979 


0.1601 


0.1375 


0.1197 


stator slot factor 


0.454 


0.477 


0.417 


0.413 


0.750 


rotor slot factor 


0.750 


0.750 


0.639 


0.568 


0.750 


envelope voluie 


23.889 


21.180 


6.201 


4.554 


2.942 


envelope weight 


60526. 1 


30880.0 


19931.3 


17243.1 


11459.0 


hysteresis loss 


78227.5 


26111.4 


80426.4 


86278.8 


54014.3 


eddy current loss 


37381.4 


18716.2 


76864.2 


103071.8 


77433.0 


stator copper loss 


l.lOE+06 


334484.0 


376335.4 


306069.2 


282567.9 


full load efficiency 


0.940 


0.980 


0.972 


0.974 


0.978 


active length 


0.142 


0.189 


0.533 


0.909 


0.852 


rotor copper loss 


8819.4 


7426.0 


17135.5 


21589.8 


27954.3 


laxiiui torque 


139102.5 


139102.5 


139102.5 


139102.5 


139102.5 


teriinal volts/turn 


103.23 


113.69 


229.60 


346.92 


212.59 


air gap volts/turn 


51.27 


72.14 


122.69 


179.75 


156.78 


Rl/turns-squared 


1.04E-06 


7.82E-07 


2.08E-06 


3.23E-06 


2.52E-06 


Xl/turns-squared 


1.50E-04 


2.32E-04 


7.89E-04 


1.67E-03 


7.40E-04 


Xi/turns-squared 


3.88E-03 


1.25E-02 


3.37E-03 


4.31E-03 


3.35E-03 


X2/turns-5quared 


6.22E-05 


1.56E-04 


4.63E-04 


1.07E-03 


9.36E-04 


R2/turns~squared 


2.84E-08 


4.44E-08 


2.97E-07 


7.69E-07 


6.73E-07 


overall length 


3.186 


3.380 


2.492 


2.594 


2.420 


overall diaieter 


2.946 


2.693 


1.697 


1.426 


1.186 
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Table 20. Characteristics of 3000 rpm induction machines 



nuiber of pole pairs 


2 


3 


4 


5 


6 


power, hp 


25775 


25775 


25775 


25775 


25775 


efficiency factor 


IBOOOO 


180000 


180000 


180000 


180000 


voluie factor 


1286.1 


1286.1 


1286.1 


1286.1 


1286.1 


shaft rpi 


3000 


3000 


3000 


3000 


3000 


primary aip-turns 


3. lOE+05 


3.62E+05 2.80E+05 


1.76E+05 


1.37E+05 


synchronous frequency 


100 


150 


200 


250 


300 


rotor radius 


0.5198 


0.6366 


0.5580 


0.3660 


0.3325 


gap disension 


0.0138 


0.0036 


0.0070 


0.0070 


0.0039 


back iron depth 


0.1739 


0.0134 


0.0178 


0.0503 


0.0388 


stator slot depth 


0.8926 


0.3785 


0.3668 


0.2274 


0.1288 


rotor slot depth 


0.1733 


0.2122 


0.1860 


0.1220 


0.1093 


stator slot factor 


0.340 


0.607 


0.545 


0.476 


0.721 


rotor slot factor 


0.750 


0.750 


0.750 


0.750 


0.750 


envelope volute 


22.412 


10.122 


8.000 


3.379 


2.065 


envelope weight 


72906.0 


19541.1 


15961.5 


13000.9 


9016.0 


hysteresis loss 


185978.6 


25156.1 


44235.6 


74610.2 


52848.0 


eddy current loss 


111088.5 


22539.3 


52845.5 


111415.1 


94701.2 


stator copper loss 


322718.0 


359218.8 221369.1 


269959.2 


226785.8 


full load efficiency 


0.968 


0.979 


0.983 


0.976 


0.980 


active length 


0.383 


0.189 


0.307 


0.818 


1.006 


rotor copper loss 


8574.1 


9089. 1 


9006.9 


18292.9 


25568. 1 


laxiiui torque 


111282 


111282 


111282 


111282 


111282 


tertinal volts/turn 


421.67 


101.44 


174.86 


361.45 


251.07 


air gap volts/turn 


119.20 


72.02 


102.80 


179.30 


200.52 


Rl/turns-squared 


1.12E-06 


9.14E-07 9.39E-07 


2.91E-06 


4.01E-06 


Xl/turns-squared 


1.31E-03 


1.96E-04 5.04E-04 


1.78E-03 


1.09E-03 


Xi/turns“squared 


4.47E-03 


6.93E-03 3.82E-03 


5.33E-03 


9.00E-03 


X2/turns-squared 


1.B2E-04 


1.72E-04 3.25E-04 


9.97E-04 


1.54E-03 


R2/turns-squared 


1.60E-07 


5.15E-08 


1.08E-07 


6.66E-07 


l.OlE-06 


overall length 


2.517 


2.750 


2.567 


2.310 


2.352 


overall diaieter 


3.210 


2.064 


1.899 


1.301 


1.008 
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Table 21. Characteristics of 3600 rpm induction machines 



nusber of pole pairs 
power , hp 
efficiency factor 
volute factor 
shaft rpa 
priiary atp-turns 
synchronous frequency 
rotor radius 
gap diiension 
back iron depth 
stator slot depth 
rotor slot depth 
stator slot factor 
rotor slot factor 
envelope volute 
envelope weight 
hysteresis loss 
eddy current loss 
stator copper loss 
full load efficiency 
active length 
rotor copper loss 
•axiiui torque 
tertinal volts/turn 
air gap volts/turn 
Rl/turns-squared 
Xl/turns-squared 
Xt/turns-squared 
X2/turns-squared 
R2/turns-squared 
overall length 
overall diaeeter 



1 


2 


3 


25775 


25775 


25775 


180000 


180000 


180000 


1286.1 


1286.1 


1286.1 


3600 


3600 


3600 


2.72E+05 


1.70E+05 


2.99E+05 


60 


120 


180 


0.4964 


0.4350 


0.5150 


0.0061 


0.0027 


0.0033 


0.3428 


0.0133 


0.0238 


1.3891 


1.0052 


0.3606 


0.1654 


0.1308 


0.1717 


0.290 


0.338 


0.585 


0.750 


0.750 


0.750 


41.461 


17.960 


6.532 


164461.5 


54127.3 


14922.3 


261300.3 


145818.0 


34381.6 


93647.7 


104519.7 


36966.2 


341655.1 


107764.3 


288857.8 


0.965 


0.981 


0.981 


0.393 


0.700 


0.247 


7562.3 


7220.5 


10787.1 


92735 


92735 


92735 


475.12 


702.61 


143.82 


140.25 


218.89 


91.39 


1.54E-06 


1.25E-06 


1.07E-06 


1.67E-03 


3.93E-03 


3.70E-04 


2.39E-02 


4.19E-02 


9.66E-03 


1.20E-04 


4.65E-04 


2.54E-04 


1.98E-07 


4.57E-07 


1.03E-07 


2.403 


2.451 


2.320 


4.469 


2.913 


1.805 



4 


5 


6 


25775 


25775 


25775 


180000 


180000 


180000 


1286.1 


1236.1 


1286.1 


3600 


3600 


3600 


2.95E+05 


2.05E+05 


1.23E+05 


240 


300 


360 


0.5305 


0.4338 


0.3144 


0.0088 


0.0079 


0.0030 


0.0238 


0.0258 


0.0367 


0.3236 


0.2792 


0.1181 


0. 1768 


0.1446 


0.1045 


0.688 


0.553 


0.750 


0.750 


0.750 


0.749 


6.536 


4.331 


1.722 


14074.1 


11549.4 


7593.8 


42814.0 


62801.7 


52473.1 


61376.6 


112537.8 


112835.4 


210564.6 


173450.1 


192320.6 


0.983 


0.982 


0.981 


0.249 


0.480 


0.964 


9954.8 


11162.4 


23436.6 


92735 


92735 


92735 


142.32 


275.51 


261.98 


94.86 


149.83 


218.01 


8.08E-07 


1.37E-06 


4.25E-06 


3.59E-04 


1.13E-03 


1.18E-03 


2.80E-03 


3.94E-03 


1.26E-02 


3.05E-04 


7.24E-04 


1.84E-03 


9.68E-08 


2.79E-07 


1.07E-06 


2.406 


2.247 


2.234 


1.773 


1.494 


0.944 
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Table 22 . Characteristics of 180 rpm induction machines 



nusber of pole pairs 
power, hp 
efficiency factor 
voluie factor 
shaft rpi 
priiary aips-turns 
synchronous frequency 
rotor radius 
gap diiension 
back iron depth 
stator slot depth 
rotor slot depth 
stator slot factor 
rotor slot factor 
envelope voluae 
envelope weight 
hysteresis loss 
eddy current loss 
stator copper loss 
full load efficiency 
active length 
rotor copper loss 
aaxiaui torque 
terainal volts/turn 
air gap volts/turn 
Rl/turns-squared 
Xl/turns-squared 
Xt/turns-squared 
X2/turns“squared 
R2/turns-squared 
overall length 
overall diaieter 



1 


0 

i. 


3 


25775 


25775 


25775 


180000 


180000 


180000 


128&.1 


1286.1 


1286.1 


180 


180 


180 


1054077 


564571.5 


569333.6 


3 


6 


9 


0.995i 


0.7236 


0.6896 


0.0025 


0.0060 


0.0195 


0.0227 


0.0240 


0.1568 


0.i714 


0.7116 


0.2536 


0.3319 


0.2412 


0.2299 


0.590 


0.563 


0.734 


0.750 


0.750 


0.750 


47.&&8 


36.837 


20.316 


217524.1 


168344.2 


98690.0 


4551.5 


16194.6 


19356.0 


81. & 


580.4 


1040.6 


4.44E+06 


1.96E+06 


3.05E+06 


0.808 


0.901 


0.854 


0.825 


2.047 


1.854 


1.08E+05 


1.29E+05 


2.08E+05 


1854700 


1854700 


1854700 


61.849 


104.401 


58.731 


29.510 


53.243 


45.967 


1.33E-06 


2.05E-06 


3.14E-06 


5.08E-05 


1.58E-04 


6.01E-05 


1.21E-02 


4.61E-03 


8.17E-04 


1.89E-05 


6.91E-05 


7.91E-05 


1.16E-07 


4.44E-07 


4.33E-07 


4.817 


4.965 


4.691 


3.384 


2.930 


2.239 



4 


5 


6 


25775 


25775 


25775 


180000 


180000 


180000 


1286.1 


1286.1 


1286.1 


180 


180 


180 


510008.9 


512534.4 


324264.5 


12 


15 


18 


0.6616 


0.7411 


0.6962 


0.0334 


0.0235 


0.0054 


0.1158 


0.1038 


0.0812 


0.2231 


0.2200 


0.1466 


0.2039 


0.1845 


0.1501 


0.750 


0.694 


0.718 


0.683 


0.696 


0.701 


19.450 


21.357 


17.600 


94554.9 


93084.9 


87558.9 


24294.6 


31609.4 


36548.6 


1741.4 


2832. 1 


3929.6 


3.15E+06 


2.59E+06 


2.09E+06 


0.850 


0.872 


0.891 


2.486 


2.159 


3.090 


2.23E+05 


1.94E+05 


2.14E+05 


1854700 


1854700 


1854700 


76.530 


76.097 


87.418 


59.124 


57.510 


77.345 


4.04E-06 


3.28E-06 


6.61E-06 


9.01E-05 


9.33E-05 


1.12E-04 


4.60E-04 


5.08E-04 


2.49E-03 


1.37E-04 


1.27E-04 


2.40E-04 


7.40E-07 


6.20E-07 


1.15E-06 


5.266 


5.217 


5.897 


2.068 


2.177 


1.859 
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Chapter Six. Nominal ship design 



6.1. Technology sensitivity analysis 

Technology sensitivity analyses, such as this thesis, 
must be able to quantitatively compare similar technologies. 
The pertinent differences must be made apparent through ap- 
propriate analysis. Inherent in the analysis must be the 
consideration of the global system complexity. Naval ships 
are extremely complex and the effects of various tech- 
nologies can be lost in the complexity. One methodology for 
technology characterization in naval ships has been proposed 
by Goddard [25]. This method has been followed to show the 
benefits of electric drive. 

6.2. ASSET 

The Advanced Surface Ship Evaluation Tool (ASSET) was 
developed over several years to be the U. S. Navy’s premier 
ship design computer program. It has its roots in HANDE, a 
hydrofoil design program developed in the 1970’ s. The ship 
design spiral is traversed in an iterative fashion until 
convergence on a number of parameters is achieved. Boeing 
Computer Services is the contractor for ASSET, under the su- 
pervision of the David Taylor Naval Ship Research and 
Development Center, at Carderock, Maryland.!® it is divided 
into large modules by ship type. These modules include 
Monohulls, Hydrofoils, and SWATH ships. 

The geometry of a particular ship is input to ASSET. 

The following were used as characteristic ship traits: full 

load displacement, certain Ship Work Breakdown Structure 
(SWBS) weight groups,! ^ endurance fuel load at 20 knots, 

16. Greenwood, R.W. , and Fuller, A.L., "Development of a 
Common Tool for Ship Design and Technology Evaluation," 

Proc. SNAME New England Section Marine Computers 1986. 

17. See Appendix E for more information on SWBS. 
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draft, maximum and sustained speed for 51,550 Installed 
horsepower, and transverse GM. From the changes In these 
parameters during the various computer runs , the effects of 
transmission choices were noted. 

The Inputs to ASSET describe the ship that Is being 
designed. The outputs cover the range of calculations pos- 
sible In structures, volume, space, machinery, propeller 
characteristics, resistance, powering, and weight. There 
are ASSET performance modules on cost, stability, hydros- 
tatics, seakeeping, manning and space but the usual syn- 
thesis output Is of more use during a technology sensitivity 
analysis . 

The descriptions of several ships are contained in an 
ASSET data base. For a particular ship, a Current Model is 
maintained that holds all of the parameters to describe that 
particular ship. In ASSET Version 2.0, over 380 parameters 
are used for each ship description. User control over most 
of these parameters is possible, or control may be given to 
the executive program which will then "design" a ship sub- 
ject to whatever constraints the user desires. 

Some intricacies of ship design are not handled well by 
ASSET. For example, the program is not able to handle 
equipment re-arrangements easily, and almost all equipment- 
level volumes are approximated from studies of past ships. i® 
For this reason, some equipment-level weights and volumes 
need to be calculated off-line and input to the program 
through its weight adjustment facility, especially if ac- 
curacy in these areas is important to the study being per- 
formed. 

The baseline ship used in these studies has a full load 
displacement of 5485 LT, carries 272 crew members, is 425 
feet long, and has a primary mission of anti-submarine war- 

18. The Enhanced Machinery Module [27], in the process of 
being made available, will improve the situation dramati- 
cally. Some of the relationships from that module were used 
in calculating electric propulsion weights. 



84 



fare (ASW) . It is armed for that purpose and has equipment 
in keeping with its size. The baseline ship is described in 
more detail below. 

6.2.1. Margins 

A naval ship design has margins in weight, vertical 
center of gravity (KG), space, ship service electrical gen- 
eration, propulsion power, accommodations, and structural 
strength which allow for equipment, mission, and system 
growth over its projected thirty year life. Without these 
margins, the ship would be difficult to modify, because 
whatever might be added in these areas would have to be paid 
for by a removal. For example, if a 50 ton radar system 
were added, the original 40 ton radar system and 10 tons of 
fuel might be removed to leave the ship at its original 
weight. With margins, the 10 tons of fuel might not be 
removed . 

Margins are typically split into Acquisition and Serv- 
ice Life allowances. Acquisition margins recognize the fact 
that ship specifications change over the design cycle and 
during construction. For example, the fourth ship built may 
have a different weapons system than the first, with a dif- 
ferent electrical requirement. If the electrical generation 
plant had to be changed during construction to accommodate 
the new weapons system, the total cost might be prohibitive. 
If an Acquisition margin is built into the original design, 
this may not occur. A Service Life margin makes allowance 
for configuration changes over the life of the ship. 

The ASSET program uses margins when synthesizing a 
ship. The margins are under operator control. The margins 
suggested by Goddard and used in this analysis are listed in 
the table. 1® 

19. Goddard, C.H., “A Methodology for Technology Charac- 
terization and Evaluation for Naval Ships," S.M. NAME AND 
O.E. thesis, MIT, 1985, p. 31. 
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Table 23. Recommended technology assessment design margins 

for a monohull surface combatant 



Acquisition Service Life 



Weights 

KG 

Space 

Electrical^ 
Propulsion powers 

Accommodations 

Strength 



12.5% of SWBS 1-7 10% 

12.5% of KG 1-7 1.0 ft. 

0 (no excess volume) 0 

20 % 20 % 



10% total EHP prior to prelim body plan 
8% prior to self-propelled model tests 
Accom = 1.1 X ship manning at delivery 
2.24 KSI of marginal stress at delivery 
(Max primary stress for hull material) 



Notes : 

a. The service life weight margin applies only to naval 
architectural limits of the ship (reserve buoyancy, 
stability, structures), not to the final design weight. 

b. In sizing the electric plant, the calculated maximum 
electric load plus these design margins shall be met with 
one generator out of service. The remaining generators 
shall not be loaded in excess of 90%. Note that the service 
life margin is not applied to SWBS group 200 which would be 
expected to remain stable over the life of the ship. 

c. Performance requirements (sustained speed, endurance 
range) are met at delivery full load displacement. 



6.3. Philosophy of effort 

The nature of this technology characterization required 
that certain limits be imposed on the total effort. (If the 
Naval Sea Systems Command were to do this study, many people 
would simultaneously be employed to investigate every 
detail.) Some items were fixed, some were allowed to float 
with the design. 

The hardest item to handle is volume. There are very 
few ways of adjusting volume as easily as weight is ad- 
justed. One way is through the use of Marginal Volume Fac- 
tors, which equate a weight penalty with every increase in 
volume. (See Howell [26].) 

The differences among transmission systems appear 
primarily in the machinery spaces and fuel tanks of a ship. 
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ASSET can handle the tankage changes without aid. However, 
it does no equipment arrangement or space analysis inside 
the machinery spaces, leaving that for arrangement experts 
to do off-line. It was decided to keep the machinery space 
volume the same in all ships of this study, no matter how 
that volume was divided into spaces. Changes in equipment 
volume will be noted in the analysis and left for the advan- 
tage of others in machinery space rearrangement. 

Moderate to high technical risk has been accepted in 
specifying equipment cooling. Current densities at the 
limits of cooling technology are used, assuming that liquid 
cooling of both the stator and rotor can be performed. 

Lower current densities would be required if such cooling 
were not possible, resulting in slightly larger, heavier, 
less efficient machines. This last statement was proven 
during the course of the thesis research, as the first 
(chronologically) current densities used were two-thirds 
larger than those listed herein. 

Other risk areas include the use of advanced vacuum 
switchgear and the assumed efficiencies of reduction gears 
and power converters. These are low risk items; the 
technology is well understood and commercially available. 



Table 24. Ship design items held constant during analysis 



Endurance speed 
Endurance ranged 
Machinery box volume 
Installed horsepower 

Payload weight, volume, and electrical 

Length 

Becun 

Ship electrical load 
Ship molded lines 
Manning 

Deckhouse and superstructure geometry 



20 knots 

5500 nautical miles 
109,670 ft3 
51,550 hp 
requirements 
425 ft. 

55 ft. 

2030 kW (24 hr. avg) 
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Table 25. Ship design items allowed to float during analysis 



Maximum and sustained speeds 
Endurance fuel load 
Ship full load displacement 
Ship draft 

Ship resistance and powering 

Ship arrangement outside machinery spaces 



Note: 

a. Endurance range was allowed to float as a comparison 
between two electrical transmission ships and the mechanical 
transmission baseline ship. It was held constant for the 
rearranged ship. 
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Most of the ship synthesis has been left to ASSET 
though parameters from Goddard were used where possible. 
ASSET designs a reasonable, generic ship with good seakeep- 
ing characteristics. As shown above, some ship characteris- 
tics were frozen to ensure transmission comparisons were not 
performed with different ships. Leaving the ship synthesis 
to ASSET allowed concentration on the specifics of the 
propulsion plant. 

6.4. Baseline ship 

The baseline ship has a mechanical transmission, i.e., 
two power trains consisting of a gas turbine, clutch and 
coupling, reduction gear, shafting, and propeller. There 
are two machinery rooms, each containing one gas turbine. 

The gas turbine used as the model is the General Electric 
LM-2500, rated at 25,775 brake horsepower. This is a very 
common marine gas turbine. It and its predecessors are 
powering the latest classes of naval combatant, such as the 
DD963, FFG7, and DDG51. 

The locked-train double-reduction gears are reversible, 
allowing the use of fixed-pitch propellers. There is no 
mechanical cross-connect allowed between the shafts. Except 
for the power level, this is the gear system being employed 
on the DDG51 class. Gears of this sort are about one per- 
cent ineff icient20 per reduction stage. Since these are 
double-reduction gears, an efficiency of 98% was used. 

An endurance speed of 20 knots has been specified. 

This is in keeping with standard fleet practices. An en- 
durance range of 5500 nautical miles permits ocean crossings 
without refueling. The lack of a cross-connect capability 
between the two shafts means at least two gas turbines will 



20. Inefficiency = 1 - efficiency. Information on stage in- 
efficiency is from a conversation with Mr. Samuel Shank, the 
author of the ASSET Enhanced Machinery Module [27]. 
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be on-line during endurance cruising. This is inherently 
inefficient; some operators alleviate the inefficiency by 
declutching one shaft and free-wheeling that propeller, 
reaching desired speeds by loading up the operating gas tur- 
bine . 

A measure of initial static stability is the ratio 
GM/Beam. GM is the vertical distance between the center of 
gravity and the metacenter of the ship. Typical values for 
this ratio are 8-10%. A lower value (6.5%) has been ac- 
cepted for purposes of comparison with the variants. A 
large ship redesign effort would have been necessary to 
bring GM/B into a better range. 

The electrical generation plant consists of three gas 
turbines, each driving a 2000 kW generator. The data used 
for the gas turbines was taken from the Detroit Allison 50 IK 
turbine-generator set used aboard the DD963 ships. 

Both the deckhouse and main hull are constructed of 
High Tensile Strength (yield stress = 50,000 psi) steel. 
Active stabilizing fins and a sonar dome are included in the 
design. The payload is listed in Table 26. A coarse layout 
of the machinery spaces is shown in Figure 20. 



Table 26. Payload for baseline and variant ships 

FFG7 Command and control suite 
Satellite, UHF, and HF communications 
SLQ32V3 electronic countermeasures 
NIXIE acoustic countermeasure 

SPS-49 two-dimensional air search and tracking radar 

SPS-55 surface search and tracking radar 

SQR-19 towed array surveillance system 

MK92 missile and gun fire control system 

Harpoon fire control system 

LAMPS III helicopters and support system 

JP-5 aviation fuel 

MK32 over-the-side torpedo system 

MK13 guided missile launching system 

MK75 76mm gun 

Close- In-Weapon-System 

Small arms 

Appropriate ammunition and reloads 



90 



Figure 20. Baseline ship subdivision and machinery arrange- 
ment 
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6.5. Backfit ship 



If the mechanical transmission of the baseline ship 
were replaced with an electrical transmission without chang- 
ing the subdivision of the ship, it would be as if an older 
ship had been updated, or "backf itted, ” with new technology. 
This is the idea behind the first variant ship. 

Propulsion motors and generators were added into the 
propulsion plant, in the original spaces. Shafting still 
runs from the machinery spaces to the propellers. A reduc- 
tion gear is still necessary for higher rpm propulsion 
motors, but electrical cross-connect may improve the en- 
durance fuel efficiency. 

Some rearrangement within the machinery space is neces- 
sitated by the backfit. The machinery spaces in the 
baseline ship are not long enough to contain the stack-up 
length of gas turbine, generator, motor, and reduction gear, 
without “folding" the power train. This may be accomplished 
by changing the design of the reduction gear or by placing 
the gas turbine and generator (which require a mechanical 
connection) side-by-side (or transversely) with the motor 
and reduction gear (another mechanical connection), using 
transmission line to electrically connect them. Since there 
are a variety of ways to rearrange the machinery box, and 
since the chosen method has no effect on the analysis, the 
rearrangement was not specified. 

The propulsion generators and motors may operate at 
different rotational speeds, and therefore different elec- 
trical frequencies, so power converters must be used between 
them. Power converters change the frequency of the power 
being transferred between the generators and motors , through 
the use of cycloconverters or thyristors. They add another 
inefficiency to the transmission. A reasonable estimate of 
the efficiency of an 18 MW power converter is 97 %. 21 

21. Professor John Kassakian, MIT, private communication. 
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6.6. Rearranged ship 



The second variant ship takes advantage of the benefits 
of machinery space rearrangement . The same machinery space 
total volume is preserved but split into five spaces. The 
propulsion motors are placed very near the propellers, 
resulting in much shorter shafting runs and an increased 
GM/B ratio. The decrease in shaft length means a decrease 
in shaft weight and more space available outside the 
machinery box. The shafts previously ran through shaft al- 
leys that may be returned to other uses. This rearrangement 
is almost certainly not the optimum one and can be improved 
in the sense of space efficiency. It is an arrangement, 
however, that can demonstrate the benefits to be expected of 
a ship designed for electric drive. Power converters again 
connect the generators and motors. Transmission line forms 
the connections, at a much lower weight than shafting. 

Many other choices in large components are possible for 
this rearranged ship. For example, three propulsion gener- 
ators and gas turbines driving two or four propulsion motors 
might have been chosen. The number of prime movers was 
retained from the mechanical baseline ship, however, to make 
the comparison of transmissions realistic. Too many chauiges 
might have obscured the fundamental differences in 
efficiency, weight and volume. 

Both geared and direct-drive propulsion motors were 
used in this variant. When geared motors are used, the 
reduction gears are also placed near the propellers. A 
coarse layout of the rearranged ship is shown in Figure 21. 

6.7. Weight and volume algorithms 

Few components in this thesis are exact commercial 
models. The weight, volume, and other characteristics are 
taken from those for which data was available. The equa- 
tions for shafting and reduction gears were taken from the 
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ASSET theory manuals [28], while the switchgear and braking 
resistor equations came from the ASSET Enhanced Machinery 
Module [27], which is not yet generally available. These 
equations represent much study by ship and equipment desig- 
ners, incorporating equipment which commercially avail- 
able. Where possible, the ASSET equations were verified 
against other studies and actual equipment [29, 30, 31]. 

For example, the machinery in the FFG7 was used as a model 
and verification for reduction gears and shafting. 

Appendix E gives more explanation, as well as presenting a 
computer program used to generate weight and volume figures. 
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Figure 21. Layout of machinery spaces on rearranged ship 
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Chapter Seven. Machine design and system synthesis 



7.1. Machine matrix 

During this study, synchronous, permanent magnet, and 
induction machines were designed at shaft speeds of 180, 
1800, 2400, 3000, 3600, and 7200 rpm. The number of pole- 
pairs for the 180 rpm machines ran from one to twenty-five. 
For the higher rpm machines, from one to six pole-pairs were 
used. If every machine could be a generator or a motor, a 
165x165 machine matrix results. From these 27,225 possible 
combinations, two were chosen and input to the two variant 
ships for synthesis in ASSET. 

These particular rpms were chosen partially because of 
the choice of the gas turbine. The LM-2500 operates at a 
full-load speed of 3600 rpm, making multiples and "nice" 
fractions of that speed desirable. A 3600 rpm, two-pole 
machine has a synchronous frequency of 60 hertz, the stand- 
ard in the United States. A four-pole machine at 1800 rpm 
is also a 60 hertz machine. If "nice" frequencies result 
from rpm choices, results may not be obscured. The low rpm 
machine is tied to propeller rpm. For the baseline ship, 
maximum speed propeller rpm is 170 rpm. A ten rpm allowance 
for "battle override" gives a requirement for 180 rpm. 

The reduction in the number of machine combinations is 
a bit more difficult to explain. First, it was observed 
that reduction gears add greatly to the weight and volume of 
the transmission and detract from its efficiency. Second, 
induction generators are notoriously difficult to control. 

It was then decided that generators would not be induction 
machines and any generators used would operate at the same 
shaft speed as the gas turbine, eliminating a possible 
reduction gear. The matrix then measured 12x165 and had 
1980 combinations. 

From this point on, the decision theory espoused in 
Schweppe and Merrill [13] was used, specifically using "knee 
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sets" to eliminate uncompetitive designs. The considera- 
tion, or Figure-of-Merit , was always minimum weight, volume, 
and inefficiency. These external characteristics are those 
"seen" by a ship design (since propulsion current and volt- 
age were specified in terms of turns -- fertile ground for 
another tradeoff study). Initially, individual motors of 
all of the higher rpms were considered together to select 
the best geared propulsion motors. Since the weight and 
volume of a reduction gear varies with shaft rpm, the weight 
and volume of the reduction gear was included with that of 
the machine to select the best machine-gear system. Direct- 
drive machines were also selected. The 3600 rpm synchronous 
and permanent magnet machines were also considered 
separately as propulsion generators. 

Once the initial selection of machines was made, the 
matrix measured 5x11, or 55 combinations. All of the com- 
binations were plotted in knee curves that showed the 
volume, weight, and inefficiency of each transmission. The 
volume and weight of shafting, braking resistors, cooling 
systems, power converters, and the inefficiency of the power 
converters were common to all combinations and were not in- 
cluded at this level. The inefficiency of any reduction 
gears was included where appropriate. Three of the gener- 
ators were synchronous machines and two were permanent mag- 
net machines. Of the motors, two were 1800 rpm synchronous 
machines, two were 1800 rpm permanent magnet machines, three 
were 1800 rpm induction machines, and two were 180 rpm 
synchronous machines. The 1800 rpm motors clearly dominated 
the higher rpm machines, largely because of the differences 
in reduction gear weight and volume. 

Since there was no single dominating combination, ten 
of the 55 combinations were selected. These ten were among 
the best at least twice on the knee curves. These ten com- 
binations were composed of only synchronous machines . 

Programs to calculate off-design-point efficiency were 
written. Each motor and generator was evaluated at the 
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power level and rpm appropriate for the sustained and en- 
durance speed conditions of the ship, using the delivered 
horsepower (DHP) and propeller rpm taken from the ASSET syn- 
thesis run on the baseline mechanical ship. Some of the 
combinations had very low endurance efficiencies. There was 
a sharp division evident between the geared (lower weight, 
much higher volume) and the direct-drive systems. 

The ten combinations, with their maximum, sustained, 
and endurance speed transmission inefficiencies, were again 
made the subject of knee curves. A simple scoring scheme 
was devised to rank the combinations according to their 
grouping on these last knee curves. If a combination was in 
the best group on a particular knee curve, it was given two 
points for that curve. If it was in the second best group, 
it received one point. If it was in neither the best or 
second best group, it received no points. When the scores 
were totaled, two combinations stood out. One was a geared 
drive system and one a direct-drive system. These two com- 
binations were used in both of the variant ships and are the 
subject of the next chapter. 

7.2. Knee curves 

Figures 22 through 24 show the knee curves for the 
propulsion generators. The first letter of the generator 
ID indicates whether it is a synchronous machine or a per- 
manent magnet machine. All of the generators are 3600 rpm 
machines. The generators selected were SB (four poles), SC 
(six poles), SD (eight poles), PB (four poles), and PC (six 
poles). Generator SA was not selected because of its poor 
showing on the volume-weight curve, even though it was com- 
petitive on the volume-efficiency curve. 

The 180 rpm, direct-drive propulsion motor curves are 
in Figures 25 through 27. They were not combined with 
geared motors because one of the points being investigated 
was whether or not geared motors were "better" than direct- 
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drive motors. The clumping of the machines necessitated 
other graphs on different scales to distinguish between the 
machines. Machines 1-25 are synchronous, 26-50 are per- 
manent magnet, and 51-56 are induction machines. The number 
within a group indicates the number of poles in the machine, 
e.g., machine 32 has (32-25)x2=14 poles. Valid designs with 
over twelve poles were not achieved for the induction 
machines. Machines 3, 8, 11, and 16 were selected for fur- 
ther work. These are all synchronous motors. 

The higher rpm motors were combined with their reduc- 
tion gears to form system knee curves. In all of these knee 
curves, machines 1-6 are 1800 rpm, 7-12 are 2400 rpm, 13-18 
are 3000 rpm, 19-24 are 3600 rpm, and 25-30 are 7200 rpm. 
Figure 28 is the volume-efficiency curve for synchronous 
machines, showing the distinct grouping of the machine-gear 
systems due to the high volume of the reduction gears. Note 
the high values of the permanent magnet machines in the 
volume-weight curve Figure 29. The clumping of induction 
motors around the low inefficiencies is shown in the weight- 
efficiency curve of Figure 30. From these curves, the 
motors on page 97 were selected. 

The initial motor and generator combinations were made 
and plotted on more knee curves (Figures 31 to 33), On 
those curves, the high-volume or high-weight nature of the 
combinations can be seen. Since the multiple-attribute 
decision theory embodied in knee curves does not say how to 
select between high-volume or high-weight, the best of each 
were selected. Combinations 1, 2, 8, 9, 12, 13, 19, 20, 30, 
and 31 were chosen. The of f -design-point efficiencies were 
calculated and all of the information was plotted. 

Figure 34 is a bar-graph of the maximum, sustained, and en- 
durance speed transmission efficiencies of the various com- 
binations, including the reduction gears, if any, and power 
converters. The final combination knee curves are sum- 
marized in the scoring scheme of Table 27, which was ex- 
plained on page 98 . Combinations 12 and 20 were chosen to 
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use in the ships of the study. 

This is a good method to choose among the possible 
machines. During the course of this thesis, the above path 
was followed through several complete iterations and a few 
partial ones. As stated in Schweppe and Merrill, knee 
curves serve very well to eliminate uncompetitive options, 
allowing concentration on the better ones. 

Table 27. Final combination knee curve scores 



Combo ID 


Firsts Seconds 


Total 


1 


3 


1 


7 


2 


3 


0 


6 


8 


1 


0 


2 


9 


0 


1 


1 


12 


6 


1 


13 


13 


4 


0 


8 


19 


1 


1 


3 


20 


3 


4 


10 


31 


0 


2 


2 


Conclusion : 


test 


combinations 


12 and 
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Figure 22. Curve of volume-efficiency for 3600 rpm gener- 
ators 



PROPULSION GENERATORS 



25775 HP, 3600 RPM 
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Figure 23. Curve of volume -weight for 3600 rpm generators 



PROPULSION GENERATORS 



25775 HP, 3600 RPM 
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Figure 24. Curve of weight-ef f iciency for 3600 rpm gener- 
ators 



PROPULSION GENERATORS 



25775 HP. 3600 RPM 
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Figure 25. Curve of volume-efficiency for 180 rpm motors 



PROPULSION MOTORS 

DIRECT--DRIVE. 180 RPM. 26775 HP 
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Figure 26. Curve of volume-weight for 180 rpm motors 



PROPULSION MOTORS 



DIRECT-DRIVE, 180 RPM, 25775 HP 
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Figure 27. Curve of weight-efficiency for 180 rpm motors 



PROPULSION MOTORS 



DIRECT-DRIVE, 180 RPM, 25775 HP 
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Figure 28. Curve of volume-efficiency for geared motors 



SYSTEM KNEE CURVE 



SYNCHRONOUS MACHINES. 25775 HP 
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Figure 29. Curve of volume -weight for geared motors 
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Figure 30. Curve of weight-efficiency for geared motors 



SYSTEM KNEE CURVE 



INDUCTION MACHINES, 25775 HP 
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Figure 31. Volume-efficiency curve for initial PM and PG 
combinations 



INITIAL PM AND PG COMBINATIONS 



25775 HP 




110 



TOTAL WEIGHT (KG) 
(Thousands) 



Figure 32. Volume-weight curve for initial PM and PG com- 
binations 



80 -r 

75 - 

70 - 

65 — 

60 - 

55 - 

50 - 

45 - 

40 - 

35 - 

30 -- 
0 



INITIAL PM AND PG COMBINATIONS 

25775 HP 



44 



54 



41 

52 



42 

53 



“1 1 r 

20 40 

TOTAL VOLUME (M^3) 




"T" 

60 



111 



TOTAL ErnCIENCr) 



Figure 33. Weight-efficiency curve for initial PM and PG 
combinations 
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Figure 34. Final combination transmission efficiencies 



TRANSMISSION EFFICIENCY 



/\^ 
0,8 -/\^ 
/\yy 

/\i 




/\^ 



IS 

/\ 

/\ 

/\ 

/\ 



/\^ 

/\^ 



/\ 

/\ 

/\p 

/\p 

/\p 

/\p 

/\p 

/\p 

/\p 



PTS 

/\p 

/\p 

/\p 

/\p 

/\p 

/\p 

/\p 

/\p 

MnM 



PT^ 



IS 

kM^ 

/\yy 

/\y 

WNp 

14\^ klMkl 



/\p 
/\p 

/\p 

/\^ 



/\p 



/\y 

/\y 



I4\t^ ]A\p 



/\y /\y /\y 



kt>k 

/\k 

4 






\/tM 






/\p 

/\y 

/\^ /\;^ 



H\k 



. y 

/\y 

/w/y /\y 



/\y /\% 



/\ 
/\ 
/\ 

/\p 

/\p 

/y 



^ -y 






l/N:^ I4\^ 






/\ 

/\ 

/\ 

/\ 

/\ 

/\ 

/\Z 

/\Z 

/\z 

/\p 

/\z 

/\p 

/\z 

/\p 

/\z 

/\p 



z^S 

z\ 

z\ 



i4\Z 



z\z 

W\p 



kNp 



z\p 

Z\Z 

z\z 

z\z 

z\z 

Z\Z 

z\z 

z\z 

z\z 

Z\Z 

zffi 



1 

2 



I 

12 



1.3 



19 



T 

20 



30 



I 

31 



[Z71 MAX 



COMBINATION ID 

1X3 sus 



WZX end 



113 



Chapter Eight. Analysis 



This chapter presents an analysis of the ships with 
electric transmissions. Standard naval architectural 
methods have been used to observe and comment on the variant 
ships, comparing them to the mechanical baseline ship. Con- 
clusions and recommendations follow the analysis. 

The names used to describe the various ships imply 
their internal arrangement, equipment, and ASSET Design Mode 
Indicator (DMI). The two DMIs used were ENDURANCE, when en- 
durance range was held constant at 5500 NM, and FUEL WT, 
when the usable fuel weight was held constant at 996.6 LT.22 
The ship names are as follows: 



MECH 23 BASELINE: Baseline, mechauiical transmission 

ELEC 23 BASELINE 12: Backfit ship, geared motors 
ELEC 23 BASELINE 20: Backfit ship, direct-drive motors 
NEW MR ELEC 12: Rearranged ship, geared motors 
NEW MR ELEC 20: Rearranged ship, direct-drive motors 
CONSTANT FUEL ELEC 23 BASE 12: Backfit ship, geared 

motors 

CONSTANT FUEL ELEC 23 BASE 20: Backfit ship, direct- 

drive motors 



8.1. Direct effects 

The direct effects of an electric transmission are the 
changes in weight cuid volume of the propulsion system, as 
well as the transmission efficiency. Included are the 
weight and volume of the propulsion motors and generators, 
transmission lines, cooling systems, switchgear, power con- 
verter, exciter, braking resistor, any reduction gears and 
their associated gear oil, and shafting. These items are 
listed in Table 31. A positive difference from the baseline 
ship means a heavier and/or larger ship. 



21. Professor John Kassakian, MIT, private communication. 

22. Not all fuel in a ship is usable. There are nooks and 
corners of fuel tanks that are inaccessible to the fuel sys- 
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Note that the only variant that has a lower direct 
weight effect than the baseline is NEW MR ELEC 12. The ac- 
cumulation of weight increases in the others makes them 
heavier, while NEW MR ELEC 12 has lower motor and shafting 
weight than the rest. Geared drive is always lighter than 
direct-drive, largely due to the high weight of the direct- 
drive propulsion motors. With respect to volume, direct- 
drive is always smaller than geared drive, because of the 
reduction gears. All electric transmissions are larger than 
the mechanical baseline, but the smallest variants, within 
motor type groups, are the rearranged ships. 

Table 28 contains the maximum, sustained, and endurance 
speed transmission efficiencies of the two generator-motor 
combinations. Note that the of f -design-point efficiency of 
the direct-drive combination is significantly lower than the 
geared combination, even though it does not have the added 
inefficiency of reduction gears. This is in large part due 
to the poor efficiency of the slowly rotating direct-drive 
motor at the endurance speed. 

Table 28. Transmission efficiencies 

Combo Maximum Sustained Endurance 

12 0.9307 0.9266 0.8817 geared combination 

20 0.9209 0.9093 0.7754 direct-drive combination 

Endurance efficiency with one generator driving two motors. 

Table 29. Propulsion generator efficiencies 

PG ID Maximum Sustained Endurance 

SC 0.9891 0.9870 0.9737 geared combination 

SC 0.9891 0.9872 0.9768 direct-drive combination 

Table 30. Propulsion motor efficiencies 
PM ID Maximum Sustained Endurance 

S4 0.9898 0.9876 0.9526 geared combination 

SL8 0.9598 0.9496 0.8184 direct-drive combination 



115 



The above tables show the efficiencies of the motor and 
generator used in two particular combinations. The 
efficiency of the motor, reduction gears, and power con- 
verter have a direct effect on the efficiency of the gener- 
ator, as they change the loading point of the generator. 
Generally, motors and generators are more efficient when 
they are loaded more closely to their design point. (The 
same is true of gas turbines.) The inefficiencies of the 
power converters and reduction gears, if any, are included 
in the transmission efficiencies. 
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Table 31. Direct volume and weight effects 



CONSTftNT CONSTANT 















FUEL 


FUEL 




HECH 


ELEC 


ELEC 


NEH 


NEM 


ELEC 


ELEC 




23 


23 


23 


HR 


HR 


23 


23 


BASELINE BASELINE 


BASELINE 


ELEC 


ELEC 


BASE 


BASE 






12 


20 


12 


20 


12 


20 


Electric Propulsion weights: 


in LT 












PHs 


0 


18.08 


102.94 


18.08 


102.94 


18.08 


102.94 


P6s 


0 


10.27 


10.27 


10.27 


10.27 


10.27 


10.27 


Trans lines 


0 


0.18 


0.18 


0.77 


0.77 


0.18 


0.18 


Cooling sys 


0 


5.98 


5.98 


5.98 


5.98 


5.98 


5.98 


Switchgear 


0 


1.56 


1.56 


1.56 


1.56 


1.56 


1.56 


Power converter 


0 


7.16 


7.16 


7.16 


7.16 


7.16 


7.16 


Exciters 


0 


3.50 


3.50 


3.50 


3.50 


3.50 


3.50 


Braking resistor 


0 


10.00 


10.00 


10.00 


10.00 


10.00 


10.00 


Reduction gears 


78.90 


41.19 


0.00 


41.19 


0.00 


41.19 


0.00 


Shafting 


69.00 


66.84 


66.84 


41.72 


41.72 


66.84 


66.84 


M29B (op fluid) 


14.7 


14.7 


0 


14.7 


0 


14.7 


0 


Direct effect 


162.60 


179.46 


208.43 


154.93 


183.90 


179.46 


208.43 


Diff fe Baseline 


0 


16.86 


45.83 


-7.67 


21.30 


16.86 


45.83 


Electric Propulsion volutes: 


in cubic 


feet 










PNs 


0 


108.86 


909.13 


108.86 


909.13 


108.86 


909.13 


PGs 


0 


61.49 


61.49 


61.49 


61.49 


61.49 


61.49 


Trans lines 


0 


20.9 


20.9 


6.39 


6.39 


0.33 


0.33 


Cooling sys 


0 


200 


200 


200 


200 


200 


200 


SHitchgear 


0 


70.2 


70.2 


70.2 


70.2 


70.2 


70.2 


Power converter 


0 


1089.73 


1089.73 


1089.73 


1089.73 


1089.73 


1089.73 


Exciters 


0 


245.02 


245.02 


245.02 


245.02 


245.02 


245.02 


Braking resistor 


0 


1422.88 


1422.88 


1422.88 


1422.88 


1422.88 


1422.88 


Reduction gears 


2731.04 


1425.75 


0.00 


1425.75 


0.00 


1425.75 


0.00 


Shafting 


517.71 


501.5 


501.5 


306.06 


306.06 


501.5 


501.5 


Direct effect 


3248.75 


5146.33 


4520.85 


4936.38 


4310.90 


5125.76 


4500.28 


Diff fa Baseline 


0.00 


1897.58 


1272.10 


1687.63 


1062.15 


1877.01 


1251.53 



8.2, Indirect effects 

Indirect effects are again composed of weights and 
volumes, but these are the ripple effects of the propulsion 
system through the ship. For exaunple, if a transmission is 
more efficient at endurance speed, it should be expected 
that less onboard fuel would be needed to achieve the same 
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endurance range as a less efficient transmission. This is 
indeed the case. Another important indirect effect is the 
change in full load displacement, which is tied to the 
powering characteristic of the ship. 2 3 The following tables 
list the characteristics and indirect effects of the various 
ship configurations. 

Every electric transmission had lower maximum and sus- 
tained speeds than the baseline ship (by about 0.43 knots), 
but also lower EHP requirements at those speeds. The lower 
EHPs are a reflection of lower drafts (less resistance). 

The lower speeds show that the transmission efficiencies of 
the variant ships are lower than the mechanical baseline 
ship. There are more components in the electrical power 
trains, hence the lower efficiencies. The speed difference 
of 0.43 knots may be regarded by some as significant; it is 
about the speed reduction to be expected by a fouled bottom. 

The endurance range of all ships except those with con- 
stant fuel load is 5500 NM. The fuel load in the others 
varies greatly, showing the benefit of electrical cross- 
connection. In the two constant fuel ships, there was an 
increase in the endurance range of 1350 and 1400 NM, respec- 
tively, for the geared and direct-drive transmissions. This 
is an indication of fuel cost savings from the electric 
transmission. If a ship refuels every three steaming days 
(receiving a third of its tank capacity), steams 100 days 
each year, and fuel is priced at about $18 per barrel, this 
represents about a $600,000 savings per ship per year. 

The initial static transverse stability of the variants 
was degraded by the change in propulsion equipment. As pre- 
viously stated, the 6.5% 6M/Beam ratio of the mechanical 
baseline ship is not as large as desired for an actual ship, 
but provided a benchmark to measure relative changes. 

tern. Typically, 95% of the onboard fuel is usable. 

23. For the same molded lines, ships with higher displace- 
ments will have greater wetted surface areas and higher 
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Table 32. General ship characteristics 





HECH 

23 

BASELINE 


ELEC 

23 

BASELINE 

12 


ELEC 

23 

BASELINE 

20 


NEN 

HR 

ELEC 

12 


NEM 

NR 

ELEC 

20 


CONSTANT 

FUEL 

ELEC 

23 

BASE 

12 


CONSTANT 

FUEL 

ELEC 

23 

BASE 

20 


LBP 


425 


425 


425 


425 


425 


425 


425 


Ship disp FL 


5485 


5297.5 


5417.1 


5234.7 


5353.4 


5484. 1 


5530.2 


OiFf Baseline 


0 


-187.5 


-67.9 


-250.3 


-131.6 


-0.9 


45.2 


Nav beai 


55 


55.04 


55.02 


55.04 


55.03 


55 


55 


Nav draft 


16.44 


16.07 


16.31 


15.95 


16.18 


16.44 


16.53 


Depth sta 10 


38 


38 


38 


38 


38 


38 


38 


Accouodations 


272 


272 


272 


272 


272 


272 


272 


6«t/B 


0.065 


0.059 


0.058 


0.061 


0.064 


0.063 


0.061 


6H1 


950.01 


973.41 


959.1 


979.76 


967.39 


950.07 


943.41 


Full load K6 


23.23 


23.87 


23.75 


23.84 


23.53 


23.34 


23.44 


LCB/LBP 


0.516 


0.514 


0.515 


0.513 


0.514 


0.516 


0.516 


LCF/LBP 


0.565 


0.565 


0.565 


0.564 


0.565 


0.565 


0.565 


Metted surface 


23347.7 


23010 


23230.5 


22886.3 


23115.7 


23346.6 


23425.8 


Cp 


0.6 


0.596 


0.598 


0.594 


0.597 


0.6 


0.601 


Materplane area 


17641 


17575.9 


17624.8 


17539.8 


17602 


17640.9 


17653.2 


Hull volute 


618504 


618504 


618504 


618504 


618504 


618504 


618504 


Dkhs volute 


121528 


121528 


121528 


121528 


121528 


121528 


121528 


Total volute 


740032 


740032 


740032 


740032 


740032 


740032 


740032 


Endur fuel vol 


44365 


36592 


39647 


36367 


39396 


44364 


44364 


Kachy box vol 


109671 


109671 


109671 


109624 


109624 


109671 


109671 
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Table 33. Powering 



CONSTANT CONSTANT 





HECH 

23 

BASELINE 


ELEC 

23 

BASELINE 

12 


ELEC 

23 

BASELINE 

20 


NEW 

NR 

ELEC 

12 


NEK 

«R 

ELEC 

20 


FUEL 

ELEC 

23 

BASE 

12 


FUEL 

ELEC 

23 

BASE 

20 


Powering: 

Vaax 


29.13 


29.02 


28.83 


29.09 


28.9 


28.82 


28.7 


EHP 9 Vaax 


34821 


33131 


32741 


33154 


32763 


33066 


32700 


Vsus 


27.94 


27.82 


27.6 


27.88 


27.66 


27.63 


27.48 


EHP 9 Vsus 


28U>3 


26684 


26165 


26700 


26181 


26636 


26136 


Endurance 


5500 


5500 


5500 


5500 


5500 


6550 


6086 


EHP 9 Vend 


&851 


6632 


4772 


6557 


6698 


6851 


6905 


HPi 


51550 


51550 


51550 


51550 


51550 


51550 


51550 


kWi 


4000 


6000 


6000 


6000 


6000 


6000 


6000 


Avg 24 hr load 


2030 


2030 


2030 


2030 


2030 


2030 


2030 
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In the constant fuel ships, the high weight of the 
propulsion motors and propulsion motors and propulsion gen- 
erators, combined with a smaller amount of vertically- lower 
fuel reduced GM/B by over one-half percent, a not incon- 
siderable amount. For the rearranged variants, GM/B 
decreased less than the two previous ships. They also have 
less low fuel, but the propulsion generators are lower than 
in the backfit ships and the propulsion motors are very much 
lower. Only a very small decrease in GM/B was seen in the 
constant fuel ships because the constant fuel load compen- 
sated for the increased high weight of the electric trans- 
mission. The longitudinal metacentric height, GMl , in- 
creased for all variants, though it seems it should have 
decreased with the decrease in waterplane area and draft. 

The machinery space volume was the same for all ships. 

There were no big surprises in the area of weight. The 
structural weight (WlOO) encloses the same volume in every 
ship, so it was about constant. The propulsion plant weight 
(W200) varied with the type of transmission. Weight groups 
W300, W400 , W500, W600 , and W700 were virtually identical in 
every ship, and the variable loads were dominated by the 
change in fuel weight. The Design and Builders Margin is a 
function of the light ship weight (summation of WlOO through 
W700), so the margin weight moved with the light ship 
weight. The miles-per-gallon figure of NM/LT of fuel shows 
the endurance efficiency of electrical cross-connection. 

Only a few comments need be made regarding Table 35. 

The structural weight fraction shows the changes in full 
load displacement, remembering that the WlOO weights were 
all about the same. The same may be said for the weight 
fraction of the W300 through W700 groups and payload weight. 
Higher propulsion plant weights in the variant ships drove 
up the W200 fraction, except for NEW12. The fuel weight 
fraction shows the same behavior as the miles-per-gallon 
figure . 
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Table 34. Ship weights 



CONSTANT CONSTANT 





MECH 

23 

BASELINE 


ELEC 

23 

BASELINE 

12 


ELEC 

23 

BASELINE 

20 


NEH 

HR 

ELEC 

12 


NEH 

HR 

ELEC 

20 


FUEL 

ELEC 

23 

BASE 

12 


FUEL 

ELEC 

23 

BASE 

20 


Height: 
DFH onboard 


1049.1 


865.3 


937.5 


860 


931.6 


1049.1 


1049.1 


Usable Fuel Ht 


996.6 


822 


890.6 


817 


885 


996.6 


996.6 


DiTT T* Baseline 


0 


-183.8 


-111.6 


-189.1 


-117.5 


0 


0 


NH per LT fuel 


5.52 


6.69 


6.18 


6.73 


6.21 


6.57 


6.11 


Payload weight 


571.2 


571.2 


571.2 


571.2 


571.2 


571.2 


571.2 


HlOO 


1684.8 


1686 


1692 


1660.6 


1666.6 


1686.2 


1692.1 


H200 


343.6 


351.8 


387.1 


326.1 


361.3 


351.9 


387.1 


H300 


236.5 


236.5 


236.5 


236.5 


236.5 


236.5 


236.5 


H400 


302.2 


302.2 


302.2 


302.2 


302.2 


302.2 


302.2 


N500 


615.4 


613.5 


614.3 


613.5 


614.2 


615.4 


615.4 


N600 


426.6 


426.6 


426.6 


426.6 


426.6 


426.6 


426.6 


H700 


95.9 


95.9 


95.9 


95.9 


95.9 


95.9 


95.9 


Loads 


1317.1 


1120.9 


1193.2 


1115.6 


1187.2 


1305.1 


1304.9 


D&B largin 


463.1 


464.1 


469.3 


457.7 


462.9 


464.3 


469.5 


Disp FL 


5485.2 


5297.5 


5417.1 


5234.7 


5353.4 


5484. 1 


5530.2 






Table 35. Naval architectural analysis indices 



CONSTANT CONSTANT 



NECH 

23 

BASELINE 


ELEC 

23 

BASELINE 

12 


ELEC 

23 

BASELINE 

20 


NEH 

HR 

ELEC 

12 


NEW 

HR 

ELEC 

20 


FUEL 

ELEC 

23 

BASE 

12 


FUEL 

ELEC 

23 

BASE 

20 


L/D 


11.18 


11.18 


11.18 


11.18 


11.18 


11.18 


11.18 


L/B 


7.73 


7.72 


7.72 


7.72 


7.72 


7.73 


7.73 


B/T 


3.345 


3.425 


3.373 


3.451 


3.401 


3.345 


3.327 


GHl/LBP 


2.235 


2.290 


2.257 


2.305 


2.276 


2.235 


2.220 


HlOO/DTl 


0.307 


0.318 


0.312 


0.317 


0.311 


0.307 


0.306 


N200/DT1 


0.063 


0.066 


0.071 


0.062 


0.067 


0.064 


0.070 


M300/DT1 


0.0A3 


0.045 


0.044 


0.045 


0.044 


0.043 


0.043 


H400/DT1 


0.055 


0.057 


0.056 


0.058 


0.056 


0.055 


0.055 


M500/DT1 


0.112 


0.116 


0.113 


0.117 


0.115 


0.112 


0.111 


H600/DT1 


0.078 


0.081 


0.079 


0.081 


0.080 


0.078 


0.077 


N700/DT1 


0.017 


0.018 


0.018 


0.018 


0.018 


0.017 


0.017 


Hfuel/DTl 


0.182 


0.155 


0.164 


0.156 


0.165 


0.182 


0.180 


Hpayload/DTl 


0.104 


0.108 


0.105 


0.109 


0.107 


0.104 


0.103 


Nld/DTl 


0.240 


0.212 


0.220 


0.213 


0.222 


0.238 


0.236 


Vib/Vtot 


0.148 


0.148 


0.148 


0.148 


0.148 


0.148 


0.148 


HPi/Dfl 


9.398 


9.731 


9.516 


9.848 


9.629 


9.400 


9.322 


HPitViax/Dfl 


273.76 


282.39 


274.35 


286.47 


278.29 


270.91 


267.53 


kNi/DTl 


1.094 


1.133 


1.108 


1.146 


1.121 


1.094 


1.085 


Nt/DTl 


0.050 


0.051 


0.050 


0.052 


0.051 


0.050 


0.049 


Vtot/DTl 


134.914 


139.695 


136.610 


141.370 


138.236 


134.941 


133.816 


NlOO/Vtot 


5.100 


5.103 


5.122 


5.026 


5.045 


5.104 


5.122 


N200/HPi 


14.930 


15.287 


16.321 


14.170 


15.700 


15.291 


16.821 


Vib/HPi 


2.127 


2.127 


2.127 


2.127 


2.127 


2.127 


2.127 


H300/kNi 


0.039 


0.039 


0.039 


0.039 


0.039 


0.039 


0.039 


Vib/(HPi+kHi) 


1.906 


1.906 


1.906 


1.905 


1.905 


1.906 


1.906 


H500/Vtot 


1.863 


1.857 


1.859 


1.857 


1.859 


1.863 


1.863 


H600/Vtot 


1.291 


1.291 


1.291 


1.291 


1.291 


1.291 


1.291 


Dfl/Vtot 


16.603 


16.035 


16.397 


15.845 


16.204 


16.600 


16.739 
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For this technology characterization, everything 
devolves to total values. What is the total effect on the 
ship, once the individual pieces are put together? Table 36 
gives the bottom line. The electric propulsion plants are 
larger and heavier than their mechanical drive cousin; 
however, the extra weight and volume are more than compen- 
sated by the savings in fuel weight and volume. If a ship 
is designed from the beginning to be an "optimized" electric 
drive ship, over 6300 cubic feet of volume and 250 LT may be 
saved. The savings might be used for other systems, to 
reduce the overall size and cost of the ship (maybe allowing 
a larger buy, since 30 ships times 250 LT is a 7500 LT 
ship) , or to extend the naval architectural limits of the 
ship design. 

If a ship is backfitted with this technology, it is un- 
likely that tank volume can be recovered. However, the dual 
benefits of increased time-on-station and better fuel 
economy are realized. In this case, the choice between 
geared or direct-drive systems can be made by selecting the 
system with the most leverage, i.e., if the ship is volume- 
limited, use the lower volume direct-drive system (since the 
shafts are already in place). 

To put the volume and weight savings in perspective, 
note that 6300 cubic feet and 250 LT translates to twenty 
Tomahawk missile cells. The ship would be volume limited, 
with about 200 LT of weight savings still unused. This is a 
significant addition to the firepower of any ship, and the 
unused weight allows for ship growth. 
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Table 36. Total differences 





HECH 

23 

BASELINE 


ELEC 

23 

BASaiNE 

12 


ELEC 

23 

BASELINE 

20 


Total voluies: 
Fuel voluie 


44365 


36592 


39647 


Fuel difl 


0 


-7773 


-4718 


Prpln voluie 


3240.75 


5146.33 


4520.85 


Prpln diff 


0.00 


1897.58 


1272.10 


Total 


47613.75 


41738.33 


44167.85 


Total diff 


0.00 


-5875.42 


-3445.90 



Total Heights: 
Dll difference 


0.00 


-187.50 


-67.90 


prpln diff 


0.00 


16.86 


45.83 


fuel diff 


0.00 


-183.00 


-111.60 


other 


0.00 


-20.56 


-2.13 







CONSTANT 


CONSTANT 






FUEL 


FUEL 


NEM 


NEH 


ELEC 


ELEC 


HR 


HR 


23 


23 


ELEC 


ELEC 


BASE 


BASE 


12 


20 


12 


20 



36367 


39396 


44364 


44364 


-7998 


-4969 


-1 


-1 


4936.38 


4310.90 


5125.76 


4500.28 


1687.63 


1062. 15 


1877.01 


1251.53 


41303.38 


43706.90 


49489.76 


48864.28 


-6310.37 


-3906.85 


1876.01 


1250.53 



-250.30 


-131.60 


-0.90 


45.20 


-7.67 


21.30 


16.86 


45.83 


-189.10 


-117.50 


0.00 


0.00 


-53.53 


-35.40 


-17.76 


-0.63 
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8.3. Closure 



8.3.1. Conclusions 

This thesis has demonstrated the usefulness of electric 
drive transmissions in reducing ship weight and volume. 
Electric drive transmissions are better than mechanical 
drive transmissions on a ship basis. They provide, besides 
the weight and volume advantages, substantial arrangement 
flexibility and the opportunity to use new technologies in 
the ship design arena. The technical risk associated with 
these different technologies is minimal, as there is much 
industrial experience with electric machines, advanced 
switchgear, and the like. If the weight and volume reduc- 
tions are reinvested in the ship design through more optimum 
arrangements and subdivision, a substantially more efficient 
ship may be realized. Such a ship could successfully com- 
pete with the best of current ships . 

Small, light, high-power motors can be designed to a 
fair degree of detail with a computer optimization scheme if 
a meaningful objective function can be devised. For a ship 
system, the objective function should contain measures for 
volume, weight, efficiency, and relative cost (if a par- 
ticular material is significantly more expensive than other 
used) . A steepest-descent scheme can be combined with a 
Monte Carlo scheme to quickly converge on the objective 
function. 

The use of electric drive, and its consequent electri- 
cal cross-connect, can reduce the endurance fuel load by as 
much as 17.5%. When used in combination with an improved 
machinery arrangement and subdivision, that percentage can 
rise to 18%. If the fuel load stays constant, the endurance 
range may increase as much as 25%. 

On both an equipment weight basis and a ship weight 
basis, systems composed of a direct-drive propulsion gener- 
ator (with the same shaft rpm as the prime mover) and a 
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geared propulsion motor are better than those systems using 
no gears. Regarding volume, a non-geared system has lower 
equipment volume but higher ship volume due to the lower en- 
durance fuel efficiency. 

Geared propulsion motors have better of f -design-point 
efficiencies than those in direct-drive systems, primarily 
due to their higher rpm. A reduction in output power (in a 
motor) of 75% means only a few percent reduction in 
efficiency for a geared motor, while the same power reduc- 
tion means a 20% or more efficiency reduction for a direct- 
drive motor. 

Permanent magnet machines do not appear attractive for 
ship propulsion systems. They are both heavier and larger 
than candidate systems using synchronous, and, to a smaller 
extent, induction machines. Their low air gap flux density 
is the main detractor. Current permanent magnet materials 
cannot develop the energy product to compete with other al- 
ternatives, even though the NdFeB magnets are now in the 
marketplace. Induction machines may be useful as propulsion 
motors, but in this thesis they did not appear so. There- 
fore, ship propulsion generators and motors should only be 
synchronous machines. 

8.3.2. Recommendations 

The same modeling approach used in this thesis should 
be taken with variable reluctance machines (VRM) . Although 
no VRMs have been built at ship propulsion power levels, it 
is not inconceivable that they could serve in such a 
capacity. 

The induction machine model used here needs refinement, 
especially in the area of limiting maximum rotor current 
density. All of the machines need an analysis of their 
transient and dynamic characteristics. 

The recent advent of liquid hydrogen temperature super- 
conducting materials may signal an era where conventional 
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machines are overshadowed by the smaller, higher-flux 
machines possible with superconducting technology. However, 
if these new materials fail to provide the required current 
density, design of satisfactory machines may not be pos- 
sible . 

Integrated electric ship service propulsion plants may 
be beneficial additions to electric drive technology. Their 
influence on the systems suggested here may be an area of 
interest for future ship designs. 
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Appendix A. Definitions of machine variables and constants. 



bd 

BETA 

br 

BR 

BRl 

BREM 

BRNGS 

BSAT 

CP 

CRHO 

cw 

D 

dcore 

DCU 

DMAG 

doa 

dr 

ds 

eaf 

ef f cy 



ew 

f req 
g 



magnet operating point flux density, Tesla. Used 
in permanent magnet machines. 

hysteresis loss factor of M19 magnetic steel. The 
figure 2.5 was used. 

air gap flux density, Tesla. Used only in per- 
manent magnet machine model . 
air gap flux density, 1.05 rms Tesla, 
residual induction, 7.89 kilogauss, of M19 mag- 
netic steel, at BSAT and T1 . 
remanence flux of NdFeB, 1.21 Tesla, 
weight percentage of rotor shaft bearings, 1.03, 
or three percent of rotor weight, 
max flux density anywhere, 1.5 rms Tesla, 
stator coil pitch. The figure 0.8 was used, 
copper electrical resistivity, 1.724E-8 ohm- 
meters . 

copper weight, kg. 

density, 7.65e3 kg/m3, of M19 magnetic steel. 

back iron depth, meters. 

copper density, 8968.0 kg/meter^ . 

density of NdFeB, 7.4e3 kg/meter^ . 

overall machine diameter, meters. 

depth of rotor slots, meters. 

depth of stator slots, meters. 

p.u. internal voltage, used in syn only. 

efficiency, defined as 

output power 

effcy = 

output power + ph + pe + i2r + i2rr 

effective weight of machine, a combination of 
weight, volume, and efficiency. Used as the ob- 
jective criteria for the optimization scheme, 
machine synchronous frequency, Hz. 
air gap dimension, meters. 
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gams 


stator geometric factor, non-dimensional. Used to 
find convergence on active length. 


gamr 


rotor geometric factor, non-dimensional. Used to 
find convergence on active length. 


GMIN 


minimum machine air gap, 0.002 meters. 


HCl 


coercive force, 0.48 oersteds, of Ml 9 magnetic 
steel, at BSAT and Tl. 


il 


induction machine primary current, amps. 


i2 


induction machine load current, amps. 


i2r 


copper loss, watts, on stator. 


i2rr 


copper loss, watts, on rotor. 


jr 


full load rotor current density, amps/m* . 


jrnl 


no load rotor current density, used in syn only. 


js 


full load stator current density, amps/m* . 


JSMAX 


maximum stator current density, 12e6 amps/meter* 


ke 


efficiency weighting factor. 


KM 


magnet material cost weighting factor, defined as 


ks[] 


$ per pound of magnet material 
KM = - - - 

$ per pound of magnetic iron 

harmonic winding factors. 


kv 


volume weighting factor. 


1 


machine active length. 


le 


rotor winding space factor, used in induction 
machine . 


ler 


combined length of rotor end windings . 


Im 


radial dimension of permanent magnet. 


loa 


overall machine length, meters. 


Ir 


ratio of rotor slot width to slot pitch. In per- 
manent magnet machines, defined as ratio of magnet 
width to rotor "slot" pitch. 


Irat 


length ratio, used in permanent magnet program to 
calculate the effect of magnet overhang. 


Is 


ratio of stator slot width to slot pitch. 
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MAX_TIP_ 

minpwr 

MU 

MUR 

nr 

NU 

overhang 

P 

PC 

pe 

PF 

PI 
ph 
PS I 

r 

rlpl 

r2pl 

rev 

relpl 

RHO 

riv 

RMAX 

rpm 

RSF 

sev 

siv 

slip 



SPEED maximum allowable rotor tangential velocity, 
200 meters per second. 

minimum required power of the machine, horsepower. 
MU of air, 4TTe-7 Henries/meter. 

NdFeB relative reversible permeability. The 
figure 1.05 was used. 

number of rotor bars in induction machine, 
anomalous loss factor of M19 magnetic steel. The 
figure of 2 was used, 
percentage of magnet overhang, 
number of machine pole-pairs. 

permanence coefficient of NdFeB. The figure 1.1 
was used. 

eddy current loss, watts. 

power factor of all machines. The figure 0.8 was 
used. 

3.141592654. 
hysteresis loss, watts. 

Trated/Tpullout for induction motor. The figure 

0.55 was used. 

rotor radius , meters . 

per length primary resistance, ohms /meter. 

per length secondary resistance, ohms /meter. 

rotor copper volume, meters^ . 

per length Thevenin equivalent resistance. 

electrical resistivity, 52 micro-ohm-cm, of M19 

magnetic steel . 

rotor iron volume, meters® . 

maximum rotor radius, 2.0 meters. 

machine shaft revolutions per minute, referred to 
as rotor speed. 

rotor slot space factor. The figure 0.35 was 
used . 

stator copper volume, meters®. 

stator iron volume, meters®. 

guessed slip of the induction machine. 
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slipl 


larger derived slip. 


slip2 


smaller derived slip. 


smax 


maximum machine slip. 


SSF 


stator slot space factor. The figure 0.35 was 
used . 


T1 


thickness, 0.014 inch, of M19 magnetic steel. 


tmaxpl 


per length maximum torque . 


vlapl 


per length Thevenin equivalent voltage. 


va 


IP+jQI , VA rating, volt-amps. 


vagpl 


per length air gap voltage. 


vol 


machine envelope volume, meters^ , with a margin. 


VOLALL 


volume allowance for frame and foundation. A ten 
percent allowance was used. 


volmag 


volume of permanent magnet material, meters^ . 


\rtpl 


per length terminal voltage. 


w 


electrical frequency, radians per second. 


wr 


width of rotor slots. 


wm 


mechanical angular velocity of rotor, radians per 
second. 


wt 


weight of copper and iron in a machine, plus a 
margin. 


WTALL 


weight margin for fraone and foundation. A ten 
percent margin was used. 


vrt_iron 


iron weight, kg. 


wtmag 


weight of permanent magnet material, kg. 


xlpl 


per length primary impedance. 


x2pl 


per length secondary impedance. 


xbeltpl 


per length belt impedance. 


xelpl 


per length Thevevin equivalent impedance. 


xmpl 


per length air gap magnetizing impedance. 


xrdpl 


per length rotor differential leakage impedance. 


xrspl 


per length rotor slot leakage impedance. 


xs 


p.u. synchronous impedance. 


xsepl 


per length stator end turn impedance. 


xslot 


slot impedance. 


xsdpl 


per length stator differential leakage impedance. 
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xsspl 

xzzpl 



per length stator slot leakage impedance, 
per length zig-zag impedance. 
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Table 37. Listing of various functions used in computer 
programs 



Ittt a pseudo randoa nuaber generator 1)1/ 


♦define 


HULTIPLIER 


25173 


♦define 


nODULUS 


3276B 


♦define 


INCREHENT 


13849 


♦define 


HODFLT 


32768.0 


double randoiO 





{ 



extern long int seed; 

seed=(HULTIPLIERt5eed+INCREHENT) Z MODULUS; 
return(seed/HODFLT); 

} 



double 5wf(n) /t stator winding factor 1/ 

int n; It hartonic order t/ 

{ 

double kp, kb; 

extern double cosO, sinO; 

kp=cos(0.3142tn); ft pitch factor, assuaes O.B coil pitch 1/ 

kb=lsin(0,5236tn))/(0.5236tnl; It breadth factor, froi 

Kirtley’s "Basic Foraulas ... " and assutes an 
electrical winding angle of 60* t! 

return(kplkb); 

} 



double abs(q) 

double q; 

{ 

if (q < 0.0) 

q= (-l.O)tq; 
return (q); 

} 



double sinh(u) 

double u; 

{ 

double expO, ans; 

ans = 0.5t(exp(u) - exp(-u)); 
return (ans); 

} 



double cosh(u) 

double u; 
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{ 

double expO, ans; 

ans = 0.5t(exp(u) + exp(-u)); 
return Ians); 

} 

double rwf(n) /I rotor winding factor, saie coiients as swfO M 

int n; 

{ 

double kp, kb; 
extern double sinO; 
kp=l; 

kb=(5in(0.5236lnl)/(0.5236ln); 

return(kptkb); 



tdefine tol 1.0e*6 
IdeHne pi 3.141592654 

double best (p,x) /t bessel 1 (order, arg) (/ 

double x; 
int p; 

{ 

double bi; 
int biaO; 

double exp(),5qrtn,bi5(); 
it (p<0) abort (" besi: negative index"); 
it (x<0) abort (" besi; negative arguient"); 
it (x>60) bi=exp(x)/5qrt(2lpi$x); 
it ((x<12) ! ibia(p,x,&bi)) bi-bis(p,x); 
return (bi); 



double bis(p,x) 
double x; 
int p; 

{ 

double tabsO; 
double xx; 
int i,tk,k; 
double bi=0; 
double t=l; 
xx=x/2.; 

for (i=l;i<p+l;i++) t=ttxx/i; 

if (t<l.0e-36) bi=0; 

else 

{ 

bi=t; 

xx=xx(xx; 

for (k=li(k<1001)W((fab5(t)-fab5(biltol))>0);k++) 
{ 

fk=k((p+k); 
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t=tlxx/fk; 

bi=bi+t; 

} 

\ 

i 

return (bi); 

} 

int bia(p,x,pbi) 
double Xytpbi; 
int p; 

{ 

int ist,fn,lt,fk; 

double xx,bi,t; 

double ^abs(),sqrt() ,exp() ; 

fn=4lplp; 

t=1.0; 

bi=l.O; 

xx=0.125/x; 

kor (k=l;(k<301bii((fabs(t)-fabs(biUol))>0)|k++) 

{ 

4k=(2lk-m(2lk-l); 

t=ttxxt(fk-fn)/k| 

bi=bi+t; 

} 

If (k==311 ist=l; 
else 

{ 

ist=0; 

bi=bilexp(x)/sqrtl2.0tpitx); 

} 

)pbi=bi; 
return list); 

> 

double besk(p,x) /) aodified bessel function Kp(x) 1/ 

double x; 
int Pi 

{ 

double bk; 

double exp(),sqrt(),kO(),ki(); 

if (p<0) abort (' besk; negative index *); 

if (x<0) abort (' besk: negative arguient '); 

if (x>60) bk=exp(-x)/sqrt(2.0tx/pi); 

else if (p==0) bk=kO(x)| 

else if (p==l) bk=kl(x); 

else { 

double gO,gl,gj; 
int j; 
gO=kO(x); 
gl=kl(x); 

for (j*2;j<p+l;j++) 

{ 

gj=2l(j-l)lgl/x+g0; 

gO=gl; 
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Ql'Qii 

) 

bk=gj; 

J 

return (bk); 

) 

double kO(ii) 
double k; 

{ 

double bk; 

double log(),sqrtl),exp(l; 
if (x<l) 

{ 

double a,b, 2 ,c,g 0 ,x 2 j,f,hj,rj; 

int j; 

b=0.5tx; 

a=.5772i566+log(b); 

c=btb; 

gO=-a; 

x2j=i; 

f=l; 

hj=0; 

for (j=i;j<7;j++) 

{ 

rj=l.0/j; 

x2j=x2jtc; 

f=flrjtrj; 

hj=hj+rj; 

g0--g0+x2jlft(hj-a); 

) 

bk=gO; 

} 

else 

{ 

double tU2],a,b,c,pa; 
int 1; 
a=exp(-xl| 
b=1.0/x; 

C'sqrtibi; 

tC0]=b; 

for (l=l;K12;l++) tm=tll-l]tb; 

pa = I.2533l4l-.1566642m01; 

pa += .088lll2Blt[l]-.09l39095ltI2); 

pa += .I344596ttl3]-.2299850tt[4]; 

pa += . 3792411415]-. 5247277ltt6]; 

pa += .55753681417]-. 426263311181; 

pa += .2l84518ltl9]-.06680977ltll0]+.009l893831tlll] 

bk=atctpa; 

} 

return Ibk); 
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double kl(x) 
double x; 

{ 

double bk; 

double log(),exp(),sqrt(); 
if (x<l) 

( 

double a,b,c, 9 l,x 2 j,f,hj,rj; 

int j; 

b=x/2; 

a=.57721566+log(b); 

c=btb; 

x2j=b; 

f=l; 

hj=I; 

gl=l.0/x+x2jl(.5+a*hj) j 
for (j=2;j<9; j++) 

{ 

x2j=x2jtc; 

rj=1.0/j; 

f=flrjtrjj 

hj=hj+rj; 

9l=gl+x2jlfl(.5+la-hj)lj); 

) 

bk=gl; 

} 

else 

{ 

double a,b,c,t[121,pa; 

int 1; 

a=exp(-x); 

b=I.0/x; 

c=sqrt(b)i 

ttO]=b; 

for (1=1;1<12;1++) t(l]=ttl-l]lbj 

pa = l.2533l4l+.449997ltt0]-.1468583Um; 
pa += .1280427lt[2]-.1736432lt[3]+.28476l8tt[4]; 
pa += -.4594342ttC5]+.6283381$tt6]-.6632295ltt7]; 
pa += .5050239lt[8]-.258l304lt[9]; 
pa += .0788000lttn0] -.010824I8U[ll]; 

bk=atctpa; 

} 

return (bk); 

} 

double besip (p,arg) /t deriv of bes I (order, arg) M 

int p; 
double arg; 

{ 

double x,y,z; 

X = besi (p-l,drg); 
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y = p t besi (p,arg)/arq; 
z = X - y; 
return (z); 

} 

double beskp(p,arQ) /t deriv of bes Kforder, arg) 1/ 

int PI 
double arg; 

{ 

double x,y,z; 

X = "besk (p-l,arg); 
y = - p t besk (p,arg)/arg; 
z = X + y; 
return (z); 

} 
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Appendix B. Synchronous Machines and General Relations 
An equivalent circuit for a synchronous machine is 



Ra Xs 




Figure 35. Synchronous machine equivalent circuit 

where Ra is the stator resistance and Xs is the synchronous 
reactance. The internal voltage of the machine, Eaf, is 
developed between the stator and rotor across the air gap. 
The current direction is shown as if the machine is a motor. 

Vta = Rala + jXsIa + Eaf 

Eaf represents a mutual coupling between the stator and 
rotor, and 



w M If 4 Po kl kf 1 r Ns Nf 

Eaf = where M = 

^2 w g p* 



If is field current, Nf is the number of series field turns 
and kf is the fundamental rotor winding factor. Since it is 
never desired to specify the number of turns on either the 
stator or the rotor, a scheme has been devised so that 
derivations are conducted in volts -per-turn, ampere-turns, 
and ohms-per-turn-squared, which results in power in watts. 
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Figure 36. Phase belt conductor area 



For the stator, the transverse area can be divided into 
phase belts. The area labelled 'A’ is the area occupied by 
phase 'a’, in one direction. Then, Hs la = Aa Ja , where Ja 
is the current density in phase 'a’ conductors. The area 
has some effective conductor area, subject to the need for 
conductor cooling passages and insulation area. Therefore, 
for stator and rotor currents, analysis yields 



Aa Ja 




Af Jf 

and If = 

Nf 



It is also desired never to specify the number of slots on 
either the stator and rotor. Accordingly, slot space fac- 
tors are defined as 



h 



slot 
Wid-t h 









pl't'ch 



Figure 37 . Slot space factors 
stator slot width 

Is 

stator slot pitch 
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Ir 



rotor slot width 
rotor slot pitch 



For a typical turbogenerator conductor bar, the copper area 
is about thirty-five percent of the conductor envelope area. 
Variables titled SSF (stator slot space factor) and RSF 
(rotor slot space factor) embody this thirty-five percent. 
The conductor area of a single stator phase in one direction 

(r+g) 2tt Is ds SSF 

Aa = and, for the rotor, 

6 

r 2tt Ir dr RSF 

Af = There is only one phase on the 

2 rotor. 

The armature resistance is 

rho It Ns 

Ra = 

Aa/Ns 

where rho is the electrical resistivity of copper and It is 
the turn length. If a stator turn can be modeled as 




Figure 38. Stator turn 



and the circumfrence is 2rr(r+g), then 



It 



21 + 2 { 



p 



} 



where -1(4/3) is (l/sin(60®)) and CP is the coil pitch. Then 
Ns2 TT (r+g) Is ds SSF 
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The stator and rotor copper losses are 



Rala^ = rho It Aa Ja^ Rflf^ = rho lt£ Af Jf^ 

The rotor is full pitched. 

The vector diagram for the synchronous machine equiv- 
alent circuit is (assuming Ra ==> 0) 




Figure 39. Synchronous machine vector diagram 

By use of the law of cosines, 

Eaf2 = Vt2 + (Xsla)2 - 2VtXsIa cos(r + 90«) 

where t is the power factor angle. If this equation is put 
into per unit form, with eaf = Eaf/Vt and xs = Xsla/Vt, then 

eaf2 = 1 + xs^ + 2 xs sin(r) 

xs can be calculated, allowing the calculation of eaf. 

There is simple linear relationship between eaf, the no-load 
rotor current density (jrnl), and the full-load rotor cur- 
rent density (jr). It is jr = (eaf) (jrnl), because Eaf is 
directly proportional to If and If is directly proportional 
to Jf . Ampere’s Law states 




If one chooses an integration path around half of the rotor, 

4 Nf If kf 

Po H 2g = Po 

TT P 
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If the no-load condition is chosen, and it is recognized 
that Br is essentially constant with increases in Vt 
(because when Vt increases, so must Eaf ) , then 



Br 



2 Mo r Ir dr RSF kf jrnl 
g P 



and jrnl 



Br_g p 

2 Mo r Ir dr RSF kf 



The synchronous reactance is a measure of flux linkage 
within the machine. It is Xs = w Ls = w (Lai + 1.5 Laa), 
where Laa is the armature single phase inductance and Lai is 
the slot leakage inductance. The factor of 1.5 is derived 
from the 120® separation between the three phases of these 
machines, as 



Xsla = w (la Lai +IaLaa + IbLab + IcLac) 



Due to the symmetry of the machine. 



Lab = Lba = Lac = Lea = Lbc = Lcb = -0.5 Laa 



Then, Xsla = w la (Lai + 1.5 Laa). 
Self-inductance is 



4 Mo kfl^ Ns^ 1 r 

Laa = 

TT g p2 

from Ampere’s Law and L = ( f lux/curront ) . If a single con- 
ductor per slot is postulated and the effects of slot teeth 
are ignored, then the leakage inductance is 



slots 

Lai = (Pself + Pmutual ) Ns* 

(pole) (phase ) 
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I 

I 



I 



where P is reluctance. For this conductor configuration, 



Mo 1 

Pself = - - ds ws 

3 


Mo 1 

and Pmutual = -- - ds ws, 

2 


with ws equal to slot width. 


Since ws is not known, we use 


the stator slot space factor, 
of slots to yield 


Is, multiplied by the number 


(slots) ws = Is TT (r+g) 


Therefore, 



5 Ns^ Po 1 ds Is TT (r+g) 
Lai = -- 



36 p 

The real power developed 


by the machine is 


Pwr = 3 Vt la pf , where pf is 


the power factor. Through the 


use of Lenz ’ Law and Ampere ’ s 
expressed as 


Law, terminal voltage may be 


2 r 1 w Br ks 

Vt = 

P 





Using our previous relation for la , the expression for Pwr 



is 




2tt r 1 w Br ks Ja 

Pwr = 

P 


SSF (r+g) Is ds pf 


Finally, winding factors 
ing breadth factor, kb, is 


need to be derived. The wind- 


sin(m n r/2) 

kbn= 

m sin(n r/2) 





where r is the electrical angle between adjacent slots, n is 
the harmonic order, and m is the number of slots per pole 
per phase. The winding pitch factor, kp , is 
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kpn = sin (n cx/2) where cx is the electrical angle 



between sides of the coils (pitch angle). For a three phase 
winding, a = 2tt p. If a 0.8 coil pitch is assumed (which 
will rid the machine of certain harmonics during balanced 
operation) then « = 180® (1 - 0.8) = 36® and 
kpn = sin(0.3142 n). Assumptions are needed to calculate 
kbn without specifying the number of stator slots or turns. 
If the winding angle is specified as did Kirtley [32], then 



kb n 



sin(n Ow/2) 
(n ©w/2) 



A reasonable electrical winding angle is 60® , since roost of 
the stator periphery will contain turns. The breadth factor 
devolves to kbn = sin(0.5236 n)/(0.5236 n) , for which the 
fundamental harmonic factor equals 0.955. The winding fac- 
tor is the product of the breadth and pitch factors. For 
the rotor, a pitch factor of one is assumed. 
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Table 38. Listing of synchronous machine design program 



iinclude "stdio.h* 
linclude “def.h* 



It prograt naie: syn.c for synchronous, round rotor lachines 1/ 

long int seed; /t start point for randoi nutber generator tl 

double b[26H26], h[ll][26], ks[8], kr[8]; 

/t b is "best" array, h is “hold" array, ks/kr are winding factors 1/ 

•ainO 

{ 

double designjointO , rnd^walkO, swfO, rwfO, ke, kv, linpwr, 
stepsize, randoiO, absO, freq, rpa; 
int p, iteration, i, j, best, print_out(), loops; 

FILE tfopenO, tfp; 

printf (•\nReading input data froi SYN.DAT . . 



fp=fopen(“syn.dat'',V'l; /t input seed for randoi nuibers 1/ 
fscanf (fp,"Xd“,liseedl; 
fscanflfp, "Xd*, ip); 
fscanf(fp, "Xlf“, iiinpwr); 
iinpwrt=746.0; 
fscanf(fp, 'Xlf, ike); 
fscanftfp, 'Xir, ikv); 
fscanf(fp, nif, irpt); 
fclose(fp); 



/t input nuiber of pole pairs 1/ 

/I input lachine power, derived ft ASSET )/ 
/t convert to watts tl 
!% CERs for Effective Weight )/ 

i\ tachine tax shaft rpt %i 



printf ("\nHow tany loops do you want? “); 
scanfCXd*, iloops); 

printf (*\n\nDoing prograt calculations . . .\n"); 



for (i=l; i < 8; i+=2) It hartonic winding factors M 

( 

ks[il=swf (i); 
kr[i]=rwf (i); 

} 



freq=rptlp/60.0; It tax electrical frequency t! 

It MAIN BODY OF THE PROGRAM t! 



for (i=l; i <= loops; ++i) 

( 

stepsize^O.l; 

iteration=0; 

design^point (iinpwr, p, ke, kv, freq); 

It put stuff in the hold array t! 

while (iteration <= 10) 

{ 
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rnd_Mdlk(iinpHr, p, stepsize, ke, kv, freq); 

/t stagger around t/ 

best=0; /I index to best EM of the lot 1/ 

tor (j=l; j<=10; ++j) 

if (hUHlSl < hlbestHlBl) 

best=j; /t find the best lachine t/ 

if (abslIhlOinSl - hlbestHlBD/hlOHlB]) < 0.005) 

/t siall iiproveient in EH t/ 

{ 

stepsize/=2.0; 

^iteration; 

) 

else /I transfers best to 0 position t/ 

{ 

for (j=l; j <= 25; ++j) 

hfOlUl = hlbestlljl; 

} 

) 

for (j=l| j <= 25; ++j) 

b[i)[j]=h[bestl[jl; /I keep the best lachine 1/ 



best=l; 

for (i=l; i <= loops; ++i) 

if IbCilElB] < blbestniBl) 

best^i; H find and keep the best of the best t/ 

•inp»tr/=746.0; H turn back into hp 1/ 

print_out(best, p, linpHr, ke, kv, rpi); /t output to disk file t/ 



fp=fopenl'syn.dat', ■»»■); /I output seed 1/ 

fprintf(fp,’Xd*, seed); 

tprintflfp,’\n7d", p); 

fprintf (tp,’\n21f iinpiir); 

tprintflfp,'\nXlf", ke); 

fprintf (fp, “\n21f, kv); 

fprintf Ifp, ■VnXlt", rpi); 

fcloselfp); 

} 

/t END OF HAIN PROBRAN; ALL THAT FOLLDH ARE FUNCTIONS t/ 



double design_point(iinpHr, p, ke, kv, freq) 

/t deteriines a randoi design point t/ 

double ainpHr, ke, kv, freq; 
int p; 

{ 

double r, jrnl, jr, js. Is, Ir, dcore, ds, dr, g, h, 1, xs, eaf, i2rr, 
va, ph, pe, i2r, vol, Ht, effcy, eH, xsl, xs5, xs7, xsal, 
c», siv, scv, riv, rev, loa, doa, t»t_iron, tind_sivl); 
extern double sqrtO, randoal); 
int c=0, d=0; 

while (d != 1) { 
while Ic != 1) 
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{ 

r=rindoiOtRnA)(; 

if ((2IPUrlfreq/p) < HAX_TIP_SPEED) /I check rotor tip speed 1/ 
break; 

} 



w=2IPllfreq; 




lr=randoi()t0.5 * 0.25; 


It rotor slot factor tl 


dr=randoi()lr/5.0; 


It slot no deeper than 20X of rotor radius tl 


dcore=(BRIr)/(BSATtp); 


It back iron depth tl 


ds=randoi( ) tO. 9tdcore; 


It slot depth ( 901 of body depth tl 


while (c != 1) 


It gap diiension tl 



{ 



Q^randoiOtfO.ltr - 6H1N) * 6H1N; 
if (g > 0) 

break; 

} 

ls=randotnt0.5 + 0.25; /I stator slot factor 1/ 

js=randoi()tJSHAX; Jt full load stator current density t/ 

jrnl = (BRJgtp)/(8lt1UIRSFtrtdrHrlkrn]); It noload rotor current density tl 

xsl=ks[l]tksni; 

xs5=(ks[51lks[51/25); 

xs7=(ks[7nks[71/49); 

xsal=(5tlstPltds/18); 

xs=(«UljslSSFI(r+9)l(rtg)JdslPllls I (12»(xsl + xs5 ♦ xs7)/(PUgJp) + xsaD) 

/(l2trtBR<ks[ll); 

It p.u. synch iipedance tl 

if (xs > 2.0) 

continue; It don’t uant xs too big tl 

else 

++d; 

eaf=sqrt(l * xstxs * 2txst0.6); It 0.6 is sin(r), pur factor angle, 
eaf is p.u. internal voltage at full load tl 

jr=eaftjrnl; It jr full load, linear uith eaf tl 

l=(iinpwrtp)/(2tBRlHlks[l]trtSSFtdstPII(r+g)UstjstPF); It active length tl 

va=2tPUrlllwlBRtksCl]tjstSSFt(r+g)ldstls/p; It va rating tl 

siv = find_siv(l, r, g, ds, dcore, Is); It stator iron voluie tl 

riv = lIPltrKr - 2ldrllr); It rotor iron voluie tl 

scv = 2tPIt(r+gUdstlsl(l + 2.3094IPll(r+g)tCP/p); 

It stator copper voluie tl 
rev = 2tPUrldrtlrt(l + 2.3094tPllr/p); It rotor copper voluie tl 

CM = (rev + scvUDCU; It total copper weight tl 

loa = 1 + 4t(r+g); It length-over-all tl 

doa = 2t(r+g+ds+dcore); It over-all -diaieter tl 

vol = VQLALLt(loatPltdoatdoa/4) ; It lachine envelope voluie tl 

Mt = HTALLtlcN * DKBRNGStriv * siv)); It lachine weight in kg tl 
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wt^iron = DKriv + siv); 



!% iron weight only t/ 



ph = 31.86225IBETAHreqlBRllHCllwtJron/D; 

!% hysteresis loss in watts, uses iron weight of wachine 1/ 

pe = (106236.9INUIBSATIBSATITllTllfreqlfreqtwtJron)/(RH0ID); 

/I eddy current loss in watts, uses iron weight of ■achine H 

i2r = 2.0ISSFIdslCRH0IPII (1 + 2.3094IPII(r+g)tCP/p) IjsIjsKr+gllls; 

!% stator copper loss in watts 1/ 
/t revised 1-12-87 M 

i2rr = 2.0IRSFtdrtrllrtjrljrlCRHOtPU (1 + 2.3094>PItr/p); 

!% rotor excitation losses, 1-12-87 %! 
i2r+=i2rr; /I total copper losses >/ 



effcy=(iinpwr)/(iiinpwr + ph + pe + i2r); 



ew=wt + kel(l-effcy) + kvtvol; /t Effective weight 1/ 



h[0][l]=js; 


hlOH2]=freqi 


hI0][3]=»*i 


h[0H4]=r; 


ht0][5]=9} 


ht0][6]=dcore; 


h[0][7]=ds) 


h[0][8]=dr; 


h[0][9]=ls; 


h[0][10]=lr; 


h[0][in=vol; 


htO][12]=Htj 


h[0][13]=ph; 


h[0][14]=pe; 


ht0][15]=i2r 


ht0][16]=va; 


h[0][17]=effcyi h[0][18]=eHi 


h[0][19]=lj 


h[0][20]=jr; 


ht0][21]=jrnl; 


hC0][22]=xsi 


ht0][23]=eaf; 


h[0][24]=loa; 


h[0H25]=doa 



/t this section just changed all the variables in the "hold" array %! 



return; 

}} 



double rnd^walkiiinpwr, p, stepsize, ke, kv, freql 

/I walks about design^point 10 tiaes t/ 
double stepsize, ainpwr, ke, kv, freq; 
int p; 

{ 

double r, jrnl, jr, js. Is, Ir, dcore, ds, dr, g, w, 1, xs, eaf, i2rr, 
va, ph, pe, i2r, vol, wt, effcy, ew, xsl, xs5, xs7, xsal, 
cw, siv, scv, riv, rev, loa, doa, wt^iron, find^sivO; 
extern double sqrtO, randoaO; 
int i=l; 



while (i <= 10) 

{ 

/I read in the walk around the design point 1/ 

js=h[0][l]Ml + stepsizeKrandotO - 0.5)); 
if (js > JSMAX) 

js = JStlAX; n reset to liait M 

w=2IPItfreq; 

r=h[0][4]l(l + stepsizel (randoaO - 0.5)); 
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if ((Nir/p) > HAXJIP.SPEED) 

continue; /t go to next try if violated 1/ 

g=h[0][S]t(l * stepsnetirandoAf ) - 0.5)1; 
if (g < 6H1N) 

g=6HlN; It reset to tlie litit t/ 

dcore=(BRtr)/(6SATtp); It lost efficient use of iron t/ 

ds=h[0][7]t(l * stepsizeMrandoiO - 0.51); 

if Ids > dcore) It can’t tiave too-deep slots t/ 

ds=dcore; It reset to the liait 1/ 

dr=h[0][81t(l ^ stepsizetlrandoiO - 0.5)); 
ls=hC0][9]t(l * stepsizetirandoaO - 0.5)); 
if (Is > 0.75) 

ls=0.7S; It reset to the litit t/ 

if (Is < 0.25) 

ls=0.25; 

lr=h[01[101t(l * stepsizeUrandoaO - 0.5)); 
if (Ir > 0.75) 

lr=0.75; It reset to the litit 1/ 

if (Ir < 0.25) 

lr=0.25; 

It cotputation section of the Malic t/ 

jrnl=(BRIgtp)/(8tMUtRSFIrtdrtlrtkr[l]l; It noload rotor current density 1/ 

xsl^lcsnitksill; 
xs5=(ks[5]lks[5]/25); 
xs7=(ks[7Hks[71/49); 
xsal=(5tlstPItds/lB) ; 

xs=(l1UtjstSSFt(r+g)l(r+g)»dstPItls t (12l(xsl + xs5 + xs7)/(PIIg»p) + xsaD) 

/(12trtBRtks[ll); 

It p.u. synch itpedance tl 

if (xs > 2.0) 

continue; /( can’t have xs too big tl 

eaf-sqrtd * xstxs * 2txst0.6); It O.b is sin(T), ptr factor angle, 

eaf is p.u. internal voltage at full load tl 

jr=eaftjrnl; It jr full load, linear Mith eaf tl 

l:(iinpHrtp)/(2tBRtHtksCl]trtSSFtdstPII(r^g)tlstjstPF); It active length tJ 

va=2IPltrm«lBRtks[l]tjstSSFI(r+g)tdstls/p; It va rating tl 



siv = find.sivd, r, g, ds, dcore. Is); 
riv = llPlirllr - 2tdrtlr); 
scv = 2IPIt(r+g)tdstlsl(l + 2.3094tPH(r+gdCP/p); 



rev = 2tPIIrtdrllrt(l + 2.3094IPI»r/p); 

CM = (rev + scvltDCU; 

loa = 1 + 4t(r+g); 

doa = 2l(r+g+ds+dcore); 

vol = VOLALUdoatPItdoatdoa/4) ; 

Mt = UTALLKcm * DdBRNGStriv * siv)); 



/( stator iron volute tl 

It rotor iron volute tl 

It stator copper volute tl 
It rotor copper volute tl 
It total copper Height tl 
It length-over-all tl 
It over-all-diateter tl 
It tachine envelope volute tl 

It tachine Height in kg tl 
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wt_iron = DKriv + siv); 



IX iron weight only M 



ph = 31.86225lBETAtfreqtBRltHCKwtjron/D; 

/t hysteresis loss in watts, uses iron weight of iachine 1/ 

pe = (106236.9tNUtBSATtBSAmmUfreqtfreqtwt_iron)/(RH0tD); 

it eddy current loss in watts, uses iron weight of tachine t/ 

i2r = 2.0tSSFIdstCRH0JPIt (1 + 2.3094tPU(r+g)tCF7p) Ijsijst ir+gHls; 

/t stator copper loss in watts 1/ 
n revised 1-12-87 M 

i2rr = 2.0tRSFtdrlrtlrljrljrtCRH0tPI< (1 + 2.3094tPUr/p); 

ft rotor excitation losses, 1-12-87 1/ 
i2r+=i2rr; IX total copper losses X! 

effcy=(iinpwr)/(iinpwr + ph + pe + i2r); 

ew=wt + ket(l-effcy) + kvivol; IX Effective weight X! 



h[i][l]=js; 


h[i][2]Mreq; 


h[i}[3]=H; 


h[i][4]=r; 


h[i}[5]<9; 


h[i][6]=dcore; 


h[i][7]=d5| 


hm[8]=dr; 


h(i][9]=ls; 


h[i][10]=lr; 


h[i][in=vol; 


h[i][12]=Ht; 


h[i][13]=ph; 


h[iH14]:pe; 


hti][15]=i2r 


h[i][16]=vaj 


h[i][17]=effcyj 


ha][18]=eH; 


h[i][19]=l; 


h[i][20]=jr; 


h[i][21]=jrnl; 


h[i]I22]=xs; 


h[i][23]=eaf; 


h[iH24]=loaj 


hti][25]*doa 



IX this section just changed all the variables in the "hold" array XI 



++i; IX go to the next hlilU XI 

) 



return; 

} 



print_out (best, p, tinpwr, ke, kv, rpi) 
int best, p; 

double Minpwr, ke, kv, rpt; 

{ 

char outfile[14]; 

FILE tfpo, tfopenO; 
int i; 

printf ("\nHhat is the nate of the file where you want the output? "); 
scanf("ls", outfile); 

fpo=fopen(outfile, “w"); 



fprintf(fpo,"Zd", p); 
fprintf (fpo,“\nXlf ", linpwr); 
fprintf(fpo,"\nllf, ke); 
fprintf(fpo,"\nXlf", kv); 
fprintf(fpo,"\nllf", rpt); 
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for li=l; i <= 25; ++i) 

fpnntf lfpo,’\nXH‘,b[best][i]); 
fprintf(fpO) *\n’); 
fclose(fpo); 

} 



double find_siv(l, r, 9 , ds, dcore, Is) 
double 1 , r, g, ds, dcore, Is; 

{ 

double one, tHO, three, four; 

one = lr+q+ds+dcore)tlr+q+ds+dcore) - lr+q)llr+q); 
tNo = 2tPlt(r+q)tds»ls; 

three = (r+q+ds+dcore)l(r+g+ds+dcore) - (r+q+dsi t(r +q+ds) 

four = IKPItone - two) + Plt4l(r+q)lthree; 
returnlfour); 

} 
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Table 39. Listing of synchronous efficiency program 



tinclude “stdio.h" 
linclude "def.h“ 



/I prograii naie: seft.c to find efficiency of synchronous •achines M 

It works with only a single sachine tl 

lainO 

{ 

FILE tfopenO, tfp; 

double r, jrnl, jr, js, Is, Ir, dcore, ds, dr, g, w, 1, xs, eaf, i2rr, 

va, ph, pe, i2r, vol, wt, effcy, ew, xsl, xs5, xs7, xsal, kstB], 
siv, riv, wt_iron, find_sivi), pieff, pijs, pirpi, dhp, rpi, iinpwr, 
ke, kv, freq; 

extern double swfO, sqrtO; 
int e=0, f, p, i; 
char infiletH]; 

for (i=l; i < 8; i+=2) It hartonic winding factors XI 

ks[il=swf (i); 

} 

X5l=ks[l]tks[l]; 

xs5=(kst5]tks[5]/25); 

X57=(k5[7]Jksl7]/49); 
xsal = (5UstPltds/l8); 



printf (“\nCalculate5 efficiency of a single iotor,\n"); 
while (e != 1) 

{ 

f = 0; 

printf ("Khat is the na«e of the input file? *); 
scanfl'Zs’, infile); 
fp = fopenlinf ile, ’r‘); 

fscanfifp, ’U‘, lip); It input nuaber of pole pairs t/ 

fscanfifp, 'Xlf, liainpHr); 

■inpur 746.0; It noH in uatts tl 

fscanfifp, ‘Zlf, bite); 

fscanfifp, ‘Zlf, bkv); 

fscanfifp, ‘Zlf, brpa); 

fscanfifp, ‘Zlf, bjs); 

fscanfifp, ‘Zlf", bfreq); 

fscanfifp, "Zlf", bw); 

fscanfifp, *Zlf, br); 

fscanfifp, ‘Zlf, bg); 

fscanfifp, *Zlf‘, bdcore); 

fscanfifp, 'Zlf", bds); 

fscanfifp, "Zlf’, bdr); 

fscanfifp, *Zlf“, bis); 

fscanfifp, 'Zlf, blr); 

fscanfifp, 'Zlf, bvol); 

fscanfifp, 'Zlf, but); 



154 






I 

§; 

i 




fscanfUp, 'Xlf, iph); 
fscanfifp, ■XH“, ipe); 
fscanf(fp, ’IH', ii2r); 
fscanfifp, "SH", iva); 
fscanf(fp, 'Iir, Iteffcy); 
f5can^(fp, 'inf*, lew); 
fscanf(^p, ’Iir, &1); 
tscanfifp, "Xlt”, lijr); 
fscan^(fp, ■Ilf, lijrnDi 
fclose(fp); 

siv = tind_siv(l, r, g, ds, dcore, Is); /> stator iron volute t/ 

riv = lIPIlrKr - 2ldrtlr); /I rotor iron volute 1/ 

wtjron = Dtlriv + siv); /I iron weight only 1/ 

while It != 1) 

{ 

printt('\nNhat is the sustained speed tachine horsepower? ’); 
scantCZlf , 4dhp); 
dhp 1= 746.0; 

printtfNhat is the sustained speed tachine rpt? '); 
scanfCIlf, iptrpt); 

freq = ptrpttp/60.0; IX tax electrical frequency M 

ptjs = jsidhp/tinpwr; H PH stator current 1/ 

x5=(HUIjstSSFIIr+g)l(r+g)ldslPUl5 > (12l(xsl + xs5 + xs7)/(Pllgtp) + xsaD) 

/(12lrtBRIi(sn]); /I p.u. synch itpedance 1/ 

eaf=sqrt(l + xsixs + 2txsl0.6); /I 0.6 is sin(r), pwr factor angle, 
eaf is p.u. internal voltage at full load 1/ 
jr=eaftjrnl; /I jr full load, linear with eaf 1/ 

ph = 31.86225»BETAIfreqlBRl»HCUwt.iron/D; 

/t hysteresis loss in watts, uses iron weight of tachine M 
pe = (106236. 9tNUIBSATIBSATtTltTltfreqlfreqlwt_iron)/(RH0ID); 

/I eddy current loss in watts, uses iron weight of tachine 1/ 
i2r = 2.0tSSFtdslCRH0IPIt (1 + 2.3094IPll(r+g)»CP/p) Ijsljsllr+gltls; 

/) stator copper loss in watts t/ 
i2rr = 2.0tRSFIdrlrllrljrtjrtCRH0IPU (I + 2.3094>PI«r/p); 

/I rotor excitation losses, 1-12-87 1/ 
pteff = dhp/ldhp + ph + pe + i2r + i2rr); 
printf('\n Sustained speed efficiency is Ilf, pteff); 



printf f \nHhat is the endurance speed tachine horsepower? '); 
scanfCIlf , Itdhp); 
dhp 1= 746.0; 

printffUhat is the endurance speed tachine rpt? '); 
scanfCIlf , iptrpt); 

ptjs = jsidhp/tinpwr; IX PM stator current 1/ 

freq = ptrpttp/60.0; IX tax electrical frequency 1/ 

xs=(HUIjslSSFIlr+g)l(r+g)tdslPllls I (12t(xsl + xs5 + xs7)/(Pltglp) + xsaD) 

/(12trtBRtlcsCl]); IX p.u. synch itpedance X! 

eaf=sqrt(l + xsixs + 2lxsl0.6); IX 0.6 is sin(t), pwr factor angle, 
eaf is p.u. internal voltage at full load X! 
jr=eafljrnl; IX jr full load, linear with eaf 1/ 

ph = 31.86225tBETAtfreqlBRltHCllwt,iron/D; 
pe = (106236. 9INUIBSATtBSATITllTltfreqlfreqlwtJron)/(RH0ID): 
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i2r = 2.0tSSFtd5tCRH0tPII (1 + 2.3094IPII(r+9)tCP/p) tj5tj5t(r+Q)ns; 
i2rr = 2.0tRSFtdrtrtlrtjrljrlCRH0IPIt il + 2.3094tPMr/p); 
paeff = dhp/(dhp + ph pe + i2r + i2rr); 
pnntf("\n Endurance efficiency is Xlf", pieff); 

printf ("\nSaiie lachine'^ *); 

5canf(”Xd% «cf); 
if (f == 0) 

continue; 
else if (f == 2) 

{ 

e = 1; 
break; 

} 

} It end of f-loop tl 

) It end of e-loop XI 

} It end of tain prograt tl 

double find^sivd, r, g, ds, dcore, Is) 
double 1, r, g, ds, dcore, Is; 

{ 

double one, tNo, three, four; 

one = (r+g-^ds+dcoreXIr+g+ds+dcore) - (r+g)l(r+g); 
two = 2IPII(r+g)td5lls; 

three = (r+g+ds+dcore)l(r+g+ds+dcore) - (r+g+ds)t(r+g+ds); 

four = IKPItone - two) + PIt4t(r+g)lthree; 
return(four); 

} 
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Appendix C. Permanent magnet machine 



The equivalent circuit for a permanent magnet machine 
is almost identical to that of a synchronous machine. The 
only difference is the source of the internal voltage, which 
develops the field flux wave that interacts with the arma- 
ture flux wave. The field flux wave is a result of per- 
manent magnets built into the rotor to develop magnetic 
poles . 

A typical magnetic circuit, combined with Ampere’s Law, 

shows 




Hmlm + Hg g - Ni 



Figure 40. Typical magnetic circuit 



We use a constitutive law of B = PH and assume that any 
steel has P = ®. If flux leakage is about zero, BmAm = BgAg 
since flux is solenoidal. Then, 

Ni - Hgg Po Hg Ag 

Im and Am = 

Hm Bm 

The magnet volume is 

„ , , (Ni - Hgg) Po Hg Ag 

Hm Bm 

Minimum magnet volume occurs when the magnet’s maximum 
energy product (MEP) , HmBm , is a maximum. If current is 
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zero and only magnitudes are used, 



Vm 



P o Hm Bm 



The load line of the magnet is developed as 



Im 


g 


Bm 




Am 


Po 


Hm 




Bm 


Po 


Im 


Ag 


Hm 


g 


Am 



The permanance coefficient, Pc, is the slope of the load 
line. On a magnet diagram 




Figure 41. Magnet operating point diagram 

A good algorithm for machine design is to adjust the 
geometric dimensions to operate at the MEP, on the load 
line. If operation at MEP is assumed, the needed slope is 
determined, the dimensions are randomly generated, PoHg is 
calculated, and the design is maximized for PoHg and mini- 
mized for magnet volume, then a search technique has been 
delineated . 

Magnets may "overhang" the active length at either end 
to account for manufacturing imprecision and to permit a 
smaller armature diameter. This overhang affects the 
developed flux. 
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Figure 42. Permanent magnet machine diagram 

The flux per pole, 0, with overhang and a given Ja and Jm, 
is the same as would exist in a configuration in which 
Ja’ =Jm’ =Ja+ND and there were no overhang. N is a non- 
linear parameter promulgated in graph form by Ireland [33] . 

A good power fit for N is 

(S’ Jm 

N = 0,38558 ( )o. 7oei3 

D 

<S Jm 

where = overhang and 0 i overhang i 0.34. Then, 

D 

Ja (S' Jm 

+ 

D D 

Y = and 0with = 0w/o Y 




D 



This f lux-with-overhang is applied to the problem as would 
be the usual flux. VHiat is the usual flux? A permanent 
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:: 










magnet hysteresis diagram shows the residual flux density to 




The magnet does not operate at that point, but rather 
on the load line. Also, there is not magnet material at 
every point along the circumference of the rotor. A Fourier 
series is a good way to find the f lux-wlthout-overhang. 

Using a developed rotor, 




Figure 44. Magnet material on a developed rotor 



PM = Bd 

wr p wr 

Ir = = 

iTr/p TT r 



® n TT X 2TTr 

B(x) = J An cos( ), T = 

n = i T p 
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PT 


niTx 


An = 


B(x) 


cos ( 


- ) dx for n = 1, 2, 3 . . 


T 


J 0 




T 


B(x) = 


Bd 


for 


0 i X i (iTr Ir/p) 




0 


for 


(TTr Ir/p) i X 1 (TTr/p) 




Bd 


for 


(TTr/p) i X i (TTr/p) ( 1 + lr ) 




0 


for 


(TTr/p) ( 1+lr ) i X i (2TTr/p) 



After integration and manipulation, 



2Bd IT Ir IT Ir 

An = [ sin( ) - cos( ) + 1] 

TT 2 2 



The equation represents only the fundamental term of the 
Fourier series. The f lux-without-overhang is 



2Bd 

Br = 

TT 



sin( 



TT Ir 




TT Ir xp 

cos( ) + 1] cos( ) 

2 2r 



As is usual in steady-state analysis, the magnitude is used. 

The next quantity to find is Bd. If the magnet 
operates on the load line, the operating point flux is 



Brem 

Bd = where Pr is the relative 

1 + pr /Pc 

reversible permeability of the magnet and Brem and Pc are as 
previously defined. With this relative permeability, the 
magnet length is Im = g Pr . The magnetic machine can now be 
specified. An end view with dimensions is given in Chapter 
Four . 
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Table 40. Listing of permanent magnet machine design program 



linclude “stdio.h" 
♦include "def.h" 



/I prograt nate: pii.c for pertanent tagnet tachines 1/ 

long int seed; /t start point for randoi nuiber generator Xf 

double b[lOl3[271, h[in[27], ksC8]; 

IX b is "best" array, h is "hold" array, ks are winding factors X! 

•ainO 

{ 

double design_point(), rnd_»alk(), swfO, ke, kv, ki, tinpwr, 
stepsize, randoiO, absO, freq, rpi; 
int p, iteration, i, j, best, print^outO, loops; 

FILE tfopenO, tfp; 



printf ("\nReading input data frot PMrt.DAT . . ."); 



fp=fopen("pia.dat","r"); fX input seed for randoi nuabers X! 
fscanf (fp,"!5d-,^seed); 
fscanf(Fp, "ltd", ip); 
fscanfifp, "Xlf", iainpwr); 

•inpwrt=746.0; 
fscanf(fp, "Xlf", ike); 
fscanf(fp, "Xlf", ikv); 
fscanfifp, "Xlf", ika); 
fscanf(fp, "Xlf", irpi); 
fcloselfp); 



iX input nuaber of pole pairs X! 

IX input aachine power, derived fa ASSET XI 
IX convert to watts XI 
IX CERs for Effective Weight 1/ 



IX aachine aax shaft rpa XI 



printf ("\nHow aany loops do you want? "); 
scanf("Xd", iloops); 

printf ("\n\nDoing prograa calculations . . .\n"); 

for (i=l; i < 8; i+=2) IX haraonic winding factors XI 

{ 

ksli]=swf (i); 

} 



freq=rpalp/60.0; IX aax electrical frequency XI 

IX MAIN BODY OF THE PROBRAH XI 



for (i=i; i <= loops; ++i) 

{ 

printf("\ni=Zd", i); 

stepsize=0.1; 

iteration=0; 

design_point(ainpwr, p, ke, kv, ka, freq); 

IX put stuff in the hold array XI 

while (iteration <= 10) 

{ 

rnd_walk(ainpwr, p, stepsize, ke, kv, ka, freq); 
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/I stagger around t! 

It index to best EH of the lot t/ 



best=0; 

for (j=l; j<=10; ++j) 

if (hljlflBl < hlbestinSl) 

best=j; /( find the best iachine 1/ 

if (absUhlOKlB] - htbest][lB])/(h[0]nB]+. 001)1 < 0.005) 

/t stall iiprovetent in EH 1/ 

/I 0.001 talies care of div by zero »/ 

{ 

stepsize/=2.0; 

++iteration; 

) 

else /t transfers best to 0 position 1/ 

{ 

for (j=l; j <= 26; ++j) 

htOHj] = hlbestltj]; 

) 

> 

for (j=l; j <= 26; ++j) 

b[i]Cj]=h[best][j]; /t lieep the best iachine M 

} 

best=l; 

for (i=l; i <= loops; ++i) 

if (blilClB] < btbestHlBl) 

best=i; /t find and keep the best of the best 1/ 

•inpHr/=746.0; /I turn back into hp M 

print.outibest, p, iinpHr, ke, kv, ka, rpi); /I output to disk file M 

fp=fopen(‘pii.dat*,*«*) ; It output seed »/ 

fprintf(fp,*Jd", seed); 

fprintf (fp,'\nXd*, p); 

fprintf(fp,*\nXlf, linpwr); 

fprintf(fp,'\nllf, ke); 

fprintf (fp,"\nXlf kv); 

fprintf(fp,'\nXlf, ki) ; 

fprintf (fp,“\nZlf’, rpi); 

fclose(fp); 

} 

n END OF MAIN PR06RAH; ALL THAT FOLLOH ARE FUNCTIONS 1/ 

double design_point(iinpHr, p, ke, kv, ki, freq) 

!% deteriines a randoi design point 1/ 

double linpHr, ke, kv, ki, freq; 
int p; 

{ 

double r, js, Is, Ir, dcore, ds, g, h, 1, xs, va, ph, itiag, 

pe, i2r, vol, Ht, effcy, ei, xsl, xsS, xs7, xsal, li, voliag, 
bd, find_lr(), c», siv, scv, riv, loa, doa, find.sivO, it_iron; 
extern double sqrtO, randoiO; 
int c=0, d=0; 

Hhile (d 1= 1) ( 
while (c != 1) 
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{ 

r=randoi()IRnAX; 

if ((2<PIIrlfreq/p) < KAX.TIP.SPEED) !\ check rotor tip speed 1/ 
break; 

} 



a=2tPllfreq; 




bd=BREH/(l + nUR/PC); 


It 


lr=find.lr(bd); 


It 


ahile (c != 1) 


It 



{ 

g=randoil)t(0.1tr - GHIN) 
if (g > 0) 

break; 



} 



+ 



air gap flux, operating point, froa lagnet 
characteristics 1/ 
find lagnet slot factor t/ 

gap diaension 1/ 

6H1N; 



la=gtnUR; 



/t radial length of aagnet, see PH-8, 12/15/86 1/ 



ls=randoa()t0.5 + 0.25; 
js^randoaDIdSHAX; 
dcore=(bdt(r+la))/(BSATIp); 
ds=randGa()t0.9ldcore; 



/I stator slot factor t/ 

/t full load stator current density t/ 
/t back iron depth t/ 

/I slot depth < 90Z of body depth 1/ 



xsl^ksfntksni; 

xs5=(kst5]tks[53/25); 

xs7=(ks[71tkst7]/fl9); 

xsal=l5tlstPltds/18); 

X5=(HUtjstSSFt(r^gila)t(r+g^la)tdstPltls 

t U2t(xsl + xs5 + xs7)/(Pltglp) + xsal))/(12t(r+la)tbdtksmi; 

It p.u. synch iapedance t/ 



if (xs > 3.0) 

continue; /t don’t aant xs too big 1/ 



else 



++d; 



/I ay escape hatch — aission coaplete t/ 



l=(ainpHrlp)/(2lbdtHtks[ntlrHa)ISSFtdslPll(r+g+la)tlstjslPF); 

It active length t! 



va=2tPltlr^la)lltMtbdtksnitjstSSFt(r^g^la)tdstls/p; It VA rating t! 



siv = find_siv(l, r, g, ds, dcore. Is, la); It stator iron voluae 1/ 
riv = llPltrtr; It rotor iron voluae 1/ 
scv = 2tPlt(r+g+la)tdstls«ll + 2.3094tPlt(r+g+la)tCP/p); 



ca = scvtDCU; 

loa = 1 + 4t(r+g+la); 

doa = 2t(r+g+la+ds+dcore); 

vol = V0LALLI(loatPltdoatdoa/4); 

at = HTALHIca + DKBRNGStriv + siv)); 
at_iron = DKriv + siv); 



It stator copper voluae t! 
It total copper aeight tl 
It length-over-all tl 
It over-all-diaaeter tl 
It aachine envelope voluae 1/ 

It aachine aeight in kg tl 
It iron aeight only tl 



ph = 31.86225tBETAtfreqlBRllHCltat_iron/D; 
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/t hysteresis loss in Halts, uses iron Height of lachine t/ 



pe = (106236.9tNUtBSATIBSATtTllTUfreqtfreqlHt_iron)/(RH0tD); 

/t eddy current toss in Halts, uses iron Height of lachine \! 

voliag=2<Pltr>ltlitlr; /t aagnet voluie M 

Httag=voltagtDKA6; /I lagnet Height t/ 

i2r = 2.0tSSFIdslCRH0tPU (I + 2.3094IPU(r+g+Ii)tCP/p) IjstjsKr+g+lams; 

/t stator copper loss in natts 1/ 
/» revised 1-12-87 M 



effcy=(ainpHr)/(ainpHr + ph + pe + i2r); 



eH=wt + kitHtaag * ketd-effcyl * kvtivol * kitvoliag); 

/t Effective Height 1/ 



h[0][l]=jsi 

h[0H6]=dcore; 

h[0][ll]=vol; 

h[0][l6]=vai 

h[0][20]=Htiag,‘ 

h[0][23]=0.0; 

h[0][26]=doa; 

n this 



h[0H2]=freq; 

h[0H7]=ds; 

h[0ni2]=Ht) 

h[0H17]=effcy; 

h[0][21]=volaag; 



h[0][31=Hj 

h[0][81=li; 

h[0][131=ph; 

h[0H18]=eH; 



h[0H4]=r; 

htOU9]=ls; 

h[0][14]=pe; 

h[0][19]=l; 



h[0][221=xs; 



/( this one lakes overhang=0 1/ h[0H241=bd; 



h[0][51=g) 

h[0][l0]=lr; 

h[0][l5]=i2r; 



h[0H25]=loa; 



section just changed all the variables in the ‘hold* array t/ 



return; 

)} 



double rnd_Halk(iinpHr, p, stepsize, ke, k.v, ki, freq) 

/t nalks about design_point 10 tiies 1/ 
double stepsize, linpnr, ke, kv, ki, freq; 
int p; 

{ 

double r, js. Is, Ir, dcore, ds, g, h, 1, xs, va, ph, pe, i2r, 

vol, Ht, effcy, en, xsl, xs5, xs7, xsal, li, voliag, Htiag, bd, 
find_lr(), find_br(), br, overhang, Irat, ch, siv, scv, riv, loa, 
doa, find_siv(), nt.iron; 
extern double sqrtO, randoiO; 
int i=l, j, check=0, ccj; 

Hhile (i <= 10) 

{ 

/t read in the nalk around the design point t/ 

js=h[0imt(l * stepsizetfrandoiO - O.SII; 
if (js > ISHAXI 

js = JSMAX; /t reset to liiit M 

H=2tPUfreq; 

r=h[0][41t(l * stepsizeKrandoiO - 0.5)); 
if ((H»r/p) > HAX.TIP.SPEED) 
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continue; /I go to next try if violated 1/ 

g^hlOHSltll -1 stepsizetlrandoiO - 0.5)1; 
if (g < 6HIN) 

g=GHlN; /I reset to the luit 1/ 



overhang=randoi 010.34; 



l=h[0][19]; 
l«=gtnUR; 

bd=BREH/(l + nUR/PC); 
lr=find.lr(bd); 
ls=h[0][9]l(l + stepsizeKrandoiO - 0.5)); 
if (Is > 0.75) 

15=0.75; 
if (Is < 0.25) 

ls=0.25; 



/t lagic nuiber *0.34' frot the book 
on perianent lagnets by Jaaes Ireland. 

The effect of overhang is to concentrate air 
gap flux, or reduce leakage, t/ 

/t new la, based on the new g M 
It aagnet operating point flux density tl 
It find nea Ir tl 



It reset to the liait tl 



for (j=l; j <= 10; ++j) It aagnet char convergence loop tl 

{ 

lrat=l/(24(r+g+la)); 
br=bdlfind_br (overhang, Irat); 

It find effect of overhang on air gap flux tl 
dcore=(brl(r+la))/(BSRTtp); It aost efficient use of iron tl 
ds=h[0][7]tdcore/h[0][&]; It keep thea in the saae ratio 

as upon exit fa design_pt 1/ 



xsl=ks[ims[]]; 
xs5=(ks[5]»ks[5]/25); 
xs7=(ks[7]lks[7]/49l; 
xsal = (541stFilds/18); 

xs=(NUtj34SSFI(r+g+la)<(r+g+la)tdslPItls I (12t(xsl * xs5 ♦ xs7)/(PUglp) 

* xsal))/(12t(rala)tbrtksm); 

It p.u. synch iapedance tl 

l=(ainpHrtp)/(2tbrtHtks[l]((rHa)tSSFtdstPlt(rtgaia)tl5tjstPF); 

It active length tl 

) It end of aagnet char convergence loop tl 

if (xs > 3.0) 

{ 

++check; 

if (check > 25) It can’t close on decent xs tl 

{ 

for (ccj=i; ccj <= 10; ++ccj) 
htccjK 181=10000000.0; 

It aake this randoa walk undesirable tl 
printfC burp*); 

break; It go to next design point tl 

) 

continue; It can’t have xs too big tl 

} 
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va=2tPIt (r+l«)tltHtbrlks[ntj5lSSFI(r+g+li)td5tl5/p; 



/I VA rating «/ 



51V = find_siv(l, r, g, ds, dcore, Is, li); !% stator iron voluie t/ 
riv = llPIlrtr; /I rotor iron voluie 1/ 
scv = 2tPIt(r+g+li)<ds»l5t(l t 2.3094tPIt(r+g+l#)»CP/p); 

/t stator copper volute t/ 



CN = scvtDCU; 

loa = 1 + 4l(r+g+lt); 

doa = 2t(r+g+lt+ds+dcore); 

vol = VOLALLtUoatPItdoatdoa/4); 



/t total copper Height t/ 
/I length-over-all t/ 

/I over-all-diaieter t/ 

/t lachine envelope volute t/ 



Ht = WTALLKch + DKBRNGStriv + siv))j 
Htiron = Dtlriv + siv); 



/t tachine Height in kg 1/ 
/t iron Height only t/ 



ph = 3l.86225tBETAIfreqtBRltHCllHt.iron/D| 

/t hysteresis loss in Hatts, uses iron Height of tachine t/ 



pe = (106236. 9»«UIBSATtBSAmitTltfreqltreq«Ht.iron)/(RH0ID); 

/t eddy current loss in natts, uses iron Height of tachine t/ 

voliag=2IPItrtH(1.0+(overhangt(2t(r+g+lt))))tlttlrj /I tagnet volute t/ 

Httag=voltaglDnA6; /t tagnet Height t/ 

i2r = 2.0tSSFIdslCRH0IPlt (I + 2.3094«Pll(r+g+lt)tCP/p) tjsIjsKr+g+ltllls; 

/t stator copper loss in natts t/ 

/t revised 1-12-87 1/ 



effcy=(tinpHr)/(tinpHr + ph + pe + i2r); 

en=Ht + ktlntsag + ket(l-effcy) + kvKvol + ktivoltag); 

/t Effective Height t/ 



h[i][ll=js; 

h[i][6]=dcore; 

h[iHll]=vol; 

hti]tl6]=va; 



h[i][2]=freq; 

h[iH7]=ds) 

h[iHl2]=Ht; 

hti]tl7]=effcy; 



h[i][3]=H; 

h[i][8]=lt; 

h[i][13]=ph; 

h[iH18]=eH; 



h[i][4]=r; 

h(iH9]=ls; 

h[i][14]=pe; 

h[iH19]=l) 



h[i}120]=Httag; hli H21]=voltag; 
h[i][23}=overhang; 



h[i][22]:xs: 



h[i][3]=g; 

htim0]=lr; 

h[i][15]=i2r 



h[iH24]=br; h[iH25]=loa; hti]t26]=doa 
/I this section just changed all the variables in the 'hold' array 1/ 



++i; H go to the next hliHl 1/ 

check=0; 

} 



return; 

} 



print_out(best, p, tinpnr, ke, kv, kt, rpt) 
int best, p; 

double tinpnr, ke, kv, kt, rpt; 

{ 

char outfilelMl; 
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FILE Ifpo, IfopenO; 
int i ; 

printf (“\nHhat is the nate of the file Hhere you want the output? "); 
scanfCZs', outtile); 

fpo=fopen(outtile, “w'); 

fprintf (fpoi'Zd’, p); 
fprintf (fpo,*\nXlf ", •inpwrl; 
fprinttffpo, ’\nllf", ke); 
fprintf (fpo,’\nXlf", kv); 
fprintf(fpo,’\nIlf", ki); 
fprintf(fpo,"\nZlf, rpi); 
for (i=l; i <= 26; ++i) 

fprintflfpo, ■XnXlf, btbestlli]); 
fprintfffpo, ’\n'); 
fclose(fpo); 

} 



double find_lr(br) /I bracketing t/ 

double br; 

{ 

extern double cosO, sinO; 
int c=0; 

double lr=0.5, top=1.0, tbr=0.0, absO, bti=0.0; 

while (c!=l) 

{ 

tbr=(2tBREH< (sin(0.5IPltlr) - co5(0.6IPIllr) + l.OH/PI; 
if (abs((br-tbr)/br) <= O.OOl) /I check for convergence 1/ 
break; 
if (br > tbr) 

{ 

bte=lr; 

lr=bti + (top- bti)/2.0; 

} 

else /» br < tbr 1/ 

{ 

top=lr; 

lr=bti + (top-bti)/2.0; 

} 

if (abs((top-lr)/top) <= 0.0005) /I no convergence t/ 

{ 

printfCXnllf, Ir); 
abortCno solution for Ir’); 

} 

} 

return(lr); 

) 
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double find br (overhang, irat) 
double overhang, Irat; 

{ 

extern double poH(); 

double expon=0. 706133501, factor=0. 385576858, enn, yyy; 

i\ enn is Ireland's N, yyy is Ireland's Y XI 

enn = factor I pow(overhang, expon); 
yyy = (Irat + overhang)/(lrat ^ enn); 
return (yyy) ; 

} 

double find_siv(l, r, g, ds, dcore, Is, li) 

double 1, r, g, ds, dcore, ds. Is, li; 

{ 

double one, tno, three, four; 

one = (r+g+ln-ds+dcore)l(r+g+li+ds+dcore) - (r+g+li) Kr'^g+li); 
two = 2tPU(r+g+li)ldslls; 

three = (r+g+li+ds+dcore)t(r+g+li+ds+dcore) - (r+g+li+ds)l(r+g+li+ds) ; 

four = IKPltone - two) + PU4t(r+g+U) tthree; 
return(four); 

} 
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Table 41. Listing of permanent magnet efficiency program 



iinclude "stdio.h" 
linclude "def.h“ 



It program name: peH.c to find efficiency of permanent magnet machines tl 
It works with a single machine tl 

eainO 

{ 

FILE tfopenO, tfp; 

double r, js, Is, Ir, dcore, ds, dr, g, w, 1, xs, It, km, sgeff=1.0, 

va, ph, pe, i2r, vol, wt, effcy, ew, xsl, xs5, xs7, xsai, ks[83, 
siv, riv, wt^iron, find^sivO, pmeff, pgeff, teff, pgjs, pmjs, 



pmrpm, geff=0.98, pceff=0,99, gr, dhp, rpi 
wtmag, volmag, overhang, bd; 
extern double s«f(), sqrtO; 
int e=0, f, p, i; 
char infilelH]; 



iinpwr, ke, kv, freq, 



siv = find_siv(l, r, g, ds, dcore, Is, Im); 
riv = llPlirlr; 
wt_iron = Dt(riv + siv); 



It stator iron volume tl 
It rotor iron volume tl 
It iron weight only tl 



printf (“\nCalculates efficiency for a single machine. \n“); 
while (e 1= 1) 

{ 

f = 0; 

printf ("What is the name of the input file? "); 
scanf(“Zs", infile); 
fp = fopenlinf ile, V); 

fscanfifp, "’Ad", Skp); It input number of pole pairs tl 

fscanfifp, "Xlf", Ifminpwr); 
minpwr t= 746.0; It now in watts tl 

fscanfifp, 
fscanfifp, 
fscanf ifp, 
fscanfifp. 



fscanfifp, 
fscanf ifp, 
fscanf ifp, 
fscanf ifp, 
fscanf ifp, 
fscanf ifp, 



fscanfifp, 
fscanfifp, 
fscanfifp, 
fscanfifp, 
fscanfifp, 
fscanf ifp, 
fscanfifp, 



"Ilf, 


Ue); 


“Ilf, 


tikv); 


■Ilf, 


tiki); 


■Ilf, 


&rpi) ; 


■Ilf, 


fejs); 


■Ilf, 


ifreq); 


■Ilf, 


(ih); 


“Ilf, 


tir); 


"Ilf, 


liQ); 


■Ilf, 


&dcore) 


■Ilf, 


tids); 


■Ilf, 


Ua); 


■Ilf, 


Ils); 


■Ilf, 


Ur); 


■Ilf, 


&vol); 


"Ilf, 


&Ht); 


■Ilf, 


&ph); 


■Ilf, 


Jipe); 


■Ilf, 


M2r); 


■Ilf, 


&va); 
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fscanfdp, 'Xir, ieffcy); 
fscanfUp, lew); 

fscanfifp, 'Xlf, U); 
fscanfifp, "Xlf, iHtsag); 
fscanfifp, ’ilf’, ivoliag); 
fscanfifp, ’XIP, ixs); 
fscanfifp, *Xlf", tiQverhang); 
fscanfifp, 'Xlf', Ibd); 
fcloselfp); 

Hhile If != 1) 

{ 

printf rVnNhat is the sustained speed eachine horsepower? '); 
scanfl'Xlf, idhp); 
dhp 1= 746.0; 

printfl'Hhat is the sustained speed tachine rpt? *>; 
scanfl'Xlf *, ipirpt); 

freq = parpstp/60.0; /t lax electrical frequency M 

pijs = jsJdhp/iinpwr; It PN stator current 1/ 

ph = 31.86225tBETAIfreqtBRUHCUwtJron/D; 

It hysteresis loss in watts, uses iron weight of lachine t! 
pe = 1106236. 9INUIBSATtBSATtmTllfreqlfreqtwt.iron)/IRH0ID); 

It eddy current loss in watts, uses iron weight of tachine t/ 
i2r = 2.0ISSRdslCRH0IPll II + 2.3094IPItlr+g+lt)tCP/p) tjstjstlr+g+lt)41s; 

It stator copper loss in watts t! 

pteff = dhp/ldhp + ph + pe + i2r); 

printfl'\n Sustained efficiency is Xlf, pteff); 

printf l'\n\nWhat is the endurance speed eachine horsepower? '); 
scanfl'Xlf, yhp); 
dhp 1= 746.0; 

printfl'Hhat is the endurance speed tachine rpt? '); 
scanfl'Xlf, Jiptrpt); 

ptjs = jsidhp/tinpwr; It PK stator current t/ 

freq = pirpttp/60.0; It tax electrical frequency t! 

ph = 31.86225tBETAIfreqtBRUHCltwtJron/D; 

pe = 1106236. 9INUt6SATtBSATtTUmfreqtfreqtwt.iron)/IRH0tD); 

i2r = 2.0ISSFtdstCRHOIPU II + 2.3094IPItlr+g+rt)tCP/p) ljs»js»lr+g+lt)tls; 

pteff = dhp/ldhp + ph + pe + i2r); 

printfl'\n Endurance efficiency is Xlf, pteff); 

printf I'VnSaie tachine? '); 
scanfl'Xd', &f); 
if If == 0) 

continue; 
else if If == 2) 

( 

e = 1; 
break; 



) 

} 

} 



It end of f-loop >/ 

It end of e-loop 1/ 

It end of tain prograt t! 
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i 




double find^sivd, r, g, ds, dcore, Is, li) 
double 1, r, g, ds, dcore, Is, 1 a; 

{ 

double one, two, three, four; 

one = (r+g+lt+ds+dcore)t{r+g+lA+ds+dcore) - (r+g+li)t(r+g+lA); 
two = 2IPII(r+9+li)ld5lls; 

three = (r+g+li+ds+dcore)l(r+g+li+ds+dcore) - (r+g+la+ds)t (r+g+li+ds) ; 

four = U(PIIone - two) + PII4l(r+g+li)Uhree; 
return(four); 

} 



172 



Appendix D. Induction machines 



The equivalent circuit for an induction machine is 
shown in the figure. 



R, X, 




Figure 45. Induction machine equivalent circuit 

The stator leakage reactance, XI, is the sum of the 
stator differential leakage reactance (Xsd) , the stator slot 
leakage reactance (Xss), and the stator end turn leakage 
reactance (Xse). In turn, Xsd is the sum of the belt and 
zig-zag reactances. 

The rotor leakage reactance, X2, is the sum of the 
rotor differential leakage reactance (Xrd), the rotor slot 
leakage reactance (Xrs), and the skew leakage (Xskew) . Xm 
is the magnetizing reactance and Rc is the core resistance. 
R1 and R2 are the stator and rotor resistances, respec- 
tively. 

A derivation of the properties of electric machines, 
using Ampere’s Law and the constitutive relation earlier 
postulated, yields the result that 

® Po 6Ns^lr 

flux = E I ks^ cos(wt) 

n = 6kl;l g ffn2p2 

k=0 
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From this important stator inductances are taken. 

Belt leakage reactance (Xbelt) is the sum of the reac- 
tances due to phase-belt harmonics of an "infinite" slot 
winding. In most machines, the most important harmonics 
present are the fifth and seventh, as the third is canceled 
in balanced operation. Then the belt reactance is 



SPowNs^lr ks5 ks7 

Xbelt = [( y + ( y ] 

TT g p2 5 7 



This is the harmonic form of the fundamental mutual reac- 
tance of Appendix B. The winding factor is again the 
product of the pitch and breadth winding factors. 

Zig-zag reactance is leakage due to all the air gap 
harmonics that would be produced if the winding had one slot 
per pole per phase. For a phase belt of one slot, with each 
slot carrying the same current and equally separated in time 
and space phase, the zig-zag reactance alone would be 
present. Belt leakage occurs because phase belts are ac- 
tually several slots wide. Zig-zag reactance has harmonic 
orders higher than seven, with the same form as Xbelt. No 
even or triplen harmonics will be present. 

The fundamental harmonic of the flux yields the mag- 
netizing reactance, Xs, which can be viewed as that required 
to "energize" the air gap. 

6 P 0 wNs 1 rks^ 

Xm = 

TT g p2 

Figure 46 shows a typical stator slot. The stator slot 
leakage reactance, summing the self and mutual reactances, 
is 



Xss = 



ns2 ^ w2 



ds ns ^ 

+ )(---- + PNp) 

2 ws 6 
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where ns is the number of stator slots and Np is the coil 
throw, in slots. Np = (CP ns)/(6 p) , a result available by 
manipulation. w2 




\AIS 



Figure 46. Stator slot diagram 



Then, 



3 w Mo 1 Ns^ d2 ds 

Xss = + )(1 + CP) 

ns w2 2 ws 

and 

(1 - Is) 2tt (r+g) 

ns = 

ws 

where we have traded the need for the knowledge of the num- 
ber of stator slots for the need to know the width of an in- 
dividual slot. A reasonable relationship between slot 
dimensions is d2 = ds and w2 = ws . Then, 



3_Mo_w_l Ns2 ds (1_+ CP) 

" Is 2tt (r+g) 

The stator end turn leakage reactance may be estimated 
by treating the two end regions as a single helically shaped 
winding. If the active region of the machine is ignored and 
the helix given air core properties, the inductance can be 
found from standard sources. 



Xse 



-16 w n Mo iS sin^ (9w/2) No^ 
le p2 



Ip’( 



piT r 
le 



)Kp’ ( 



piT r 
le 



) 



where le is the combined length of both end windings. 
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Presume the winding radius to be r, with helix pitch 
©w = tt/3 and 

le = 2nr tan(Tr/3) Np /ns 

Ip ’ and Kp ’ are the first derivatives with respect to their 
arguments of the hyperbolic Bessel functions Ip and Kp . 

When the three phases are summed, a multiplier of 1.5 will 
be realized. Finally, using the previous results for Np , 



Xse 



BlwMoNs^r ptrr pirr 

Ip»( )Kp>( 

PTT^ CP tan(TT/3) le le 



Stator resistance is R1 = (rho It Ns)/(As/N8), as in 
Appendix B, leading to 



^ cp/p)_ 

SSF ds TT (r+g) Is 

When the actual calculation is performed, 1 is not known. A 
guess-and-iterate scheme is used. Iteration continues until 
convergence on 1 is achieved. 

Rotor resistance uses a similar scheme, but the 
presence of rotor bars and end rings instead of turns 
changes it somewhat. A model of a rotor bar is below. 



w2 




Figure 47. Typical rotor bar configuration 



Induction motor transformer models provide a way to 
find rotor resistances and inductances. The flux density 
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produced by the stator and rotor is 



Po 

B = (Fs + Fr ) 

g 

The stator mmf is Fs = -j 



3 Ns ks ws 
Ti 

TT p 



where Ti is stator 



current and ws is slip frequency. Rotor mmf is 
Ws r* Br 

Fr = -j 

p2 Zs 

where Zs is the rotor surface impedance. The air gap volt- 
age, or voltage across Xm, is defined in terms of the flux 
density and rotor mmf as 



Vag = -j 



w 

Ws 



2 p 1 Ns ks 



Fr 



If rotor mmf is now identified with rotor current referred 
to the stator winding, 



3 Ns ks 



Fr = -j l2 , where I 2 is rotor current. Then, 

TT p 

Zs w 6 1 Ns2 ks2 
Vag = I 2 

Ws TT r 



Separating Zs into its real and reactive parts and using a 
rotor surface model to describe the relation between rotor 
electric field amplitude and rotor surface current yields 



12 1 Ns2 ks2 

R2 = rslot 

nr 

12 1 Ns2 ks2 

Xrs = xsiot 

nr 

SPowlNs^rks* 1 1 

Xrd = Z {( )2 + (— )2 } 

p+k^t^nr p-k*nr 



TT g 
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If magnetic diffusion is ignored, end ring resistance 
can be calculated by comparing losses in the rings and slot. 
The ratio of current densities is found by the ratio of the 
areas. This is then squared and multiplied by the ratio of 
volumes. When summed, 

12 1 Ns^ ks^ nr r wr 

R2 rslot [1 + ] 

nr TT 1 ler p* 

rslot = CRHO/(dr wr), nr is the number of rotor bars 
and ler is the end ring length, approximated as 
ler = 2 tt (r - wr/4 - ds/2). The rotor bar width, wr, is 
found by specifying nr and observing that (nr wr) = (2 tt Ir) . 

Rotor skew leakage arises when the rotor slots are 
skewed angularly along the axial length to prevent rotor 
cogging. Then, flux does not fully link the bars. When the 
effect is integrated over the rotor, it is seen that 

Xskew = Xm [1 - (2 sin( skew/2 )/skew)2 ] with the amount 
of skew measured in radians. When typical values of skew 
are input, we see that (Xskew/Xm) ~ 0.5%. This is a negli- 
gible effect and will be ignored. 

From the previous rotor bar model, it is seen that 
xsiot = w Po (d2/w2 + dr/2wr). Assume that d2 = wr/4 
and w2 = wr/4. Then xslot =wPo (1+ dr/2wr) . 

Fitzgerald et al [23] state that only small errors 
result if Rc is omitted. Therefore, the core branch may be 
omitted . 



Once the components of the equivalent circuit have been 
calculated, the designer must turn to power and torque con- 
siderations. The internal mechanical power of the machine 
is 



P 



(l_-_slip) 

slip 



3 122 R2 



The air gap voltage has been previously defined. The ter- 
minal voltage, Vt, may be found from Vag by means of a volt- 
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age divider and a Thevenin equivalent circuit developed. 




Figure 48. Induction machine Thevenin equivalent circuit 



R1 + j(Xl + Xm) 

Rel + jXel = (R1 + jXl) in parallel with jXm. Then operat- 
ing point torque is 

1 3 Vla2 (R2/slip) 

T = 

w (Rel + R2/slip)2 + (Xel + X2)2 

Torque is a maximum when the power delivered to (R2/slip) is 
a maximum. By matching load and Thevenin impedances, the 
power is a maximum and a slip-at-maximum- torque is found 



R2_ 

” (Rel2 + (Xel + X2)2 )o*5 

and the corresponding torque is 

1 1.5 Via* 

Tmstx ~ — 

" w Rel + (Rel2 + (Xel + X2)2 )0- = 

Typical induction motors have the ratio between maximum and 
operating point torque as 



which is used to find the machine active length to use in 
the circuit component calculations. The equation for rated 
torque can be manipulated to yield a quadratic expression 
for operating point slip, or a Newton’s method convergence 
can be used to find operating point, slip. Convergence on 
slip and active length through Newton’s method is used to 
generate the machines of this thesis. 

Finally, rotor copper losses are (1 - slip)P; stator 
copper losses are found using Ii and R1 . 
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Table 42. Listing of induction machine design program 



tinclude 'stdio.h* 
linclude 'def.h* 



/t proQrat pane: ind.c for induction aachines, 4/19/87 1/ 

long int seed; /I start point for randoa nuaber generator 1/ 

double bC261[331, h[ll][331, xs[421; 

/I b is "best" array, h is "hold" array, xs are winding factors 1/ 

aainO 

{ 

double design^pointO, rnd^walkO, swft), ke, kv, ainpwr, 
stepsize, randoaO, absO, freq, rpa, ks[431; 
int p, iteration, i, j, best, print^outO, loops, flag; 

FILE tfopenO, Ifp; 

printf ("\nReading input data froa IND.DAT . . ."); 



fp=fopen("ind.dat","r"); 
fscanf(fp,"Xd",!tseed); 
fscanflfp, "Xd", ip); 
fscanflfp, "Ilf", iainpHr); 
ainpwrt=746.0; 
fscanflfp, "Xlf", Lke); 
fscanflfp, "Xlf, bkv); 
fscanflfp, "Xlf", brpal; 
fcloselfp); 



/t input seed for randoa nuabers M 

It input nuaber of pole pairs t/ 

It input aachine power, derived fa ASSET 1/ 
It convert to watts tl 
It CERs for Effective Weight 1/ 

It aachine aax shaft rpa tl 



printf rVnHow aany loops do you want? "); 
scanfl"):d", Moops); 

printf r\n\nDoing prograa calculations . . .\n"); 



for li=l; i < 42; i+=2) It haraonic winding factors 1/ 

{ 

kslil=swf li); 

xsfil = IksliltkstiD/liti); 

} 

freq=rpatp/60.0; It aax electrical frequency tl 



It HAIN BODY OF THE PROGRAM tl 



for li=l; i <= loops; ++i) 

{ 

printfl"\nZd", i); 
stepsize=0. 1; 
iteration=0; 
flag - 0; 

design_pointlainpwr, p, ke, kv, freq, Lflag); 

It put stuff in the hold array tl 

if Iflag == 1) 

{ 

hlOltlBl = 10000000.0; 
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best = 0; 
} 



Hhile ((iteration <= 101 ik (flag 1= D) 

r 



rnd_Halk(ainpHr, p, stepsize, ke, kv, freq); 

It stagger around tl 

best=0; It index to best EM of the lot 1/ 

for (j=l; j<=10; ++j) 

{ 

if (hCj][18] < 1.01 
continue; 

if (h[j1(18] < h[best][18]) 

best^j; It find the best lachine tl 

) 

if (abs((h[0H18] - h[be5t]I18])/htO]ll8]l < 0.0051 
It saall iiproveient in EM 1/ 

{ 

stepsize/=2.0; 

<-+iteration; 

) 

else It transfers best to 0 position tl 

{ 

for (j=l; j <= 32; ++j) 

h[0][j] = hfbestlfj]; 

} 



} 



for (j=l; j<= 32; ++j) 

b(i][j]=h[best][j]; It keep the best lachine tl 

} 



best=l; 

for (i=l; i <= loops; ++i) 

{ 

if (b(iH18] < 1.0) 
continue; 

if (b[iH18] < b[best](18]) 

best=i; It find and keep the best of the best tl 

} 

■inpwr/=746.0; 

print_out(best, p, linpHr, ke, kv, rpa); 

fp=fopen(’ind.dat’, V); 
fprintf (fp,')ld', seed); 
fprintf(fp,'\nld‘, p); 
fprintf (fp,‘\nnf*, iinpNr); 
fprintf(fp,‘\nllf, ke); 
fprintf (fp,*\n)llf*, kv); 
fprintf(fp,’\nXlf, rpi); 
fclose(fp); 

} 

It END OF RAIN PR06RAN; ALL THAT FOLLOM ARE FUNCTIONS tl 
double design_point(iinpNr, p, ke, kv, freq, fla) 



It turn back into hp tl 
It output to disk file tl 
It output seed tl 
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/t deteraines a randoi design point t/ 



double oinpHr, ke, kv, freq; 
int p, tfla; 

{ 



double r, Is, Ir, dcore, ds, dr, Nr, g, n, nb, 1, i2rr, 
rl, r2, xbelt, tu, xss, xse, xl, xrd, xrs, 
xa, vag, vt, via, rel, xel, taax, radical, sb, ila, ilb, 
x2, va, ph, pe, i2r, vol, at, effcy, ea, le, tainpar, 
sa, sc, ca, siv, scv, riv, rev, loa, doa, r2saax, ttaax; 
double ler, fl, f2, f3, slip, saax, jr, js, 
il, i2, xslot, at_iron; 

extern double sqrti), randoaO, cost), sinO, coshi), sinhO, vratO, 
floor!), besipO, beskpO, tanO, find_siv(), safO; 
int cslip=0, nr=7l, ccr=l, ccg, ccs=l, k; 

a=2IPl(freq; /I synchronous frequency in rad per sec </ 

aa=a/p; /t aech angular velocity in rad per sec tl 

1 = 10.0; It initial guess on length 1/ 

taax = ainpar/laatPSI); It pull-out torque tl 



ahile leslip != 1) 
{ 



r=randoa(IIRMAl(; 

if (!2iPIIrlfreq/p) < MAX_TIP_SPEEI)) It check rotor tip speed tl 
break; 

} 



lr=randoa()l0.5 + 0.25; 
dr=randoal)lr/3.0; 
dcore=IBR»r)/(BSATtp); 
ds=randoa I ) tO. 9tdcore; 
9=randoal)t(0. Ur - BMIN) 
if (g < S«IN) 

g = eniN; 

ls=randoa()t0.S * 0.25; 



It rotor slot factor tl 

It slot no deeper than 33Z of rotor radius tl 

It back iron depth tl 

It slot depth < 90Z of body depth to start tl 
+ GrtlN; 



It stator slot factor tl 



fl=0. 

f2=0. 

f3=0. 



for (k=l; k < 100; ++k) 
{ 



fl = (l./(p + kinr)) t II. /(p + kinr)); 
f2 = (l./(p - kinr)) I (l./(p - kinr)); 
f3 += fl + f2; 

} 

ahile leer <= 5) It start convergence loop on r tl 

{ 

if Ices > 5) It ae’ve had trouble aith jr tl 

{ 



r 1= jr/JSNAX; 

if Ir > INAX.TlP.SPEEDIp/a)) 
r = «AX_TlP_SPEEDIp/a; 
if Ir > RHAX) 



r = RMAX; 
Ir 1= jr/JSHAX; 
if Ur > 0.75) 

Ir = 0.75; 
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) 



wr = (BlPIllrlir - dr/2. ) )/(4tnr + Plllr); /I rotor slot width M 
ler = 2tPlt(r - wr/4.0 - dr/2.0); It end ring length 1/ 

CCS = 1; 

while (CCS <= 51 /I start convergence loop on ds and dr 1/ 

{ 

ccg = 1) 

Mhile (ccg <= iO) /t start convergence loop on length t/ 

{ 

rl = (6tCRH0l(I + 2.3094»Plt(r+g)ICP/p))/(SSFIdstPI»lr+g)IIs); 
r2 = (6txsmiCRH0tlMl + (2trtrtIr)/(lerlplplll/(Pltrtlrtdr); 
xbelt = (6tHlHHUIrl(xst5] + xst7])/(PIIglptpl); 
xrz = (4tHtItHUIrt(x5(in + xstl3] + xsCl?] + X5tl9] + xs[23] + 
xs[29] + xs(31] t xst37] + xs[41])l/(PItg»plp); 
xss = (3tHUt»IHd5l(l+CP))/(Ist2tPItlr+q)); 
le = (PllrlCPItan(PI/3))/(3lp); 

xse = (27t»tll«Ulrlrl(-beskp(p, ptPllr/leltbesip(p, ptPItr/le)))/ 
(PItletptp); 
xl = xbelt + X 2 Z + xss + xse; 
xrd = l4tnutNmxsmtrtf3)/IPllg); 

xrs = (4txs[ntHlHHUt(Hr+0.5ldr))/(lrlPll(r - (dr/2.0) - (wr/8.0111 
x2 = xrd + xrs; 

XI = (6tHUtHtllrtxsm)/(Pltgtptp); 

vag = 2trl6RtwtltsMt(l)/p; /t air gap voltage M 

rel = (xBtx«trll/(rltrl + (xl+xB)l(xl+xill; 

/I thev equiv resistance 1/ 

xel = (xBllrllrl + xltxl + xltx«l )/(rltrl + (xl+xi)t(xl+xil); 

/t thev equiv inductance 1/ 

r2s«ax = sqrt(reltrel + (xel+x2)l(xel+x2)l ; /I r2 at siax t/ 
vt = vag/vrat(r2siax, x2, xl, xi, ri); 

/t teriinal voltage at siax M 
via = (vtlx«l/sqrt((xl + xalKxl + xi) * (rltrl)); 

It thevenin equivalent voltage t/ 
ttaax = (l.Stvlatvial/(Hit(reUsqrt(reitreH(xeHx2)t(xei+x2l ))); 

It test iaxituB torque t/ 
if (absUtiax-ttiaxl/ttax) <= 0.005) 

break; It we have convergence 1/ 

1 1= tiax/ttiax; It reset 1, reiterate 1/ 

++ccg; 

if (ccg > 10) 

{ 

Ifla = 1; 

printf('\n flag set on length'); 
return; 

) 

} It end convergence loop on length 1/ 

siax = r2/sqrt(rellrel + (xel+x2)t(xel+x2)i; 

slip = siax/3.0; It starting point for slip converge tl 

cslip = 1; 

while (cslip <= 20) 

{ It start convergence loop on slip tl 
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tiinpwr = (3lvlalvla»r2/slip)/((rel+r2/slipl» (rel+r2/slip) + 
(xel+x2IKxel+i(2)); 

i-f (abs( (ainpKr-tainpwrI/ainpKr) <= 0.005) 

break; It « have convergence 1/ 

slip 1= tiinpHr/ninpwr; 

++C5lip; 

} It end convergence loop on slip 1/ 



i2 = vagtsqrt((r2/slip)t(r2/slip) + x2tx2)/((r2/slip)l(r2/slip) + x2lx21; 
jr = (3ti2)/(PUlrldrlrlRSF); It rotor current density M 

i< (jr > JSMfllO 
{ 

dr jr/JSHA)!; 
if (dr > (r/3)i 

dr = r/3; It reset to li»it t! 

} 

ila = (r2/slip)l(xi+x2) - x2lr2/slip; 
ilb = x2t(xi+x2) + (r2/slipll(r2/slip); 

il = vagtsqrt(ilatila ^ i Ibtilb)/(xi((x2tx2 * (r2/slip)t(r2/slip) ) ) ; 
js = (3.0lil)/lPll(r+g)lds»ls»SSF); 
if (js > JSHA)!) 

ds 1= js/JSHAX; 

if ((js <= JSMAX) kk (jr <= JSMAX)) 
break; 

++ccs; 

} It end convergence loop on ds and dr tl 

if (CCS <= 5) 
break; 

++ccr; 

if (ccr > 5) 

r 

L 

Ifla = 1; 

printf(*\n flag set on r*); 
return; 

) 

) It end convergence loop on r tl 

It calculations tl 



snax = r2/sqrt(rellrel + (xel+x2)t(xel+x2)); 



siv = find_siv(l, r, g, ds, dcore, Is); It 
riv = UPltrKr - 2»drllr); It 
scv = 2IPlt(r+g)lds»lst(l + 2.3094tPlt(r+g)ICP/p); 

It 

rev = wrtdrKnrIl + 2tler); It 
cw = (rev + scvlIDCU; It 
loa = 1 + 4t(r+g); It 
doa = 2l(r+g+ds+dcore); It 
vol = VOLALLMloalPIIdoaldoa/4); It 
(*t = NTALLKcm + DKBRNBSIriv + sivl); It 
Nt_iron = DKriv + siv) ; It 
ph'= 31.8i225IBETAIfreqlBRllHCll«tJron/D; It 
■achine tl 



stator iron volute tl 
rotor iron volute tl 

stator copper volute tl 
rotor copper volute tl 
total copper Height tl 
length-over-all tl 
over-all -diateter tl 
tachine envelope volute 1/ 
tachine Height in kg tl 
iron Height only tl 

hysteresis loss in natts, uses iron Height of 



pe = (106236. 9»NUIBSATIBSAT»TllTllfreqlfreqlHtJron)/(RH0tD); 
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/t eddy current loss in xatts, uses iron weight of aachine t/ 
i2r=3.0ti Ui Url; !% stator copper loss in watts 1/ 

i2rr = sliptainpwr; /f rotor copper loss in watts 1/ 

vt = sqrt((vag+illrl)t(vag+iUrl) + (i UxDKi Utl) ) ; IX rated voltage X! 
va = 3tvt>il; IX VA rating at terminals X! 

effcy=(«inp«r)/(iinpwr + ph + pe + i2r + i2rr); ew=wt + keXl-effcy) + kvtvol; 
fective weight X! 



h[0Hll=il; h[0K2]=freq; h[0H3]=w; 
h[0][4]=dcore; h[0K7]=ds; ' h[0K8]=dr; 

h[0Hll]=vol; hC0H121=wt; h[0H13]=ph| 

h[0][16]=va; hC0]tl7]=effcy; h[0]C18]=ew; 
h[0H21]=saax; h[0H221=tiax; h[0H231=vti 
h[0K24]=rl; h[0H27]=xl; h[0K28]=xa; 

h[0K31]=loa; h[0H321=doa; 

IX this section just changed all the 



h[0][41=r; 

h[0][91=ls; 

h[0}[141=pe; 

htOK191=l; 

h[0H24]=slip; 

h[0H291=x2; 



h[0H51=g; 

hCOltlOlMr; 

h[0K151=i2r; 

ht0H20]=i2rr; 

h[0H251=vag; 

ht0]t301=r2; 



variables in the 'hold' array X! 



IX Ef- 



return; 

} 



double rnd_walk(ainpwr, p, stepsize, ke, kv, freq) 

IX walks about designjoint 10 tiaes XI 
double stepsize, ainpwr, ke, kv, freq; 
int Pi 

{ 

double r, Is, Ir, dcore, ds, dr, wr, g, w, wa, 1, i2rr, ila, ilb, 
rl, r2, xbelt, xzz, xss, xse, xl, xrd, xrs, r2saax, ttaax, 
xa, vag, vt, via, rel, xel, taax, radical, sb, tainpwr, 
x2, va, ph, pe, i2r, vol, wt, effcy, ew, le, cw, siv, scv, 
sa, sc, riv, rev, loa, doa, jr, js, 
ler, fl, f2, f3, slip, saax, il, i2, xslot, wt_iron; 
extern double sqrtO, randoaO, cosO, sinO, coshO, sinhO, vratO, 
besipO, beskpO, tanO, find.sivO, swfO; 
int d, i=l, ccg, ccr, ccs, nr=71, k, cslip; 
fl=0.; 
f2=0.; 
f3=0.; 

for (k=l; k < 100; ++k) 

{ 

fl = (l./(p + ktnr)) » (l./(p + ktnr)); 
f2 = U./(p - ktnr)) I (l./(p - ktnr)); 
f3 += fl + f2; 

) 

while (i <= 10) IX ten steps around the design point XI 

{ 

IX read in the walk around the design point XI 
w=2tPItfreq; 

r=h[0H4]t(l + stepsizetlrandoaO - 0.5)); 
if ((wtr/p) > HAX.TIP.SPEED) 

continue; IX go to next i-loop if violated XI 

g=h[0]C5]t(l + steps! zet (randoaO - 0.5)); 
if (g < eniN) 
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q=GI1IN; /I reset to the li»it »/ 

dcore=(BRtr)/(BSfiTtp); /I lost efficient use of iron </ 

ds=ht0]t7]|(] + stepsizetirandoiO - 0.5)); 

dr=htO]tB]l(l + EtepEi:e»(rando#() - 0.5)); 

if (dr > r/3.0) 

dr = r/3.0; IX reset to the liiit 1/ 

ls=h[0][9]l(l + stepsizeKrandoiO - 0.5)); 
if (Is > 0.75) 

ls=0.75; IX reset to the hut 1/ 

if (Is < 0.25) 
ls=0.25; 

lr=h[0]tl0]<(l * stepsize»(randoi() - 0.5)); 
if (Ir > 0.75) 

lr=0.75; IX reset to the liiit X! 

if (Ir < 0.25) 
lr=0.25; 

IX lech angular velocity in rad per sec X! 
IX length starting point XI 
IX pull-out torque XI 



Hi=M/p; 

1 = h[0][19l; 
tiax = ht0]t221; 



CCS = 1; 
ccr = 1; 

while (ccr <= 5) 

{ 

if (ccs > 5) 

( 

r t= jr/JSHAX; 

if (r > (HAX_TlP_SPEED»p/w)) 
r = «AXJlP.SPEEDtp/w; 
if (r > RMAX) 

r = RHAX; 

Ir t= jr/JSHAX; 
if (Ir > 0.75) 

Ir = 0.75; 

} 

wr = (BtPltlrKr - dr/2.))/(4tnr + Pltlr); IX rotor slot width XI 
ler = 2tPU(r - wr/4.0 - dr/2.0); IX end ring length XI 

CCS = 1; 

while (ccs <= 5) IX start convergence loop on ds and drt/ 

{ 

ccg = 1; 

while (ccg <= 10) IX start convergence loop on length 1/ 

{ 

rl = (6ICRH0M1 + 2.3094IPlt(r+g)tCP/p))/(SSFtdstPlt(r+g)ns); 
r2 = (OtxsmtCRHOmd + (2trtr»lr)/(lertptp)))/(Pl»rtlrtdr); 
xbelt = (AtwtUHUtrt(xs[5] + xst7])/(Pl»gtptp)); 
xiz = (6twtltHU»r«(xs[lll + xs[13] + xstl71 + xsllRl + xst231 + 
xs(291 + xst3ll + XSC37] + xs[4ll))/(PUglptp); 
xss = (3tHUtwtlldst(l+CP))/(lst2IPll(r+g)); 
le = (PltrlCP»tan(PI/3))/(3lp); 

xse = (27twmMUtrtrt(-besl(p(p, ptPltr/le)tbesip(p, ptPltr/le)))/ 
(PUletptp); 
xl = xbelt + X 22 + xss + xse; 



IX first tiie thru, no adj in r XI 

IX start convergence loop on r XI 
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xrd = (6mU>Nlllx5[l]>rtf3)/(Pllg); 

xrs = (6lx5[mHnif1Ut(wr+0.5ldr))/(lrlPII(r - (dr/2.0) - (wr/8.0))); 
x2 = xrd + xrs; 

xn = (6IMUIwtllrlx5[l3)/(PUglptp); 

vag = 2>rtBRtwtltswf (D/p; IX air gap voltage X! 

rel = (xfilxalrD/irltrl + (xUxs)l(xl+xa)); 

IX thev equiv resistance 1/ 

xel = (xtKrllrl + xllxl + xllxi) )/(rUrl + (xl+xi)t(xl+xi)); 

!% thev equiv inductance 1/ 

r2siax = sqrt(rellrel + (xel+x2)t(xel+x2)); /t r2 at saax t! 
vt = vag/vrat(r2siax, x2, xl, x«, rl); 

!% terainal voltage at saax %! 
via = (vt>xt)/sqrt((xl + xa)t(xl + xa) + (rllrl)); 

!% thevenin equivalent voltage 1/ 
ttaax = (1.5tvlatvla) / (Hit (rel+sqrt(reltrel+(xel+x2) l(xeUx2) ) ) ) ; 

/t test iaxiaua torque t/ 
if (abs((taax-ttiax)/taax) <= 0.005) 

{ 

ccq = 1; 

break; /I we have convergence 1/ 

} 

1 1= tia*/tt»ax; It reset I, reiterate 1/ 

++ccg; 

) It end convergence loop on length 1/ 

if (CCS > 5) 
break; 

siax = r2/sqrt(rellrel + (xel+*2)t(xel+)(2)); 
slip = s*ax/3.0; 

It starting point for slip converge 1/ 

cslip = 1; 
while (cslip <= 20) 

( It start convergence loop on slip t! 

tiinpwr = (3lvlalvlalr2/slipl/((rel+r2/slip)l(rel+r2/slip) + 
(*el+x2)l(xel+x2)); 

if (abs((ainpHr-tainpHrl/ainpMr) <= 0.005) 

break; It we have convergence 1/ 

slip t= tainpMr/iinpHr; 

++cslip; It no lore than 20 tries t! 

} It end convergence loop on slip t! 

i2 = vaglsqrt((r2/slip)l(r2/slip) + *2t*2)/((r2/slip)»(r2/slip) + x2»x2); 
jr = (3li21/(PltlrldrtrJRSF); It rotor current density t! 

if (jr > JSMfiX) 

{ 

dr »= jr/JSNAX; 
if (dr > (r/3)) 

dr = r/3; It reset to liiit t! 

) 

ila = (r2/slip)t(xi+x2) - x2lr2/slip; 
ilb = x2t(xi+x2) + (r2/slip)l(r2/slip); 

il = vaqlsqrt(ilatila + ilblilb)/(xiJ(x2Jx2 + (r2/slip)<(r2/stip))); 
js = (3.0till/(Pl»(r+g)ldstlstSSF); 
if (js > JSHAX) 
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ds M js/JSMflK; 

i< Ujs <= JSNAXI W (jr <= JSHA)()I 
break; 

++CCS; 

} /} end convergence loop on ds and dr M 

il ((CCS <= 5) i: (ccg > 10)) 
break; 

++ccr; 

} /} end convergence loop on r t/ 



if ((ccg >10) !! (ccs >5) I! (cslip >20) !! (ccr > 5)) 

{ It no design convergence >/ 

for (d=i; d <=10; t+d) 

hCdlElB] = 10000000.0; 
printfCburp •); 
break; 

} 



It calculations t! 

siax = r2/sqrt(reltrel + (xel+i(2)»(xel+x2)); 



siv = find_siv(l, r, g, ds, dcore, Is); It 
riv = ItPlirtlr - 2tdrtlr); It 
Ecv = 2tPIt(r+9)tdstlst(l + 2.3094tPIt(r+gHCP/p); 

It 

rev = Hrldrllnrll + 2ller); It 
CN = (rev + scvlIDCU; It 
loa = 1 + 4l(r+g)j It 
doa = 2t(r+g+ds+dcore); It 
vol = VOLALLI(loatPUdoatdoa/4); It 
«t = HTALLKch + DKBRNBSIriv + siv)); It 
Ht_iron = DKriv + siv); It 
ph’= 31.B6225tB£TAtfreqlBRllHCUi(t_iron/D; It 
■achine tl 



stator iron volute tl 
rotor iron volute tl 

stator copper volute tl 
rotor copper volute tl 
total copper Height tl 
length-over-all tl 
over-all-diateter tl 
tachine envelope volute tl 
tachine Height in kg tl 
iron Height only tl 

hysteresis loss in natts, uses iron Height of 



pe = (106236.9tNUtBSATtBSATtTltTltfreqtfreq»Ht.iron)/(RH0»D); 

It eddy current loss in natts, uses iron Height of tachine tl 
i2r = 3.0tiltiltrl; It stator copper loss in natts tl 

i2rr = slipttinpHr; It rotor copper loss in natts tl 

vt = sqrt((vag+iltrlU(va 9 +iltrll + (illxl)t(iltxi)); It rated voltage tl 
va = 3>vtlil; It VA rating at the tertinals tl 

effcy=(tinpHr)/(tinpnr + ph + pe + i2r + i2rr); en=Ht + ket(l-effcy) + kvtvol; It Ef- 

fective Height tl 



h[iHl]=il; 

h[il[61=dcore; 

h[i][ll]=vol; 

h[i][161=va; 

h[iH21]=stax; 

hti][26]=rl; 

h[i][31]=loa; 

It this 
++i; 

) 



h[il[21=freq; h[il[3]=n; 
h[i][71=ds; hti][B]=dr; 
h[i][12]=nt; h[i][131=ph; 

h[i]tl71=effcy; h[i][18]=eH; 
h[i][22]=ttax; h[il[231=vt; 
hti][271=xl; h[i][20]=xt; 

h[i][321=doa; 

section just changed all the 



h[i][4]=r; 

h[i][91=ls; 

h[i][l4]=pe; 

h[i)[191=l; 

h[i][241=slip; 

h[i][29)=x2; 

variables in the 



h[i)[5)=g; 

h[i][10]=lr; 

h[i][15]=i2r; 

h[i][201=i2rr; 

h[i][2S]=vag; 

h[iK30]=r2; 

hold" array tl 



return; 

} 
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pnnt_out (best, p, ainpwr, ke, kv, rpa) 
int best, p; 

double ainpMr, ke, kv, rpo; 

{ 

char outfileCU]; 

FILE Ifpo, IfopenO; 
int i; 

printf (“\nWhat is the nate of the file where you want the output? "I 
scanf CXs", outfile); 

fpo=fopen(outfile, "h"); 



fprintf(fpo,-Zd% p); 
fprintf (fpo,"\nXlf“, linpHr); 
fprintf(fpo,-\nXlf, ke); 
fprintf(fpo,“\nZlf“, kv); 
fprintf(fpo,"\nZlf“, rpi); 
for (i=l; i <= 32 ; ++i) 

fprintf (fpo,"\nXB“,b[best](i]); 
fprintf(fpo, "\n"); 
fclose(fpo); 



double find_siv(l, r, 9, ds, dcore, Is) 
double 1, r, 9, ds, dcore. Is; 

/ 

k 

double one, two, three, four; 

one = (r+9+ds+dcore)l(r+9+ds+dcore) - (r+9)l(r+9); 
two = 2IFII(r+9)ldslls; 

three = (r+9+ds+dcore)l(r+9+ds+dcore) - (r+9+ds)l(r+9+ds); 

four = IKPIlone - two) + PII4l(r+9)lthree; 
return(four); 

} 

double vrat(r2siax, x2, xl, xi, rl) ft ratio of va9/vt ti 
double r2saax, x2, xl, xa, rl; 

{ 

double a, b, c, d, vrat; 
extern double sqrtO; 

a = r2sfaxl(xi+xl) + rll(x2+xi); 
b = rllr2siax - xilx2 - xlt(x2+xi); 
c = alr2s§ax - blx2; 
d = blr2saax ^ alx2; 



vrat = (x«lsqrt(clc + dld))/(ala + bib); 
return(vrat);} 
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Table 43. Listing o-f induction efficiency program 



linclude ■stdio.h" 
linclude "def.h" 



/I progras naae; ieff.c to find efficiency of induction lachines </ 

/I works with a single nachine )/ 
iainO 
{ 

FILE tfopenO, Ifp; 

double r, Is, Ir, dcore, ds, dr, g, w, 1, i2rr, va, ph, pe, i2r, vol, 

wt, effcy, ew, siv, riv, wt^iron, find^sivO, pteff, pgeff, teff, 
parpa, geff=0.98, pceff=0.99, gr, dhp, rpt, linpwr, ke, kv, freq, 
sgeff=l.O, wr, w®, rl, r2, xl, xi, vag, vt, via, rel, xel, radical, 
x2, sa, sb, sc, slip, slipl, slip2, saax, il, 12, tiax; 
extern double sqrtO; 
int e=0, f, p, i; 
char infilefUl; 

printf (*‘\nCalculates efficiency for a single lachine.Xn"); 



while (e != 1) 



{ 

f -- 0; 

printfl'Khat is the naie of the input file? ’); 
scanfl’Zs", infile); 
fp = fopenlinfile, *r"); 

fscanflfp, ’Zd”, lip); /t input nuiber of pole pairs 1/ 

fscanflfp, "Xlf", litinpwr); 

iinpMr t- 746.0; /t now in watts tl 

fscanflfp, "Xlf, like); 

fscanflfp, "Xlf, likv); 

fscanflfp, "Xlf, lirpi); 

fscanflfp, "Xlf, Ml); 

fscanflfp, "Xlf", lifreq); 

fscanflfp, "Xlf, tiH); 

fscanflfp, ‘Xlf, lir); 

fscanflfp, 'Xlf, tig); 

fscanflfp, ‘Xlf, Mcore); 

fscanflfp, "Xlf“, lids); 

fscanflfp, ‘Xlf, Mr); 

fscanflfp, "Xlf, Ills); 

fscanflfp, ‘Xlf, Ur); 

fscanflfp, "Xlf, livol); 

fscanflfp, ‘Xlf, Mt); 

fscanflfp, ‘Xlf, liph); 

fscanflfp, ‘Xlf, ipe); 

fscanflfp, ’Xlf, U2r); 

fscanflfp, "Xlf, liva); 

fscanflfp, "Xlf, Mffcy); 

fscanflfp, "Xlf", Itew); 

fscanflfp, "Xlf, U); 

fscanflfp, "Xlf, U2rr); 

fscanflfp, "Xlf", (siax); 
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fscanHfp, ‘Xlf, SiUaxI; 
fscanfifp, “XH"' Wt); 
fscanfifp, *Xir, islip); 

^5can^(fp, "XH“, ivag); 
fscanfifp, "Xlf, Irl); 
fscanfifp, *Xlf“, kl); 
fscanfifp, 'Ilf, ki); 
fscanfifp, ’Xlf“, k2)j 
fscanfifp, "21f, k2); 
fcloseifp); 

while if != 1) 

{ 

printf i‘\nNhat is the sustained speed aachine horsepower? '); 
scanfi'ZH", idhp); 
dhp t= 746.0; 

printf f°Khat is the sustained speed aachine rpa? '); 
scanfi'Xlf, (parpa); 

freq = parpalp/60.0; /t PH frequency tl 

via = ivtlxal/sqrtiixl + xaltixl + xa) + frllrDI; 

It thevinin equivalent voltage 4/ 

rel = ixalxalrll/irUrl + ixl+xa)tlxl+xa)); 

/t thev equiv resistance t/ 

xel = ixatirltrl + xUxl + xltxa))/irl4rl + ixl+xa)tlxl+xa)); 

/t thev equiv inductance 4/ 

sa = rel4rel + ixel+x2)41xel+x2); /4 pieces of slip quadratic 4/ 
sb = 24rel4r2 - i3414vla4vlalr2/dhp) ; 
sc = r24r2; 

radical = sb4sb - 44sa4sc; 
if iradical < 0.0) 

aborti*\nGot a negative radical in the slip eqn.'); 
slipl = i-sb + sqrt(radical))/i24sa); 
slip? = i-sb - sqrtiradical))/i24sa); 
if iislipl < 0.0) islip2 < 0.0)) ' 

aborti‘\nTwo negative slips.*); 

It now we will use the saallest positive slip 4/ 
else if Kslipl < 0.0) U (slip? > 0.0)) 
slip = slip2; 

else if ((slip! > 0.0) kk (slip? < 0.0)) 
slip = slipl; 

else 

slip = (slipl > slip?) ? slip? : slipl; 

It tslip=ain(slipl, slip?) 4/ 

if (slip > saax) 

abort(*\n slip is aore than saax.*); 

siv = find_siv(l, r, q, ds, dcore, Is); It stator iron voluae 4/ 

riv = 14F'14r4(r - ?4dr41r); It rotor iron voluae tl 

wt_iron = 04 (riv + siv); It iron weight only 4/ 

ph = 31.86?254BETA4freq4BR14HC14wt_iron/D; 

It hysteresis loss in watts, uses iron weight of aachine tl 
pe = (106236. 94NU4BSAT4BSAT4T14T14freq4freq4wt.iron)/(RH04D); 

It eddy current loss in watts, uses iron weight of aachine tl 
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i2 = 5qrt((5liptdhp)/(3tr2t(l-5lip))); 
il = i2 + vag/xa; 



/I load current >/ 

/t i2 plus lagnetizing current 1/ 



i2r=3.0lilliltrl; /I stator copper loss in watts t/ 

i2rr = sliptdhp; IX rotor copper loss in watts 1/ 

pteff = dhp/ldhp + ph + pe + i2r + i2rr); 
printf(“\n Sustained efHciency is Xlf“, pieff); 

printf ("\n\nWhat is the endurance speed lachine horsepower? "); 
scanfrXlf", 4dhp); 
dhp t= 746.0; 

printf ("lilhat is the endurance speed tachine rpt? M; 
scanfi'Zir, ipirpa); 

freq = parpatp/60.0; /I PH frequency X! 

sb = 2treltr2 ~ (3tUvlatvlatr2/dhp); 
radical = sbtsb - 4tsatsc; 
if (radical < 0.0) 

abort C\n6ot a negative radical in the slip eqn.M; 
slipl = (-sb + sqrt(radical))/(2tsa); 
slip2 = (-sb - sqrt(radical))/(2tsa); 
if ((slipl < 0.0) (slip2 < 0.0)) ’ 

abort(“\nTwo negative slips."); 
else if ((slipl < 0.0) U (slip2 > 0.0)) 
slip = slip2; 

else if ((slipl > 0.0) kk (slip2 < 0.0)) 
slip = slipl; 

else 

slip = (slipl > slip2) ? slip? : slipl; 
if (slip > saax) 

abortC\n slip is aore than saax."); 
ph = 3l.86225IBETA<freqlBRltHCUwt.iron/D; 
pe = (106236.9INUtBSATtBSAmmi<freq<freqtwtJron)/(RH0tD); 
i2 = sqrt((sliptdhp)/(3tr2t(l-slip))); It load current XI 

il = i2 + vag/xa; IX i2 plus aagnetizing current XI 

i2r=3.0tiltillrl; IX stator copper loss in watts XI 

i2rr = sliptdhp; IX rotor copper loss in watts XI 

paeff = dhp/ (dhp + ph + pe + i2r + i2rr); 
printf(“\n Endurance efficiency is Xlf", paeff); 

printf ("\nSaBe aachine? M; 

5canf(-):d", S(f); 
if (f == 0) 

continue; 
else if (f == 2) 

{ 

e = 1; 
break; 

} 

} IX end of f-loop XI 

} IX end of e-loop XI 

} IX end of aain prograa XI 
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double find_siv(l, r, g, ds, dcore, Is) 



double 1, r, g, ds, dccre. Is; 

{ 

double one, two, three, four; 

one = (r+g+ds+dcore)t(r+ 9 +ds+dcore) - (r+g)l(r+g); 
tHO = 2IPII(r+g)ldsUs; 

three = (r+g+ds+dcorelKr+g+ds+dcore) - (r+g+dslKr+g+ds); 

lour = IKPIlone - two) + PIt4l(r+g)tthree; 
return (lour ) 5 
} 
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Appendix E. Weight and volume algorithms 



The ability to correctly characterize a new naval ship 
technology depends in part on the ability to calculate the 
weight and volume associated with that technology. Algo- 
rithms for this purpose were taken from a variety of 
sources, including the Advanced Surface Ship Evaluation Tool 
(ASSET) theory manuals. The ASSET algorithms are the result 
of data analysis for naval ships that have been constructed, 
as well as studies for other ship designs. 

The Ship Work Breakdown Structure (SWBS) categorizes 
all ship weights. The general categories are: 



WlOO 


Ship structures 


W200 


Propulsion plant 


W300 


Electric generation plant 


W400 


Command and control equipment 


W500 


Auxiliaries and distributed 
systems 


W600 


Outfit and furnishings 


W700 


Ship armciment 


WFOO 


Ship Loads 



Within W200 are several sub-groups that pertain to 
electric propulsion. They are: 



W235 
W235. 1 
W235.2 
W235. 3 

W235.4 
W235.5 
W241. 1 

W242 

W243 

W244 

W245 



Electric propulsion devices 
Propulsion motors 
Propulsion generators 
Transmission lines and 
propulsion cables 
Cooling systems 
Switchgear 

Locked- train -double -reduct ion 
reduction gears 
Propulsion clutches and 
couplings 
Shafting 

Propulsion shaft bearings 
Propeller weight 
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Jipil 






ASSET allows SWBS groups to be adjusted in weight, 
which allows technology sensitivity analyses, such as this 
thesis, to be performed. Only a few of the W235 sub-groups 
need to be calculated outside of ASSET and adjusted within 
ASSET. These are W235.1, W235.2, W235.4, W235.5, W241.1, 
and W243. The W235.1 and W235.2 weights and volumes were 
calculated as part of the machine design. The rest of the 
needed W235 weights were calculated in "wt.c", a copy of 
which follows. These algorithms are the U. S. Navy 
standard, and were verified against actual ships and com- 
ponents . 

Shafting and transmission line weights are dependent on 
motor and generator positions within a ship. Their weights 
and volumes were calculated from ASSET equations, using the 
layouts of the baseline and variant ships. 

Where no volume equation was found in the ASSET 
documentation, or where the result of such an equation was 
unrealistic, a literature search generally found enough ac- 
tual equipment to permit a relationship to be empirically 
determined. A linear scaling of such volumes provided 
adequate results. 



196 



Table 44. Listing of off-line weight and volume program 



tinclude ‘stdto.h" 

tdehne K 150.75 /( gear hardness factor M 

It prograi naae: Mt.c, to find the Heights of SNBS groups t/ 



■aini) 

{ 

FILE Ifp, tfopenO; 

double xg, xa, zg, za, pa, ng, na, gra, grg, np, dhp, gn, ds, q, ne, 
h235, h2353, h2354, n2355, h241=0.0, h243, h298=0.0, xprop, 
upc, vpc, ppc, next, vexc, Hbrk, vbrk, pa, v235, 
v2353=0.0, v2354=0.0, v2355, v241, v243; 
extern double sqrtO, powO; 



printf (“VnThis prograt calculates SUBS 200 Heights for screen output.'); 
printf ('\n\ntlhat is the LCG of the propulsion generator(s)? '); 
scanfCXlf, bxg); 

printfCWhat is the LCG of the propulsion aotoris)? '); 
scanfCXlf, ixa); 

printfl'Hhat is the VC6 of the propulsion generator(s)? '); 
scanfCXlf, izg); 

printf (“What is the VC6 of the propulsion aotoris)? '); 
scanfCXlf, !iza); 

printf CWhat is the nuaber of propulsion generator (s)? '); 
scanfCXlf , 4ng); 

printf CWhat is the nuaber of propulsion aotoris)? '); 
scanfCXlf , 4na); 

printfCWhat is the nuaber of gas turbines aboard? ')*, 
scanfCXlf , tne); 

printfCWhat is the rated horsepoHer of each gas turbine? '); 
scanfCXlf, 4pa); 

printfCWhat is the propeller rpa? '); 
scanfCXlf, 4np); 

printf CHoh auch horseponer is delivered to each propeller? '); 
scanfCXlf, 4dhp); 

printfCWhat is the LCG of the propeller? '); 
scanfCXlf , bxprop); 

printfCWhat is the gear ratio at the propulsion aotor? '); 
scanfCXlf , 4gra); 

printfCWhat is the gear ratio at the propulsion generator? '); 
scanfCXlf , 4grg); 



It W23S.3 Transaission lines tl 

h2353 = 0. 0000091 (xa-xg+zg-za+27) I (pat746/30000); It LT, enhanced tl 
v2353 = 0.06545l(xa-xg+zg-za+27); It cubic feet, enhanced tl 



It snag tl 



It W23S.4 Cooling systeas tl 
h2354 = 0.26lpatng/2240i 
v2354 = lOO.OIng; 
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It LT 1/ 



/I K235.5 Switch gear 1/ 

w2355 = 0.261(1)9 + 2ln#l; H enhanced ASSET, LT 1/ 

v2355 = 45.0<w2355; /I switchgear voluee, ft'S, enhanced </ 



ppc = 0.000939523lpa/n«; 
wpc = nit3.5jlppc/12; 
vpc = ntl540tppc/12; 



/I power converter rating, MH, enhanced ASSET 1/ 
/I weight pwr conv, LT, enhanced ASSET 1/ 
/I vol pwr conv, LT, enhanced ASSET 1/ 



pa = pal746/l000000; It aotor rating in HH </ 

wexc = (ne * na)tpowMpi/30) , 0.3); It weight oT exciters, LT, enhanced t/ 

vexc = 70lwexc; It vol of exciters, ft'’3 tl 

wbrk = 0.26tpatna; It weight of braking resistors, LT, enhanced t/ 

vbrk = 37lpatni; It vol of resistors, ft'‘3, enhanced 1/ 



w235 = w2353 + w2354 + w2355 + wpc + wexc + wbrk; 
v235 = v2353 + v2354 + v2355 + vpc + vexc + vbrk; 

It H241 locked-train double reduction gears tl 

if (grg 1= 1.0) It there are pg gears, btw gt and pg tl 

w241 = (1.57tpatpow((grg+l), 3.0)tng)/(3600tgrgtgrgtgrg4K); 
if (gra != 1.0) It there are pa gears, btw pa and propeller tl 

w241 *= (l,S7tdhptpow((gra^l), 3,0)tna)/(nplgratgra(K); It LT tl 

v24l = W241I34.612; It cubic feet, ratioed fa FFG7 tl 

It W243 Shafting tl 

ds = 2.152lpow((4.22347tdhp/np>, 0.333); It shaft diaa, sq inches tl 
w243 = 1.5708ldstdsl2l(xprop - xa - l.0)/2240; It LT tl 

v243 = 0,0109Udsldst(xprop - xa - 6.0); It ft‘'3 tl 



It H29B.1 LTDR operating fluids, additon to wt in ASSET tl 
if (grg != l.O) 

w298 = 0.27lpalng/2240; 
if (qra != 1.0) 

w298 += 0.27tpalni/2240; It LT tl 



printf C\n\nM235.3 Transaission lines 
printf (*\nH235.4 Cooling systeas 
printf ('\nH235. 5 Switchgear 
printf (’\n Power converters 
printf ("\n Exciters 
printf (“\n Braking resistors 
printf (■\n\nW235 Electric propulsion 
printf (’\nH241 Reduction gears 
printfC\n«243 Shafting 
printf (’\nM2?8 Sear operating fluid 
} 



Z7.2H LT X7.2lf fr3*, w2353, v2353); 

Z7.21f LT Z7.21f ft^3’, w2354, v2354); 

X7.2H LT 17.2H 1^3', w2355, v2355); 

Z7.2lf LT X7.2H ft*3‘, wpc, vpc); 

Z7.2H LT I7.21f fr3‘, wexc, vexc); 

X7.2lf LT X7.21f ft*3‘, wbrk, vbrk); 

Z7.21f LT Z7.21f ff'3, less PSs and PHs', w235, v235) 
X7.2lf LT Z7.21f ft‘'3‘, w241, v24l); 

Z7.2H LT Z7.21f ft"3”, w243, v243); 

X7.2H LT\n\n', w298); 
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Appendix F. Advanced Surface Ship Evaluation Tool output 

The output of ASSET is in text and graphic form. The 
total text output for any particular synthesis run is more 
than thirty pages. Following are several graphic outputs of 
ASSET, showing the mechanical and electrical transmission 
ships used in this thesis. 
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Figure 49. Hull isometric view of all ships 
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Figure 50. Body plan of all ships 
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Figure 51. Plan view of subdivision in mechanical baseline 
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Figure 52. Plan view of subdivision in rearranged electrical 
ship 
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